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ABSTRACT 


The  content  of  this  report,  for  the  second  semi-annual  period  of 
contract  support,  comprises  three  main  topics.  The  first  treats  of  the 
relationships  borne  in  glasses  between  microstructure  and  electrical 
properties.  One  study  examines  a  theoretical  model  that  relates  the  aperio- 
dicity  in  glass  lattices  to  electron  transport.  Another  explores  the 
structural  mechanisms  involved  in  the  nucleation  of  crystal  phases  in  glasses. 

The  second  main  topic  deals  with  the  measurement  of  parameters  of 
material  structure,  demonstrating  the  measurement  techniques  by  applying 
them  to  various  materials:  glass  ceramics,  nickel  alloys,  magnetic  films, 
and  degenerate  semiconductors.  Discussion  centers  on  analysis  by  small- 
angle  x-ray  scattering,  by  electron  microprobe,  and  by  transmission  elec¬ 
tron  microscopy.  The  errors  inherent  in  the  measurements  and  the  theoretical 
basis  for  the  analysis  of  data  receive  special  attention. 

The  third  main  topic  concerns  semiconductors  and  semiconductor  devices. 
One  study  reports,  for  degenerate  n-type  germanium  and  silicon,  the  theo¬ 
retical  variation  of  the  density  of  quantum  states  as  a  function  of  impurity 
concentration.  In  another,  the  conventional  theory  of  pn  junctions  is  ex¬ 
tended  by  removing  the  assumption  of  non-degeneracy.  Finally,  for  both 
devices  and  bulk  material,  the  sources  of  noise  are  subjected  to  theoretical 
and  experimental  study. 


SUMMARY 


The  content  of  this  report,  for  the  second  semi-annual  period  of 
contract  support,  comprises  three  main  topics.  The  first  treats  of  the 
relationships  borne  in  glasses  between  microstructure  and  electrical  prop¬ 
erties.  The  second  deals  with  the  measurement  of  parameters  of  material 
structure,  demonstrating  the  utility  of  the  measurement  techniques  by 
applying  them  to  various  materials:  glass  ceramics,  nickel  alloys,  magnetic 
films,  and  degenerate  semiconductors.  The  third  main  topic  concerns  semi¬ 
conductors  and  semiconductor  devices,  with  emphasis  on  degenerate  germanium 
and  silicon  and  on  the  sources  of  noise  in  both  devices  and  bulk  material. 
GLASSES :  To  explain  electron  transport  phenomena  in  glasses  containing  tran¬ 
sition  metal  oxides,  a  theoretical  model  is  developed  that  takes  into  account 
the  lack  of  periodicity  of  the  glassy  lattice  and  the  strong  electron-lattice 
interaction  expected  in  the  transition  metal  oxides.  This  model  is  termed 
the  Nearly  Crystalline  Ionic  Model.  As  a  test  of  the  adequacy  of  this  model, 
comparisons  are  made  between  measured  and  predicted  data  for  AC  and  DC  con¬ 
ductivity.  The  comparison  shows  that  two  special  cases  of  the  model  —  one 
postulating  cells  of  small  disorder  surrounded  by  a  small  region  of  large 
disorder,  the  other  postulating  cells  of  disordered  materials  of  a  high  tran¬ 
sition-metal-oxide  content  surrounded  by  a  small  disordered  phase  of  low 
transition-metal-oxide  content  —  give  results  in  qualitative  agreement  with 
experiment.  The  effects  of  heat  treatments  on  glass  properties  are  attributed 
to  either  phase  separation  or  ordering  of  the  lattice. 

A  second  study  explores  the  structural  mechanisms  involved  in  the  nucle- 
ation  of  crystal  phases  in  a  glass.  It  presents  information  concerning  the 
rates  of  bulk  crystallization  of  glasses  and  concerning  the  relationship 


ii 


between  the  mechanisms  of  nucleation  and  growth.  Binary  L^O-SiC^  and 
Na20-Si02  glasses,  serve  as  the  experimental  vehicles. 

MEASUREMENT  TECHNIQUES  AND  APPLICATIONS:  Of  the  measurement  techniques 
discussed,  analysis  by  small-angle  x-ray  scattering,  by  electron  microprobe, 
and  by  transmission  electron  microscopy  receive  the  primary  attention. 

Special  consideration  is  given  to  the  theoretical  basis  for  the  analysis  of 
data  and  to  the  errors  inherent  in  the  measurements.  This  work  —  coupled 
with  collateral  efforts  concerning  the  measurement  of  electronic  parameters 
of  materials  —  undergirds  a  major  objective  of  our  research  program:  the 
correlation  of  material  structure  and  electronic  properties. 

During  this  second  semi-annual  period,  x-ray  research  has  centered  on 
the  determination  of  particle-size  distribution  by  small-angle  x-ray  scat¬ 
tering.  The  theory  underlying  this  method  is  first  discussed.  Then  the 
method  is  applied  to  dispersion-hardened  nickel  alloys  containing  small 
spherical  particles  of  thoria.  The  x-ray  analysis  yields  results  for  the 
particle  size  distribution  in  good  agreement  with  published  results  deter¬ 
mined  by  transmission  electron  microscopy. 

Moreover,  the  application  of  this  x-ray  technique  to  the  determination 
of  impurity  clustering  in  degenerate  semiconductors  is  considered,  as  is  the 
structural  analysis  of  permalloy  thin  films  by  the  electron  microprobe. 

In  transmission  electron  microscopy,  the  analysis  must  take  account  of 
sample  thickness  to  avoid  introducing  considerable  error  in  the  evaluation 
of  structure.  When  the  ratio  of  sample  thickness  to  particle  diameter  (t/b) 
exceeds  2,  error  may  result.  A  reporting  is  made  concerning  the  magnitude 
of  this  error.  Partially  devitrified  Li20-2Si02  8lass  was  examined  by  both 
transmission  and  reflection  measurements,  and  the  results  thus  yielded  were 
compared  with  results  of  published  x-ray  studies.  This  research  demonstrates 
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that  quantitative  stereological  parameters  can  be  overestimated  by  as 
much  as  300%  when  measured  in  transmission  if  corrections  for  thickness 
effects  are  not  made.  On  the  other  hand,  the  proper  corrections  give  param¬ 
eter  values  that  differ  by  less  than  10%  from  the  values  obtained  from  reflec¬ 
tion  measurements  and  x-ray  analysis. 

SEMICONDUCTORS  AND  SEMICONDUCTOR  DEVICES:  Accounts  are  given  of  two  studies 
concerning  degenerate  semiconductors.  The  first  presents,  for  degenerate 
n-type  silicon  and  germanium,  the  theoretical  variation  of  the  density  of  quantum 
states  as  a  function  of  impurity  concentration.  The  second  study  extends  the 
conventional  theory  of  pn  junctions  by  removing  the  assumption  of  non-degeneracy. 
In  the  analysis,  an  approximation  is  utilized  for  the  Fermi-Dirac  statistics 
and  the  dependence  of  the  quantum  density  of  states  on  impurity  concentration 
is  neglected.  Calculations  using  the  parameters  of  silicon  illustrate  the  pre¬ 
dicted  behavior. 

Research  concerning  noise  is  reported  dealing  with  two  subjects.  One 
concerns  1/f  noise  present  in  planar  bipolar  transistors.  The  other  relates 
to  a  novel  theoretical  formulation  of  recombination-generation  noise  in  semi¬ 
conductors. 
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I.  Introduction 


The  general  objective  of  this  research  program  is  to  establish  at 
the  University  of  Florida  a  center  of  competence  in  solid-state  materials 
and  devices,  with  primary  emphasis  to  be  placed  on  materials  and  devices 
not  used  in  the  conventional  silicon  integrated-circuit  technology. 

From  the  efforts  expended  in  establishing  and  maintaining  this  center 
of  competence  have  evolved  technical  findings:  technical  findings  con¬ 
cerning  such  materials  as  glass  ceramics,  semiconducting  glasses,  magnetic 
films,  and  degenerate  semiconductors;  concerning  devices  made  from  these 
materials;  concerning  measurement  techniques;  and  concerning  methods  of 

fabrication.  The  first  of  the  findings  in  these  various  areas  appeared  in 

1 

Scientific  Report  No.  1  .  To  that  report  the  reader  is  also  referred  for 
a  more  detailed  statement  of  research  objectives  than  given  here  and  for  a 
discussion  of  the  means  to  be  used  in  achieving  these  objectives. 

The  present  report  sets  forth  the  major  findings  of  the  second  semi-annual 
period  of  contract  support.  In  the  presentation  to  follow,  Section  II  des¬ 
cribes  the  results  of  research  concerning  glasses,  Section  III  reports  the 
progress  relating  to  measurement  techniques,  and  Section  IV  treats  of  the 
findings  concerned  with  semiconductors  and  semiconductor  devices. 
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II.  Glasses  (L.  L.  Hench) 


This  sectioa  deals  first  with  a  theoretical  model  relating  the  aperio- 
dicity  in  glass  lattices  to  electron  transport,  and,  then,  with  a  study  of 
the  nucleation  of  the  crystal  phases  of  glass. 

A.  ELECTRON  TRANSPORT  PHENOMENA  IN  TRANSITION  METAL  OXIDE  GLASSES  (H.  F. 
Schaake  and  L.  L.  Hench) 


1 .  INTRODUCTION 

The  phenomenon  of  electronic  conduction  in  glasses  containing  the 
transition  metal  oxides  has  received  considerable  experimental  attention 
in  recent  years.*  Although  a  variety  of  explanations  have  been  advanced 
to  describe  the  electron  transport  process,  essentially  all  the  theoretical 
models  have  been  crystalline.  It  is  the  purpose  of  this  paper  to  develop 
a  noncrystalline  model  which  will  result  in  a  qualitative  understanding  of 
the  transport  process  itself,  and  hopefully,  shed  some  light  on  the  struc¬ 
ture  of  the  glass. 

A  brief  review  of  the  electron  transport  process  in  the  crystalline 

transition  metal  oxides  will  serve  as  a  basis  for  the  remainder  of  this 
3-5 

paper.  In  the  crystalline  state,  each  transition  metal  ion  is  surrounded 
by  six  oxygen  ions  in  octahedral  coordination  (Figure  1).  The  presence  of 
the  oxygen  ions  causes  the  d-electron  energy  levels  of  the  metal  ions  to  split. 
The  wave  function  of  an  electron  in  any  one  of  these  levels  is  highly  localized 
inside  the  octahedron  formed  by  the  oxygen  ions.  There  will  be  a  small  amount 
of  overlap  of  the  wave  function  on  one  metal  ion  with  the  wave  function 
on  a  neighboring  metal  ion,  causing  the  formation  of  a  very  narrow 
band.  Depending  on  the  particular  oxide,  these  bands  will  be  either 


*  For  reviews,  see  Mott  and  Mackenzie 
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completely  filled,  partially  filled,  or  completely  empty.  We  shall 
confine  our  attention  from  hereon  to  the  latter  case.  The  Insertion 
of  donor  states,  either  by  the  addition  of  impurities,  or  as  a  result 
of  nonstoichiometry  caused  by  the  evolution  of  oxygen  during  forming, 
can  result  in  the  presence  of  electrons  in  the  normally  empty  band. 
Consequently  the  oxide  will  be  a  semiconductor. 

These  excess  electrons  are  usually  thought  of  as  being  associ¬ 
ated  with  a  certain  ion;  the  electron  concentration  being  specified 
by  the  Me^^/Me^11*^*  ratio.  Although  the  motion  of  the  electron 
through  the  lattice  is  via  band  diffusion,  because  the  wave  function 
of  the  electron  in  the  lattice  remains  highly  localized  about  the 
metal  ions,  the  electron  transport  process  is  frequently  described 
as  a  hopping  of  the  electron  from  ion  to  ion. 

The  narrowness  of  the  conduction  band  implies  that  the  velocity 
with  which  the  electron  moves  through  the  lattice  will  be  extremely 
small.  As  a  result  of  the  low  velocity  and  the  highly  polar  nature 
of  the  lattice,  the  lattice  will  distort  due  to  the  field  of  the 
electron.  The  combination  of  lattice  distortion  and  electron  is 
treated  theoretically  as  a  quasi-particle  known  as  a  polaron,  and 
behaves  similarly  to  the  "undressed"  electron,  but  with  a  higher 
effective  mass,  and  with  a  polaron  bandwidth  narrower  than  the  pure 
eleotron  bandwidth.  Due  to  the  strong  coupling  with  the  lattice, 
the  polaron  will  be  greatly  affected  by  lattice  vibrations  (or  phonons) , 
and  as  the  temperature  increases,  an  increasing  portion  of  the  total 
conduction  will  be  due  to  a  thermally  activated  hopping  process. 

If  the  velocity  of  the  electron  is  slow  enough,  and  the 
electron-lattice  coupling  large  enough,  the  lattice  in  the  immediate 
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vicinity  of  the  electron  will  be  able  to  react  to  the  field  of  the 

electron,  effectively  trapping  it.  The  resultant  combination  of  the 

electron  and  lattice  distortion  is  known  as  a  small  polaron.  At  low 

temperatures  the  small  polaron  will  move  by  band  diffusion,  although 

with  a  very  high  effective  mass ,  and  an  extremely  small  polaron  band- 

-2  -A 

width  (of  the  order  of  10  to  10  that  of  the  pure  electron  bandwidth) . 
As  the  temperature  increases,  the  effective  mass  of  the  polaron 
increases,  and  the  bandwidth  decreases.  At  sufficiently  high  tempera¬ 
ture  (of  the  order  of  the  Debye  temperature),  however,  transport  is 
predominantly  by  a  thermally  activated  hopping  process,  in  which 
phonons  are  emitted  and  absorbed  as  the  small  polaron  moves  from  site 
to  site. 

There  is  no  agreement  at  this  time  as  to  which  mechanism  is 
present  in  each  of  the  crytalline  transition  metal  oxides,  although 
recent  experimental  evidence  for  NiO  favors  the  polaron,  rather  than 
the  small  polaron. 1 

Turning  now  to  the  situation  in  the  glassy  state,  the  theo¬ 
retical  picture  is  much  less  clear.  Mott1  has  recently  reviewed  the 
current  theory  of  electrons  in  a  disordered  lattice.  There  appears 
to  be  general  agreement  that,  if  the  lattice  is  sufficiently  disordered, 
and  if  the  binding  of  the  electrons  to  the  lattice  ions  is  sufficiently 
great,  then  all  electron  states  in  the  solid  will  be  localized. 
Localization  means  that  each  wave  function  will  have  a  maximum  about 
some  site  in  the  lattice,  and  that  the  envelope  of  the  wave  function 
on  the  average  will  decay  exponentially  with  distance  from  that  site. 
Conduction  by  band  diffusion  is  impossible;  thermally  activated  hopping 
with  the  emission  and/or  absorption  of  phonons  is  then  required  for 


-5- 


the  electron  to  diffuse  through  the  lattice.  We  note  that  the  tran¬ 
sition  metal  oxide  containing  glasses,  as  a  result  of  the  Inherent 
tight  binding  of  the  d-electrons  may  be  expected  to  exhibit  localized 
states. 

A  rigorous  treatment  of  the  electron  transport  process  in 
a  three-dimensional  glass  is  completely  impossible  at  the  present 
time.  This  is  a  result  of  the  following: 

(1)  The  lack  of  a  detailed  model  for  the  glassy  state.  There 
are  two  conflicting  models  of  the  glassy  state  at  the  present  time:5 
the  microcrystalline  model  of  the  Russian  School,  and  the  random  model 
of  Zachariesen,  which  is  generally  held  in  the  West.  One  point  of 
agreement  by  both  schools  of  thought  is  that  the  anion-cation  complex 
present  in  the  crystalline  state  (the  metal  ion  with  oxygen  ions  in 
octahedral  coordination  for  the  transition  metal  oxides)  is  maintained 
in  at  most  a  slightly  distorted  state.  This  point  can  be  inferred 
from  x-ray  diffraction  data,6  and  is  supported  by  infrared  studies.7'8 

(2)  The  construction  of  electronic  wave  functions  for  a  three- 
dimensional  disordered  lattice  is  a  difficult  problem,  particularly 

in  the  absence  of  a  specific  model  for  a  glass.  Only  recently  have 
some  papers  appeared  in  an  attempt  to  construct  such  wave  functions.8 

(3)  There  has  been  little  success  in  constructing  phonon 
wave  functions  in  a  disordered  three-dimensional  lattice.10  Classical 
calculations  do  indicate  that  high  frequency  modes  are  localized  in 
specific  regions  in  the  glass.7’8 

We  shall  therefore  limit  our  attention  to  a  one-dimensional 
one-electron  model.  Here  the  tight  binding  model,  which  assumes  that 
wave  functions  in  the  solid  may  be  approximated  by  a  linear  combination 


of  the  wave  functions  in  the  isolated  atoms  (or  metal-oxygen  complexes 
in  our  case) .  We  shall  furthermore  use  as  a  framework  for  our  calcu¬ 
lations,  Holstein’s  small  polaron  approach.5  The  justification  for 
this  step  is  twofold: 

(1)  Even  though  the  small  polaron  picture  may  not  be  appli¬ 
cable  to  the  crystalline  state,  if  there  is  sufficient  disorder,  it 
may  be  applicable  to  the  glassy  state.  In  the  glassy  state,  the  wave 
functions  will  be  localized  if  there  is  sufficient  disorder.  Locali¬ 
zation  implies  that  the  lattice  outside  the  region  where  the  electron 
is  localized  will  be  highly  polarized,  since  that  portion  of  the 
lattice  will  always  see  the  field  of  the  electron.  Furthermore,  the 
lattice  inside  the  region  of  the  localized  electron  will  be  somewhat 
polarized  as  a  result  of  the  exponential  decay  of  the  electron  wave 
function.  If  the  situation  in  the  crystalline  state  is  marginal; 
that  is,  near  to  assuming  the  small  polaron  configuration,  then  with 
sufficient  disorder  it  may  in  fact  assume  the  small  polaron  configuration. 
This  possibility  is  further  enhanced  by  the  greater  average  separation 

of  the  metal  ions  in  the  glass  in  comparison  to  the  crystalline  state, 
and  the  presence  of  other  species,  resulting  in  decreased  overlap 
of  adjacent  wave  functions. 

(2)  Holstein's  theory  can  be  qualitatively  extended  to  the 
normal  polaron  problem. 

In  Chapter  II  we  shall  therefore  cast  Holstein's  model  into 
a  form  suitable  for  an  ionic  lattice,  and  take  into  consideration 
the  lack  of  periodicity.  In  Chapter  III  we  shall  apply  the  results  of 
our  calculations  to  the  observed  AC  and  DC  conductivity  data  recently 
obtained  on  potassium  vanadium  phosphate  glasses.  Some  additional 
implications  of  the  model  shall  also  be  discussed. 
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2.  THEORETICAL 


Nearly  Crystalline  Ionic  Model 


In  accordance  with  our  previous  discussion  we  shall  develop 
the  following  model,  which  we  shall  refer  to  as  the  Nearly-Crystalline 
Ionic  Model.  Consider  a  one-dimensional  non-periodic  lattice  consisting 
of  alternating  metal  and  oxygen  ions  (Figure  2).  We  shall  assume  that 
the  forces  acting  on  the  ions  are  harmonic.  In  place  of  the  unit  cell 
which  would  concern  us  if  we  were  discussing  a  crystalline  model,  we 
shall  use  the  concept  of  a  metal-oxide  complex,  consisting  of  a  metal 
ion  and  the  two  adjacent  oxygen  ions.  We  shall  consider  our  chain  to 
have  N  complexes  or  2N  atoms.  The  analogy  with  a  three-dimensional 
complex  consisting  of  a  metal  ion  and  6  oxygen  ions  in  approximately 
octahedral  coordination  should  be  obvious. 

In  the  following  discussion  we  shall  use  the  following  notation: 
x  will  denote  the  electron  coordinate,  u^  the  displacement  of  the 
i^  ion  from  its  equilibrium  position  in  the  n*"*1  complex.  We  shall  use 
the  notation  ur  to  denote  the  triplet  of  displacements  (un_i»  ung»  un]^ 
and  the  notation  u  to  denote  the  entire  set  of  displacements  of  the  2N 
atoms.  Our  other  notation  will  follow  closely  that  used  by  Holstein/ 

The  Schrbdinger  equation  for  one  electron  in  the  lattice  will  be: 


H^ 


-tT 

i  3t 


(1) 
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Figure  2.  Nearly  Crystalline  Ionic  Lattice. 


where  the  Hamiltonian  will  be  given  by: 


H  ■  H  +  H,  .  +  H,  , 

electron  lattice  int 


The  electronic  portion  of  the  Hamiltonian  is 


electron 


-Is 72  +  £  "i  u„i> 

3x  n,i 


where  (x-R^,  Uni)  denotes  the  potential  of  the  iL  atom  of  the 
n^  complex  which  is  located  at  R  ^  from  the  origin. 


The  lattice  vibrational  portion  of  the  Hamiltonian  will  be 


given  by: 


lattice 


I'  - 
n,i 


2  2  1  \ 

It  3  .  /nil 

“i  3u2 .  8  M  V 

ni 


In  this  expression,  g  I  ^ 1  is  the  force  constant  acting  between  the 

th  th  *  ^  ^ 

ifc 1  atom  of  the  n  complex  and  the  j  atom  of  the  H  complex.  The 


primed  sum  indicates  that  the  summation  goes  over  each  atom  once;  it 
is  required  due  to  the  equivalence  of  u^  ^  and  un  ^ 

The  interaction  portion  of  the  Hamiltonian  H^nt  describes  the 
interaction  of  the  electron  with  the  lattice  vibrations  and  is  contained 
in  Heiectron  due  to  the  explicit  dependence  of  the  on  the  displacements. 

We  now  wish  to  solve  equation  (1).  To  do  so,  we  shall  use  the 
tight  binding  approximation  which  is  appropriate  for  the  narrow  bands 
of  the  crystalline  transition  metal  oxides.  We  consider  first  the 
n^  isolated  complex.  The  electronic  wave  functions  will  satisfy  the 


following  SchrBdinger  equation: 


2  2 

SliL+  z 

2m  ,  2  . 

3x  i 


U1  uni\  *n 


E  ( u  )  <p 
n  n  Yn 
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Both  (J>n  and  En(un)  will  depend  on  the  displacements.  We  have  chosen 
to  explicitly  display  this  for  the  En(un).  The  wave  function  for 
the  entire  solid  will  be  given  by: 


ij.'  ■  2  a(u)  4>  (4) 

m  m 
m 

Inserting  equation  (4)  into  equation  (1)  and  making  use  of  equation  (3), 
we  find: 


-ft  3a„(u) 

_11  z  — E -  * 

i  3t  vm 
m 


=  Z  E(u) a  ( u) $  +  Z' 

n  n  n  n  . 
n  n,i 


2  2 
IT  3^ 


2M^  a  2 
i  3u  . 

ni 


x  a  (u)  + 
m 


n,i 


uni  utj  am<u) 


m 


+  Z' 
n,i 

m 


-n 

2M, 


3  2 


?  3am(u)  , 

2  3u  J  3u  ,  am(u)  „  2 


m 


ni  ni 


3u 


ni 


+  Z'  U.  (x-R  , ,  U  ,) 
.  i  ni’  ni 
n,i 

m^n 


x  *m  am(u)  (5) 

We  now  multiply  both  sides  of  equation  (5)  by  <p*  and  integrate  over 

P 

all  space  with  respect  to  the  electronic  coordinate  x.  In  doing  so, 
we  shall  come  across  integrals  of  the  form: 


T 

mp 


d>  dx 
m 


(6) 


Due  to  the  localized  nature  of  the  <j>  wave  functions  which  allows  us 
to  use  the  tight  binding  theory,  we  shall  follow  Holstein5  and  approx- 
image  equation  (6)  by: 


T  =6 
mp  mp 


(7) 
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Equation  (5)  becomes: 


3a„(u) 
JL _ P 

i  at 


v  vvu)  +  „e; 

n,i 

P.j 


si.  3  S(u)  .  Ini) 
2Mi  3a2,  8pJ. 


x  u  .u  .a  (u) 
ni  pj  P 


+  Z* 
n,i 
m 


It 

2M, 


2 /*i 


3a  (u) 


m 


p  3u 


dx 


ni 


a%  ' 

-  '♦}  ~f  dx 

3uni 


+  ui  ^X-Rni*  Uni) 

n,i 

m^n 


x  K  dx  am  <«) 


(8) 


Consider  the  last  sum  in  equation  (8).  Because  of  the  localized  nature 
of  the  three  functions  under  the  integral  sign,  we  can  approximate  the 
integral  as  zero  unless  two  functions  refer  to  the  same  site,  and  the 
third  refers  to  a  neighboring  site.  Therefore,  this  sum  becomes: 


Z  6 

e+1 


5  .  f <f>*  I  U.  (x-R  .  ,  u  .)  4  dx  a 

m,p+e  p  i  mi  mi  Tp  p 


(u) 


+  Z  5  ftp*  Z  U  (x-R  ,  u  .)  <p  dx  a  (u)  (9) 

e+1  P.m+e  pi  1  pi  pi  ym  mv 

Next,  consider  the  second  sum  of  the  right  hand  side  of  equation  (8), 
which  contains  derivatives  of  the  wave  functions  with  respect  to  the 
nuclear  displacements.  These  terms  should  be  small  (i.e.,  the  displace¬ 
ments  affect  the  wave  functions  only  slightly) ,  and  we  shall  therefore 
follow  Holstein5  and  set  this  entire  term  to  zero.  With  these  approxi¬ 
mations,  equation  (8)  becomes: 
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E ' 
F.i 


fi  _3_ 
i  3t 


2M 


+  8 


i  3u 


m 


pi 

11 


u  ,  u„ .  -  E  (u  ) 
pi  ?.j  n  n 


an(u) 


e=+l 


6  J  (u) 

n,itH-e  mn 


a  (u) 
m 


(10) 


where : 


Vu> 


'  /**  t  u 

-+1  “  i 


u«i)( 


dx 


(ID 


J 

mn 


(u) 


(.* 

n 


U.  (x-R 


ni’ 


Uni>' 


m 


dx 


(12) 


Inspection  of  equations  (11)  and  (12)  shows  that  WR(U)  represents  the 

perturbation  on  the  energy  of  the  isolated  complex  En(un)  due  to  the 

neighboring  complexes  while  >lmn(u)  is  an  overlap  integral  which  will 

allow  for  transitions  between  the  two  adjacent  sites. 

Equation  (10)  is  similar  to  that  derived  by  Holstein.  There 

are  two  important  differences  ,  however.  First  we  see  that  the  terms 

E  (u  )  +  W  (u)  and  J  (u)  are  in  general  different  for  each  site, 
n  n'  n  mn  ° 

due  to  the  effects  of  disorder.  Secondly  the  lattice  vibration  inter¬ 
action  is  written  in  a  more  general  form. 

A  general  solution  to  equation  (10)  is  still  impossible  without 
further  approximations.  We  shall  therefore  make  the  following 
simplifications : 


Mx  =  M2  =  M 

=  g  (constant,  between  nearest 
neighbors  only) 


(13a) 

(13b) 
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E(u)  +  W  (u) 
n  n  n 


A-B  u  A  -  C  (u  .  -  u  n) 
n  n  n,0  n  n,l  n-1,1 


2  3 

-  C  (u  ..  _  -  u  ,  n)  -  C  (u  , ,  ,  -  u  „  ..) 
n  n+1,0  n-1,0  n  n+1,1  n-2,1 


A  -Bu  .  •  I  C  (u.  - 
n  n  n,0  n  £  6(£) 


(13c) 


J  (u)  ■  J  (constant)  (13d) 

run  run 

Approximations  (13a)  and  (13b)  allow  us  to  establish  analytic 
normal  coordinates  for  lattice  vibrations.  The  problem  reduces  to  2N 
equal  masses  interconnected  by  springs  of  force  constant  g.  In  an 
actual  glass  the  nearest  neighbor  force  constants  should  be  slightly 
different  due  to  anharraonic  terms  in  the  lattice  potential.  This  would 
reduce  the  problem  to  a  nonanalytic  problem  in  which  the  short  wave¬ 
length  modes  become  somewhat  localized.  By  assuming  the  force  constants 
are  all  equal,  we  shall  therefore  be  assuming  a  strictly  harmonic 
approximation.  Assuming  the  masses  equal  restricts  our  phonon  spectrum 
(Figure  3);  however,  in  consideration  of  the  large  differences  between 
one-  and  three-dimensional  spectra,  and  in  light  of  the  success  of  the 
Debye  approximation,  which  is  similar  in  its  effects  as  our  approxima¬ 
tion,  we  expect  qualitative  agreement  between  our  model  and  reality. 

The  approximation  of  the  energies  as  being  linearly  related  to 
the  displacements  of  the  ions  appears  to  be  standard  in  polaron  theory. 
It  allows  the  problem  to  be  recast  into  displaced  normal  coordinates, 
as  we  shall  see.  Inclusion  of  higher  order  terms  would  recouple  the 
phonon  modes  which  the  normal  coordinates  decouple,  and  hence  shift 
the  phonon  frequencies.  The  approximation  (13c)  therefore  corresponds 
to  a  series  expansion  of  the  energy  about  zero  displacement  and  the 
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ignoring  of  nonlinear  terms.  The  term  An  reflects  that  different  sites 

have  different  energy.  The  Bn  term  yields  the  shift  of  energy  due 

to  the  displacement  of  the  metal  ion  of  the  site  on  which  the  electron 

is  located.  The  yields  the  effect  of  displacement  of  nearest 

2 

neighbor  oxygen  ions,  Cn  of  nearest  neighbor  metal  ions,  and  so  on. 

We  shall  see  that  the  Bn  term  yields  essentially  meaningless  results, 
but  we  include  it  as  it  is  similar  to  Holstein's  approximation. 5 

Finally,  the  approximation  (13d)  states  that  we  shall  ignore 
the  effect  of  displacements  on  the  overlap  integral,  and  is  justified 
in  view  of  the  smallness  of  J 

mn 

With  the  approximations  of  equation  (13),  equation  (10)  may 
be  rewritten. 


-h  _3_ 

i  at 


E ' 
P,i 


Vi-l5' 


A 


n 


B  u  .-EC  (u.  -  ux/n.) 
n  n,0  ^  n  £  6(£) 


E  6  . 

E-fl  n»mfG 


x  J  a  (u) 
nm  m 


(14) 


We  transform  to  normal  coordinates  by  setting: 

\l/2 


ni 


q^  sin  (2n+i)k  + 


(15) 


In  doing  so,  we  find  the  frequencies  given  by: 


2  2  2  1 
-  u)Q  (sin  ^  k); 


2# 

VM 


(16a) 


with : 


k  ■  where  £  is  an  integer  in  the  interval  (16b) 

N  -2N  <  £  <  2N 
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Substitution  of  equation  (15)  into  equation  (14)  yields 


["  i  h  "  l  (  '  2M  +  \  M  “k  qk  )"  An  +  Bn  1  (N}  f  qk  sin(2nk  +  f>) 


+  c„  £  <|)1/2  \[  sin[(2n+l)k  +  -J]  -  sin[(2n-l)k  + 


+  C2  Z  (1}1/2  qk[  sin[(2n+2)k  +  -J]  -  sin[(2n-2)k  +  -J]j 
Ic  N 


+  •  • 


•]s 


(q)  -  Z  6  J  a  (q) 

t  ^  n,mfe  nm  m  n 


(17) 


Collecting  terms  involving  q^,  and  completing  the  square  to  obtain  a 
term  of  the  form  (q^  -  q£n^)  we  find: 


where : 


L  J  6  .  a  (q) 

e=+l  mn  m  M 


(18) 


,(n) 


i  1/2  i  r 

=  u  2  s^n(2nk  +  -^)  +  Z  Cn  sin  2£k 

L  o  n 


and: 


x  cos(2nk  + 


*3 


(19) 


i?  = 


A  +  E  — ^  B  sin(2nk  +  7)  +  Z  C£  sin  2£k 

b  n  k  2NMu2  L  n  4  £  n 


x  cos(2nk  + 


*] 


(20) 
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Equation  (20)  represents  the  polaron  binding  energy.  We  now  treat  J 

mn 

as  a  small  perturbation.  To  obtain  our  first  order  results,  we  set 


J  “0.  We  obtain  the  equation: 
mn  n 


?  avP(d) 

i  3dt 


2  2  2 
r  TT  3  1  M  2  ,  (p)/  p.  . 

I  2M  2  +  2  M  “k  <qk  *  qk  }  x  <q) 

kL  3qk  I 


Equation  (21)  is  the  equation  for  a  set  of  decoupled  displaced  harmonic 
oscillators.  The  solution  is : 


Mw.  1/2  .  \ 

XP(q)  -  \  *N(  (qk  -  qk(P))l  exp  - 


where : 


%  "  J  h  “k  <\  +  2> 

k 


♦  I,  ♦  E<» 


n  i / o  1/2  /j 

>N(s)  -  (2  N!  TT±/Z)  e~“  U  Hn(z) 


(23a) 


(23b) 


In  equations  (23) ,  is  the  occupation  number  of  the  k  mode  with 

frequency  and  the  H^z)  are  the  Hermite  polynomials. 

We  now  "turn  on"  the  perturbation  J  .  We  look  for  solutions 

mn 

of  the  form: 


an(q)  -  I  C(p,N)  a  (n,q) 
p,N  P>N 


with 


N 

ap  ,N^n,q3  '  6np  *  (q)  exp  -  — 


(24a) 


Inserting  equation  (24)  into  equation  (18)  we  find: 


.  5  v  iMi  a  (n  q)  exp _ 

i  1  3dt  p,NC  ,q;  P  fi 


L  6  J 

,  ,  n,mfe  mn 
e=+l 


x  C(p' ,N ' )  a  ,  , (m,q)  exp  - 


ltEN’ 
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Multiplying  both  sides  of  equation  (25)  by  the  complex  conjugate  of 


equation  (24a),  and  integrating  with  respect  to  all  q,  ; 


~  7  -C-^-N'  =  z  ( pN  |  V  | p * N * )  C(p\N’)  exp 

1  3t  P'.N’ 


it. 

■ft 


*  <eS-en’> 


(26) 


with: 


<P | v|p'N' )  -  ^  Jpp,  5pipl+£  ,fc  /  ^ 


-h 


L\ 


1/2 


v 

(M\ 

ft 

k 

L\ 

h  -  qkP) : 


To  evaluate  equation  (27),  we  expand  4>, 


1/2 


N, 


ft 


1/2 


d9i 


(27) 


(qk  -  qkp  h 


power  series  about  q£P^  -  q£P  ^  =  0.  We  find: 


in  a 


N, 


Majk  \l/2 
~ft 


1  -  q<P ' } 

qk 
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ii 

,-n 

J 

1/2 


<\  -  4P)> 


9*N  (z) 
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9z 


I  Mas. 
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1  * 


1/2 


<4p,)  -  <p> 


3  *N  <z> 
+  1  — JL- 
2  .  2 
9z 


-  4P)>2  * 


Mu>, 


z= 


f"“k  1/2(0  >>)  1 

hr  <qk  -  qk  ) 


(q 


(p’) 


(28) 


and  making  use  of  the  relation: 


9*  (z) 

K 

9z 


Nk|1/2 


<z>  ’ 

k-1 


'N, +1 
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1/2 


%  (Z) 
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(28a) 
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Equation  (28)  becomes: 


/Mwjtl/2  /_.n 

V,  hr  > 
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Mu,  1 1/2 


'  -  <> 
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2  wk  +  2>  \  hr)  (qk  -  Ik  >  l~)(qk  - 


/  \  2 

')  +  higher  order  terms 


(28b) 


Evaluating  the  integral  appearing  on  the  right  hand  side  of  equation 
(27),  and  using  the  orthonormality  relation: 


f  ^  (z)  <t>N.(z)  dz 

-°°  k  k 


Nk*Nk 


(28c) 


We  find: 
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We  can  evaluate  (q^  -  from  equation  (19): 


<4P' 


’’  -4P)>  -  (s^tMv  si"<2','k  +  l) 


+  I  C  ,  sin  2£k  cos  (2p'k  + 


*] 


-  ~2  [Bp  sin(2pk  +  |) 

Mai.  *- 


+  Z  C£  sin  2£k  cos(2pk  +  -j) 

£  P  4  -I 


(30) 


To  proceed  further,  we  need  to  place 

B  ,  =  B  -  B 

P  P 

c£,  =  c£  =  c£ 

p  p 

This  assumption  is  equivalent  to  saying  that  the  lattice-electron 
coupling  is  the  same  in  all  complexes.  While  not  exactly  true,  we 
expect  the  coupling  to  be  nearly  the  same  in  all  complexes.  Using 
equations  (31),  we  can  rewrite  equation  (30): 

a/2 


(31) 


,  ( P ' )  (pK  /  1\  '  2  sin  ke  r,  r,,  ,  ,  „  .  tt^ 

(qk  }  =  In)  ~72 — LB  cos  [(k(p+e)  +  4} 


Mi 


\ 


-  Z  C  sin  2£k  sin  [k(p+e)  + 

£ 


i»] 


(32) 


where  p  -  p  +  e 


In  the  following  discussion  we  shall  deal  with  two  cases: 
Case  I  B  =  0 

|2 


sin^  k 
-ho)k 


£  "1' 

Z  C  sin  2£k 

-  £ 


(33a) 
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Case  II 


(I  C  sin  2£k)  =  0 

SL  £ 

2  2 
B  sin  k 

MajkK 


(33b) 


We  make  use  of  equations  (33),  (32)  and  (29)  and  rewrite  equation 
(27)  as: 


( pN | V | p '  N ' ) 


Z 

e=+l 


J  ,  6  , 

PP  P.P+ 


+ 


l| 

2  JYk 


x  sin^[k(p+e)  +  -|] 


1 2  Nj'. 

VNi  + 


1/2 


.1/2 


x  sin[k(p+e)  +  |]  6  , 

k’  k+1 


2(Nk+l) 


1/2 


*,k1/2  sln[k(p+c)  +  fj  6  j  (34) 

for  Case  I.  For  Case  II  the  sine  terms  are  replaced  by  cosines  of 
the  same  argument. 


Solution  for  Crystalline  Lattice 

In  anticipation  of  difficulties  which  we  shall  encounter  when 
we  treat  the  disordered  case,  we  shall  first  treat  the  crystalline 
lattice  first.  In  doing  so  we  find  that  approximations  (13b)  and 
(31)  become  exact,  and  that  all  of  the  An's  of  equation  (13c)  become 
equation,  making  all  the  E^^'s  of  equations  (20)  and  (23a)  equal. 

We  restrict  our  attention  to  the  case  when  all  N^'s  remain  unchanged. 
Using  equation  (34)  ,  equation  (26)  becomes 
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H  3C(p, 
i  3t 


Z  J  ,6  ,  ,  7i.  1  -  U  N.  +  ^  y, 

e=+l  PP  P»P  +c  k  [  lN/l  k  2  k 


X  sin2[k(p+e)  +  -|]  C(p' ,N) 


We  note  that  each  of  the  terms  inside  t'.ie  product  bracket  will  be 
nearly  unity,  since  the  N  term  is  much  greater  than  N^.  Hence  we 
can  approximate 


1  -  a  =  e 


(36a) 


"k  (1  ‘  ak> 


(36b) 


Equation  (35)  therefore  becomes 


i  3t 


Z  J  ,6  ,,  e  S(T)  C(p',N)  (37) 

£=+1  PP  P,P  +e 


where 


S  (T)  =  Z  — ^  sin2  [k(p+e)  +  |]  (Nfc  +  |) 

k 


Since  for  every  k  there  is  a  -k 


S<T)  -  I;r<Nk  +  i) 
k 


(38a) 


We  can  now  replace  by  (exp  etlu^  -  1)  ,  where  6  =  1/kT.  Equation 

(38)  becomes 


y  enu. 

S(T)  =  2  "2N  coth  ~ 

k 


(38b) 
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itE 

Assuming  a  time  dependence  of  exp  (— — )  of  the  C(p,N),  equation  (37) 
becomes 

E  C(p,N)  =  Je"S(T)  C(p'+1,N)  +  C(p-1,N)  (39) 

where  we  have  made  use  of  the  fact  that  J  ,  =  J  ,  =  J  for  the 

P.P+1  P.P-1 

crystalline  case.  The  solution  to  equation  (39)  will  be 

C(p,N)  =  ei0p  (40) 

because  of  the  periodicity.  The  energy  will  be  given  by: 

E  *  2J  cos  a  exp  [— S (T) ]  (40a) 

The  C(p,N)  assures  that  our  solutions  are  of  the  required  Bloch  form. 


Solution  for  Disordered  Lattices 


We  are  now  against  one  of  the  most  formidable  problems  in 
contempory  theoretical  physics,  and  one  for  which  there  is  not  a 
complete  methematical  solution  as  of  now.  We  are  therefore  forced 
to  anticipate  the  results  of  a  theoretical  solution. 

We  can  prove  that  in  our  problem  localized  states  may  be 
formed.  This  can  be  seen  by  writing  equation  (26)  for  the  case 
when  N  =  N' : 


H  i_ 

i  8t 


C(p,N) 


exp 


EP  C(p,N)  exp 


+  Z  (pN|v|p'N)  C(p'N)  exp 
P* 


(26a) 
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Therefore,  by 


By  taking  N  =  N'  ,  I  (N^  +  =  1  (N^  + 

k  k 

equation  (23a)  and  equation  (20),  equation  (26a)  becomes 


ti  _3_ 

i  at 


C(p,N) 


exp 


Ap  C(p,N) 


exp 


i  tA  \ 

+  l  (pN | V  j  p ’ N)  C(p'N)  exp  -  — ^  (26b) 

p-  71  I 

Equation  (26b)  is  identical  to  that  of  Anderson, : 1  and  we  may  use 

<A>  S(T) 

the  results  of  his  theorem,  that  if  —j—  e  is  sufficiently  large, 
all  states  will  be  localized  (in  the  sense  that  nonthermally  activated 
diffusion  is  absent).  Here  <A>  represents  the  deviation  of  the  An> 

The  problem  of  solving  for  these  localized  states  is  consider¬ 
ably  more  complicated,  however.  To  circumvent  this  problem,  we  note 
that  in  the  development  up  to  equation  (35),  no  use  was  made  of  the 
energy  differences  on  different  sites,  other  than  to  carry  it  in  the 
time  development  exponential.  Furthermore,  we  note  ;hat  until  this 
point,  the  local  site  variable  p  appears  only  through  the  electron- 
phonon  coupling  term,  in  the  displaced  normal  coordinate  q£^ . 

We  therefore  propose  to  handle  the  disorder  in  the  following 

manner : 

1.  We  begin  by  taking  the  crystalline  lattice  as  our  first 

approximation,  that  is,  the  J's  are  all  taken  equal,  and 

the  A  's  all  equal, 
n  n 

2.  At  this  point  we  introduce  disorder  into  the  problem. 

The  energy  difference  in  the  An's  will  now  appear  in  the 
solution  for  our  C(p,N)'s.  We  can  also  allow  our  J's  to 
be  different;  the  differences  will  affect  our  solution 
for  C(p,N) 
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3.  At  this  point  we  can  follow  Mott,1  and  take  as  the  solution 


for  C(p,N)  to  be: 

-Y  |S  -S 

Cr(p,N)  =  F(r,p)  e  r'  P  r 


(41) 


that  is,  a  state  localized  on  site  r  located  at  S^.,  The 
justification  for  this  assumption  has  been  summarized 
by  Mott: 

a.  The  solutions  to  the  random  Kronig-Penney  model  and 
the  random  6-well  model  are  of  this  type  provided 
there  is  sufficient  disorder . 1 2  * 1 ? 

b.  The  wave  functions  of  impurity  atoms  in  semiconductors 
are  of  this  type.  The  theory  of  impurity  semicon¬ 
ductors  shows  that  the  wave  functions  of  the  impurity 
atoms  are  localized  due  to  random  fields.  14  Further¬ 
more,  the  impurity  conduction  is  found  to  fit 
Anderson's  criteria.  We  should  point  out  that  our 
current  problem  is  akin  to  impurity  conduction 
experimentally . 

We  shall  anticipate  testing  the  validity  of  this  assumption 
by  comparison  of  our  theory  with  experimental  results. 

Inspection  of  the  results  of  the  disordered  Kronig-Penney  and 
disordered  6-well  models  shows  that  F(r,p)  will  be  highly  sensitive 
to  local  structure;  the  important  point  is  that  the  envelope  of  Cr(p,N) 
will  fall  off  exponentially  on  either  side  of  the  localized  center. 

The  physical  meaning  of  equation  (41)  is  that  in  the  disordered 
lattice  the  small  polaron  will  not  be  confined  to  one  site  Even 
though  adjacent  sites  have  different  energy,  the  polaron  will  tunnel 
to  adjacent  sites,  or  equivalently,  will  "resonate,"  in  the  sense 
used  by  Pauling, with  them. 
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As  in  the  crystalline  case,  the  C  (p,N)  will  be  defined  only 

on  a  site.  Also,  since  equation  (41)  will  represent  the  steady  state 

£ 

solution  of  the  problem,  the  C  (p,N)  must  be  orthonormal;  i.e., 


Z  F*(r'p)  exp 
P 


F(r,p)  exp 


-Y, 


=  6 


rr 


(42) 


To  each  localized  state  there  will  be  an  energy  Vr  which  is 
the  counterpart  of  E  in  equation  (40a)  for  the  crystalline  case.  V 
will  now  contain  the  energy  shifts  due  to  disorder. 

Now  we  wish  to  calculate  the  transition  rate  for  a  polaron 
in  a  localized  state  at  r  to  jump  to  a  localized  state  at  site  r' , 
This  problem  is  not  equivalent  to  a  calcuation  for  a  jump  from  a 
polaron  on  site  r  to  a  polaron  on  site  r' ,  because  the  eigensolutions 
to  our  problem  are  spread  out  due  to  tunneling. 

The  total  wave  function  for  a  polaron  in  a  localized  state 
on  site  r  will  be  given  by 


anr(N) 


Z  F(r,p)  exp  ”Yr  S 
P 


S 

r 


6 

np 


1/2 


Ui,  -  q 


(p). 


it 

exp 

~  -h 

A  + 


(N, 


hu). 


+  E,  +  V 
b  r 


(43) 
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Also  of  interest  is  the  fact  that 


/  E  a*r(N)  ar  (N')dq  = 
n  n 

n 

E  F*( r ,p)  F(r ' ,p) 
P 

exp  j 

x  Is  -  S  , 

|  P  r' 

exp  j 

1  \  ( 
[_vr!Sp  -  Sc  j 

r 

\ 

•h 

-i  1 

r 

- 

1/2 


-  <’k  -  4P)) 


Mai,  \  1/2 


S) 


(’k  -  <1kP>) 


dqk  exp 


it 

■h 


,  iyr  -  Vr,  -  I  ANkb„k| 


^  =  6  i  6.... , 

I  rr  NN 


(44) 


by  virtue  of  equations  (42)  and  (28c).  In  other  words,  a  polaron 
cannot  hop  from  one  localized  state  to  another  without  hopping  from 
one  site  to  another  at  the  same  time. 

To  return  to  the  problem,  we  wish  to  computer  the  transition 
rate  from  a  state  localized  at  r  to  one  localized  at  r' .  We  take 
as  our  wave  functions : 


an(q)  *  I  Q  (r',N')  a^'(N')  (45) 

n  r.#N.  n  n 

Substituting  equation  (45)  into  equation  (17),  with  a  reminder 
that  all  A^'s  are  equal,  multiplying  both  sides  by  a*r(N) ,  integrating 
over  all  q^,  we  find: 

3Q  (rN) 

-T — -  =  1  (nN|v|n'N')  Q  .(r'N')  F*(r,n) 

1  9t  n’.N'.r'  n 


X  F(r'n')  exp  -yr  Sr 


x  exp 


(46) 
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v» 


where  (nN|V|n'N')  is  given  by  equation  (34).  The  transition  rate 


from  state  r  to  r'  via  a 

jump  from  n  •*  n  +  1  is  given  by  time  dependent 

perturbation  theory: 

T*  * 

W  (nN-m'N ' 

)  =  H(r,r' ,n,n' ) 

(nN | v|n'N' ) \2 

ti 

1 

*  n 

Z  huJk(N^  -  Nk)  +  AE  J 

(47) 

where: 

AE  =  V  -  V 

_ i 

(47a) 

r 

r 

• 

lx  t\ 

1  “ 

cos  — 

(2(x) 

In  | 
xt  12 

(47b) 

ti 

9  /  \  1  ,t  ixt/ft  .  . 

—  fi(x)  =  -j  f  e  dt  (47c) 


-Y  |S  -S 

-y  , 

s  ,-s  , 

H(r,r',n,n')  = 

F*(r,n)  F(r'n')  e  rl  n  r 

r 

n  r 

(47d) 


Substituting  equation  (34)  into  (47),  we  find 


W 


r>r ' 


(nN-^n'N) 


^2  1  j2  6n  n'  +  r  H(r’r' »n»n') 

H  e=+l  n,n  £ 


* ;  \  [ 

1  "l 

l)Yk  (Nk  +  2?  sin2[k(n+e>  +  -J) 

“t  l 

- 

/ 

+  2Y, 


‘N.  -iw,  t  N.  +1  iw,  t  1 

j 

k  .  k  k  _  k 

sin  [k(n+E )  +  -^]  > 

N  N 

iAEt/fi  , 
x  e  dt 


(48) 
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ifiiaaii 


4 

Completing  the  square  in  equation  (48)  and  ignoring  the  term  in  sin  (  ) 

2 

since  it  is  of  order  1/N  ,  and  making  use  of  equations  (36a),  (36b) 
and  replacing 

Y'-fV  1  12 

W  (n+n')  -  —j  I  J  5  ,  H(r,r',n,n') 

IT  E+l  n»n  -£ 


by  (exp  gtlu^  -  1)‘ 


An  expression  valid  at  low  temperatures  will  be  of  interest 
in  discussing  the  AC  conductivity.  Inspection  of  the  second  exponential 
in  equation  (49)  shows  that  at  low  temperatures  it  will  remain  finite, 
and  hence  we  expand  it  in  a  power  series  and  obtain: 
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Now  letting  t  approach  which  is  valid  if  we  are  considering  only 

low  frequency  applied  fields;  the  terras  (sin  xt/x)  will  be  replaced 

by  27t6(x).  The  first  term  in  the  brackets  of  equation  (50)  vanishes. 

This  term  represents  transitions  which  do  not  involve  phonon  emission 

or  absorption  and  has  been  considered  previously.  The  second  term 

(single  sum)  represents  one-phonon  processes;  the  third  term  (double 

sum),  two-phonon  processes,  and  higher  order  terms  will  represent 

multi-phonon  processes.  At  low  temperatures,  only  those  modes  where 

<<  kT  will  be  activated,  and  from  equation  (33a)  we  see  that  the 

for  these  modes  is  small.  Hence  we  can  confine  our  attention  to 

2 

single  phonon  processes,  as  two  phonon  processes  are  of  order 
We  can  therefore  write: 
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m  £M 


One  Phonon  Emission: 


Wr  r  (n+n+1)  =  l-  ■  exp  [— 2S ( T)  ]  H(r ,r* ,n,n+l) 

* 


v  ^“k 

_  Yk  e  *  .  I  AE  . 

x  1  ~~ - 6 1  r  +  \ 


iN  *K 


- 1 


(51a) 


One  Phonon  Absorption: 


Wr+r  (n-mfl)  =  ^r~  exp  f-2S(T)  ]  H(r,r\n,n+1) 
■h 


2  ttJ 


Yk  i  Stimk  \-l  E 

lr\e  .  ~ 11  6  \r-\) 


(51b) 


In  considering  transport  through  a  random  lattice  we  should  have  at 
equilibrium  an  equal  number  of  emissions  as  absorptions.  With  this 
criteria  we  can  combine  equations  (51a)  and  (51b)  to  obtain: 

.2 


,r+r 


W  (n+n+1)  = 


2ttJ 


exp  j^-2S(T)j  H(r,r'  ,n,n+l) 


X 


l 


Y 

N 


k  .v 
-  coth 


Stlw, 
_ k 

2 


(52) 


To  determine  the  total  transition  rate  from  r  to  r' ,  we  must  sum  over 
all  n,  as  equation  (52)  gives  only  the  transition  rate  from  r  to  r* 
via  a  jump  from  site  n  to  site  n  +  1.  As  can  be  seen  from  equation 
(52)  the  only  n  dependent  term  is  H(r ,r' ,n,n+l) .  Let  Lhe  sites  r 
and  r*  be  separated  by  a  distance  R;  i.e.. 


S  -  S 


=  R 


(53) 
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Therefore  we  can  compute: 


Z  H(r ,r ' ,n,n+l) 
n 


Z 

n 


F*(r,n)  F(r' ,n+l)  e 


-Y  S  -S 

r  n  r 


X 


2 


(54) 


If  we  now  assume  that  most  of  the  contribution  to  the  sum  comes  from 
the  region  between  the  two  centers,  which  is  reasonable  if  R  is  great 
enough,  since  on  either  side  of  the  centers  we  will  have  the  product 
of  two  exponentially  decreasing  functions,  equation  (54)  can  be 
rewritten  as 


Z  H(r ,r ' ,n,n+l) 
n 


Z 

n 


F*(r,n)  F(r' ,n+l) |  e 


-2( y  -Y  ,)S 
r  r  n 


-2(Yr+Yr.)R  ±2(Yr-Yr»)a 

x  e  e 


(54) 


where  +a  is  the  distance  between  site  n  and  site  n  +  1.  If  we  are 
looking  at  an  average  process,  we  can  write  Yr  =  Yr»  =  Y* 

Z  H(r,r' ,n,n+l)  =  Z  |F*(r,n)  F(r',n+l)|2  e  (55) 

n  n 

There  is  no  way  of  evaluating  the  sum  on  the  right  hand  side  at  this 
time.  We  note  that  it  contains  as  its  terms  the  relative  probabilities 
that  the  jump  will  occu^  via  site  n,  rather  than  another  site  n*.  If 
R  is  large  enough,  we  should  be  able  to  approximate  it  by  a  constant 
Z,  of  order  unity.  With  this  approximation,  equation  (52)  can  be 
rewritten  as: 


W(r-Mr') 


2ttJ2  „  -2S(T) 

^Ze 


-4yR 

e 


Y,  Bhw. 

Z  —  coth  — 
k 


wk 
(56) 
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When  a  field  is  applied,  we  must  substitute  A  for  AE  in 
equation  (56): 


A  -  Vr  -  Vr,  -  eF  (57) 

Equation  (56)  can  be  evaluated  for  Case  I.  Then  from  equations  (33a) 
and  (16a)  we  can  write,  considering  coupling  from  only  the  nearest 
neighbor  oxygen  ions: 


Yk 


C2  sin2  k  sin2  2k  ,  32  C2  k  s  64  C2  ^k 

Moj2  tioj  sin"3  -r  k  Muj2  ftco  Mu)2  "fto)  Uo 

o  o  2  o  o  o  o 


(58) 


where  use  has  been  made  that  <<:  1*  Inserting  equation  (58) 

into  equation  (56),  replacing  the  sum  by  an  integral,  and  using 
relation  (57); 

W(r*r*>  -  Hi!  2  e-2S(I>  i  / 

■ft  0  Mo)  1110 

o  o 


0),  8ho), 

X  —  coth  -m— — 
u)  2 

o 


6|*±  *k 


3u,k 


CD 


(59) 


and  integrating  leads  to: 


W(r+r* )  -  ii!  z  e-2S(I>  ^  coth  H 

■ft  Mw  (Hu  ) 

o  o 


(59a) 


Equation  (59a)  is  in  the  same  form  as  that  developed  by  Miller  and 

Abrahams14  for  the  case  of  impurity  conduction  at  low  concentration, 

_2S (T) 

with  the  addition  of  the  term  e  '  ,  and  the  absence  of  a  term 
3/2 

proportional  to  R  .  Our  Z  term  could  in  fact  have  R  dependence; 
however,  it  should  vary  slowly;  the  main  variation  will  be  via  the 
exponential.  Making  use  of  typical  values5  J  ■  10  eV,  e  K  '  = 

10  2,  C2/Miij2 'fto)o  ■  5,  a>o  ■  lO3^,  equation  (59)  reads: 
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W(r+r' ) 


12 

2  x  10  AZ 


-4yR 


coth  ■— 


which  is  the  same  order  of  magnitude  as  that  given  by  PolTack  and 
Geballe16  for  Miller  and  Abraham's  expression.14 

We  note  that  the  Case  II  y^  becomes  infinite  as  •+■  0;  in 
fact  as  can  be  seen  by  inspection  of  equations  (13c)  and  (33b)  a 
uniform  translation  of  the  solid  will  give  rise  to  a  large  energy 
shift.  We  discard  this  case  as  meaningless. 


High  Temperature  Approximation 

We  will  require  an  expression  for  the  transition  rate  at  high 
temperatures  in  our  discussion  of  the  DC  conductivity.  The  integral 
3n  equation  (49)  has  been  evaluated  by  Holstein5  for  AE  ■  0.  We  shall 
follow  his  method  here.  First  we  introduce  the  transformation: 


t 


if  +  t 


(60) 


which  transforms  the  integrand  of  equation  (49)  to: 
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N 

coth  „ 
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iAET 


(61) 


The  T-integration  extends  from  -t  co  +t  on  the  real  axis,  and  two 

components  perpendicular  to  this  axis  from  -t - 2~  to  ancl  from  t 

ififi 


to  t - 5r".  Combining  the  latter  two  under  one  integral  we  obtain: 

L2 
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W(r-*-r')  -  ^  2  exp  (-4yR)  exp  (-  +  ^ 


(62) 
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j*  t 


where: 
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The  second  integral,  equation  (64),  contributes  little,  following 
Holsteins's  analysis  (that  is,  the  term  inside  the  brackets  still 
determines  the  convergence  of  the  integral) .  Equation  (63)  diverges 
when  AE  ■  0,  and  needs  to  be  corrected  by  removing  those  transitions 
for  which  N  =  N' .  Equation  (63)  therefore  becomes: 
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We  now  wish  to  let  t  -*<*>.  The  first  exponential  of  equation 
(65)  will  rapidly  converge,  while  the  second  will  cause  the  integrand 
to  oscillate.  Because  of  the  rapid  convergence,  we  can  replace 
cos  i kT  (1  -  u^x^/2): 
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Rearranging  terms,  we  find: 


=  exp  [-2S (T) ]  exp  Z  - 


AE2 


x  /  exp  I  -  Air  -  12 


liAEr  , 

-  exp  I— r —  dr  (67) 


where: 
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Evaluating  the  integral  and  letting  t  -»•  °°: 
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By  use  of  equation  (67a),  and  on  changing  the  sums  to  integrals: 


(c)  1  11  pn<i). 
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Combining  the  three  exponentials,  and  noting  that: 
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and  that  if  W  <<  kT; 
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(70c) 


(71a) 


(71b) 


In  the  crystalline  case,  Z  ■  1,  y  -  0,  AE  -  0,  and  equation  (72)  becomes 
identical  to  that  given  by  Holstein5  (his  equation  contains  an  error; 


the  tt  is  to  the  first  and  not  the  second  power).  Furthermore,  if  AE 

2 

is  small  in  comparison  to  E  ,  the  AE  term  in  equation  (72)  may  be 

cl 

neglected. 
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3.  DISCUSSION 


Review  of  Theoretical  Model 

Before  proceeding  to  apply  the  results  of  Chapter  II  to  exper- 
mental  conductivity  data,  it  is  worthwhile  to  review  the  calculations 
of  that  chapter,  with  specific  emphasis  on  the  assumptions  which  have 
been  made,  and  their  implications. 

Our  starting  equation  (10)  is  very  similar  to  that  for  the 
crystalline  lattice.5  The  essential  difference  arises  from  the  energy 
of  the  electron  in  each  of  the  complexes  being  different  due  to  the 
disorder  of  the  lattice.  We  then  approximated  the  ionic  masses  to 
be  equal,  and  that  all  nearest  neighbor  force  constants  were  equal. 

Due  to  the  anharmonic  nature  of  the  interatomic  bonds,  this  approxi¬ 
mation  is  not  strictly  true.  Inclusion  of  these  effects  will  result 
in  some  localization  of  the  vibrational  wavef unctions ,  and  hence  to 
a  nonanalytic  problem.  Mathematically,  localization  is  equivalent 
to  the  normal  vibration  coordinates  [equation  (15)]  no  longer  being 
normal.  If  the  anharmoniticity  is  not  too  great,  the  strictly 
harmonic  approximation  that  we  have  used  should  give  reasonable 
agreement  with  reality.  In  addition,  because  we  have  taken  the  ionic 
masses  to  be  equal,  the  phonon  modes  in  the  range  0  <_  k  _<  —  will  be 
accoustical  modes,  while  those  in  the  range  _<  k  <_  tt  will  be 
optical  modes. 
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An  analytic  solution  to  the  strictly  harmonic  one-dimensional 
problem  is  possible  since  the  solution  of  the  vibrational  modes  will 
be  a  plane  wave  expansion.  In  three-dimensions,  the  eigenstates  of 
the  lattice  vibrations  are  no  longer  plane  waves  due  to  the  disorder. 

It  is  for  this  reason  that  the  three-dimensional  problem  has  not  been 
solved  analytically. 10 

We  next  assumed  that  the  energy  of  the  electron  was  linearly 
related  to  the  lattice  displacements.  It  is  at  this  point  that  we 
can  distinguish  between  a  polaron  and  a  small  polaron  treatment.  In 
the  polaron  problem,  the  main  electron-lattice  interaction  would  come 
from  more  distant  neighbors  -  that  is,  the  electron  would  be  moving 
sufficiently  fast  so  that  the  lattice  in  the  vicinity  of  the  electron 
would  not  have  time  to  react  to  the  electron's  field.  In  the  small 
polaron  case,  the  electron  velocity  will  be  such  that  the  lattice 
in  the  vicinity  of  the  electron  will  react,  and  hence  the  nearest 
neighbor  interaction  will  be  important.  The  inclusion  of  distant 
interactions  enables  us  to  treat  the  polaron  case  in  our  problem. 

We  shall  come  back  to  this  point  later  on  in  this  section. 

The  assumption  that  the  overlap  integrals  J^n  are  independent 
of  the  lattice  displacements  [equation  (13d)]  removes  the  effect  of 
thermally  assisted  tunneling.  Thermally  assisted  tunneling  results 
from  the  increased  overlap  of  the  electronic  wave  functions  due  to 
lattice  vibrations  which  cause  adjacent  metal  ions  to  move  closer 
to  each  other.  Provided  that  the  wave  functions  do  not  vary  appreciably 
across  nearest  neighbor  complexes,  this  approximation  should  be  good. 

Next,  to  solve  the  wave  equation  (we  are  treating  the  coef¬ 
ficients,  a^,  as  wave  functions),  we  took  as  a  first  approximation 
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J  =  0.  In  completing  the  square  in  equation  (18)  (which  is  the 
mn 

same  thing  as  using  a  variational  principale  o  determine  the  ground 

state)  and  transforming  to  displaced  normal  coordinates,  we  allowed 

the  electron  to  interact  with  the  lattice  and  distort  it.  The  effect 

of  the  non-zero  value  for  the  J  was  treated  as  perturbation.  The 

mn  v 

Jmn's  gave  rise  to  the  matrix  elements  ( pN | V  j  p  * M ' )  which  allow  the 
polaron  to  move  from  site  to  site. 

It  is  at  this  point  that  we  have  narrowed  our  attention  to 
the  small  polaron  case.  In  our  solution,  we  have  assumed  that 

is  important,  as  it  will  be  for  the  small  polaron.  However,  by 

£ 

setting  the  C  =0  for  a  certain  number  of  inner  neighbors,  we  can 

deal  with  the  polaron  case.  Unfortunately,  we  have  no  criteria  to 

£ 

indicate  how  far  out  to  set  the  C  =0;  this  will  vary  depending  on 
the  material  considered,  and,  as  noted  in  Chapter  I,  on  the  amount 
of  disorder  present.  In  order  to  establish  such  a  criteria,  we  would 
need  to  solve  the  wave  equation  [equation  (18) ]  in  the  zeroeth  order 
for  both  J  and  the  electron- lattice  interaction.  Such  a  solution 
appears  to  be  possible  using  Anderson's  Green's  function  method;11 
however,  since  the  solution  would  require  numerical  methods,  it  shall 
not  be  considered  here. 

In  order  to  solve  the  problem,  we  assumed  that  the  interaction 

£ 

parameters  B^,  Cn  were  equal  for  all  values  of  n.  In  doing  so,  we 
have  assumed  that  the  lattice  deformation  on  different  sites,  and 
consequently  the  deformation  or  binding  erergy  is  the  same  for  all 
sites  for  a  given  set  of  phonon  occupation  numbers.  The  result  of 
nonequal  binding  energies  will  be  to  further  increase  the  energy 
differences  between  sites. 
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The  solution  to  the  crystalline  problem  reaveals  a  polaron 
bandwidth  of  2J  exp  [-2S(T)],  that  is,  one  which  decreases  exponentially 
with  temperature.  In  the  disordered  lattice,  we  apply  Anderson's 
theorem  and  find  that  the  criterion  for  localized  states  is  that  the 
parameter  (<A>/J)  exp  [S(t)J  must  be  sufficiently  large.  Here  we 
notice  that  the  criterion  is  temperature  dependent,  and  that  the 
probability  of  localized  states  increases  with  increasing  temperature. 
Furthermore,  the  presence  of  the  exponential  term,  which  is  always 
greater  than  unity  for  both  the  polaron  and  small  polaron  cases,  means 
that  localized  states  are  more  probable  in  the  polaron  problem  than 
in  the  pure  electron  problem,  and  that  the  probability  for  localized 
states,  for  a  given  amount  of  disorder,  increases  exponentially  as 
the  electron-lattice  interaction  increases.  In  addition,  there  is 
a  possibility  that  when  disorder  is  small,  there  will  be  band  states 
at  low  temperature,  but  that  these  states  will  become  localized  as 
the  temperature  increases. 

To  solve  the  problem  for  the  disordered  lattice,  we  invoked 
Mott's  conjecture1  that  the  wave  functions  decay  exponentially  as 
exp  (~  Yr|Sn-srl)*  In  other  words,  the  polaron  in  a  state  localized 
at  r  is  not  confined  to  a  single  metal  ion;  there  is  a  probability 
of  finding  it  on  nearby  ions  also.  With  our  eigenstates  thus  estab¬ 
lished,  we  determined  the  transition  rate  from  a  state  localized  on 
site  r  to  one  localized  at  r' .  Such  a  transition  cannot  be  vertical; 
that  is,  the  polaron  cannot  emit  or  absorb  a  phonon  on  one  site.  Instead 
the  transition  is  diagonal;  that  is,  the  polaron  simultaneously  jumps 
from  one  site  to  another  with  the  absorption  or  emission  of  phonons. 
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Because  of  the  overlap  of  the  eigenfunctions,  it  is  possible 
for  the  polaron  to  travel  a  large  distance  in  a  single  hop.  The 
transition  will  therefore  take  place  via  a  tunnel-hop-tunnel  mechanism. 
The  importance  of  these  transitions  follows  from  the  fact  that  nearest 
neighbor  energy  states  may  be  too  different  for  the  transition  to 
occur  easily;  but  a  transition  between  more  distance  neighbors,  where 
the  energy  difference  may  be  more  favorable,  can  occur. 

We  should  emphasize  one  fact  about  multiphonon  processes. 

These  transitions  must  take  place  in  a  single  hop,  or  via  multiple 
hops  in  which  each  intermediate  state  is  an  allowed  state.  The  proba¬ 
bility  of  a  multiple  hop  will  be  equal  to  the  product  of  the  proba¬ 
bilities  for  each  single  hop. 

We  shall  conclude  this  section  by  noting  that  the  y  appearing 
in  the  exponent  of  the  expression  for  the  transition  rate  [equations 
(59a)  and  (72)]  is  not  a  constant,  but  is  different  for  each  pair 
of  states.  This  statement  follows  from  Anderson's  Green's  function 
method, 1 1  which  will  be  treated  in  the  next  section.  This  "disorder" 
parameter,  y,  will  tend  to  zero  as  the  lattice  becomes  ordered;  hence, 
for  localized  states  in  regions  of  small  disorder,  it  will  be  small; 
while  in  regions  of  extensive  disorder,  it  will  be  large. 

AC  Conductivity 

The  AC  conductivity  of  a  disordered  solid  will  be  a  bulk  effect. 
The  applied  AC  field  will  cause  the  energy  levels  of  the  complexes 
to  fluctuate,  and  the  polaron  may  therefore  hop  from  one  site  to  another. 
Since  the  field  is  varying,  an  electron  does  not  have  to  travel  through 
the  lattice  as  is  the  case  for  DC  conduction,  but  need  merely  change 
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sites  in  order  to  contribute  to  an  AC  current.  The  result  is  that 
localized  states  which  do  not  contribute  to  the  DC  conduction  are 
brought  into  play,  causing  an  increase  in  the  observed  conductivity. 

As  the  frequency  increases ,  more  hops  occur  in  a  given  period  of  time 
(provided  that  the  transition  time  is  short  enough)  ,  and  the  net 
current  increases. 

Pollack  and  Geballe16  have  derived  an  expression  for  the 
difference  between  the  AC  and  DC  conductivity  (o^)  ,  which  is  assumed 
to  take  place  via  hopping  between  pairs  of  sites  with  an  energy 
difference  AE.  Their  expression  is  applicable  in  the  present  situation, 
by  taking  the  total  number  of  pairs  to  be  N(l-c),  where  N  is  the  total 
number  of  metal  ions,  and  c  is  the  concentration  of  electrons 
(■  Me^n^+/Me^n+'^+)  ,  and  inserting  a  term  f'(y)  dy  and  integrating 
with  respect  to  y  to  account  for  the  variation  in  the  disorder 
parameter  y: 

1  e^  2 

Re  [doAC(r ,AE ,w) ]  =  —  dp(r,AE)  N(l-c)  r 


X  / 


-1 

T 


x  ^  f 1  ( y )  dy 

.,22  .2  [AE 

1  +w  t  cosh  I  _2kx 


(73) 


The  term  dp(r,AE)  is  the  number  of  sites  of  energy  difference  AE  at 
a  distance  r,  which  we  shall  define  as: 


dp(r ,AE) 


4-rrr^  (Nc)P(r)  f(AE)  dr  d(AE) 


(74) 


where  P(r)  is  the  pair  distribution  function  of  metal  ions,  and  f(AE) 
is  the  probability  of  a  localized  state  with  energy  difference  AE.  This 
latter  function  intuitively  should  be  a  function  close  to  a  Gaussian 
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distribution  function.  Inserting  equation  (74)  into  equation  (73) 
and  integrating  over  r  and  AE,  we  find: 


2  2  2 

Re(0AC}  =  12  n2  c(1_c)  It  fff  w4t-1  7~ ~i2~2 

1  +  0)  T 


P(r)  f(AE)  f'(Y) 


.2  /  AE 
cosh  "2kT 


dr  d(AE)  dy 

where  x  is  obtained  from  equation  (59a): 
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exp  (2s(T)]  exp  4yR 


-  tanh|2kfl  <76) 

[The  term  A  is  taken  into  account  in  deriving  equation  (73).] 

Equation  (75)  is  nonanalytic  and  has  not  been  evaluated  even 
for  the  case  when  P(r),  and  the  terras  in  y  and  AE  can  be  omitted. 

We  can  discuss  the  qualitative  results,  however. 

The  first  thing  we  note  is  that  if 

1 

x  <<  — 

U 

2 

in  all  cases,  then  the  AC  conductivity  is  proportional  to  u  .  However, 
x  is  a  function  of  r,  AE,  and  y.  Pollack  and  Geballe  considered  the 
case  when  we  can  ignore  the  effect  of  y  (i.e.,  take  it  equal  to  a 
constant)  and  P(r)  [i.e.,  we  are  in  the  range  where  P(r)  i?  equal 
to  1],  They  find  that  the  AC  conductivity  is  approximately  proportional 
to  wn  where  n  is  some  constant  less  than  2.  This  behavior  arises 
from  the  circumstance  that  as  w  increases,  mere  hops  are  made  causing 
an  increase  in  conductivity,  but  the  longer  transition  times  cease 


2  2 

to  contribute  to  a  conductivity  increase  (i.e.,  — -  ^  2  '  ^  for  tliese 

1  +  0)  1 

long  transition  times) . 

The  case  when  all  factors  must  be  taken  into  consideration 

is  much  more  complex,  and  therefore  is  best  discussed  by  considering 

experimental  data.  We  consider  the  data  obtained  by  Bench  and  Clark,17 

and  Hench  and  Kinser18  for  a  potassium  vanadium  phosphate  glass. 

They  observed  that  at  low  temperatures,  the  AC  conductivity  increased 
0  8 

proportional  to  w  *  at  low  frequencies,  but  that  the  exponent  n 
gradually  increased  to  a  value  of  2  at  a  frequency  of  10^  Hz.  The 
variation  in  the  exponent  of  id0  with  a  change  in  the  frequency  can 
be  attributed  to  either  the  pair  distribution  function  P(r),  or  to 
the  y  distribution  function  f'(Y).  The  pair  distribution  function 
can  be  calculated  from  x-ray  diffraction  data,  and  is  generally  found 
to  differ  from  unity  only  for  the  first,  second  and  third  nearest 
neighbor  metal  ions.  The  distance  between  these  pairs  of  ions  is 

O 

a  maximum  of  approximately  10  A,  and  therefore  the  induced  energy 
shift  caused  by  the  field  is  of  the  order  of  10  electron  volts. 

This  value  is  much  less  than  kT,  and  hence  these  transitions  between 
second  and  third  nearest  neighbors  can  be  accomplished  by  the  phonons 
themselves,  and  should  contribute  little  to  the  AC  conductivity. 

We  are  led  to  the  conclusion  that  the  variation  in  the  exponent 
of  id  must  be  due  to  the  disorder  parameter,  v.  To  understand  how  this 
can  be,  we  must  consider  a  qualitative  solution  to  the  wave  equation 
[equation  (26)].  According  to  Anderson11 

a,,(E)  =  lim  i  o( j |W | k)  (77) 

o-*0 
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where 


4 


(j|w|k)  = 


io  +  E  -  E. 

J 


+  t: — - F"  V-  ( i  |  W  |  k)  . 

«  10  +  E  ‘  EJ  n  (78) 


The  V.  in  orr  case  is  just  the  matrix  element  (jN|V|eN).  The  energy 
J  & 

E  is  the  energy  of  the  eigenstate  V  ,  and  the  energy  is  the  energy 

of  the  polaron  on  the  complex  j  (including  the  perturbation  energy). 

The  coefficient  a.,  (E)  is  the  wave  function  at  site  k  for  a  polaron 
J  k 

localized  on  site  j.  Equation  (77)  is  solved  by  an  iteration  technique. 
From  inspection  it  is  seen  that  the  solution  for  k  distant  from  j  will 
involve  the  product  of  matrix  elements  in  the  numerator,  and  a  product 
of  energy  differences  in  the  denominator.  If  now  there  are  highly 
ordered  regions  in  the  lattice  where  the  energy  levels  are  nearly 
the  same,  then  the  wave  function  will  not  decrease  rapidly.  It  will 
decrease  rapidly  in  those  areas  where  disorder  is  high.  Thus  if  we 
may  postulate  the  structure  of  the  glass  as  being  composed  of  partially 
ordered  regions  separated  by  regions  of  large  disorder,  the  transition 
times  will  not  vary  continuously  with  r,  but  will  come  in  "clumps," 
changing  sizeably  as  the  disordered  regions  are  crossed.  A  similar 
behavior  occurs  if  instead  of  ordered  and  disordered  regions  we  postu¬ 
late  disordered  regions  of  one  phase  (metal  oxide  rich)  surrounded 
by  disordered  regions  of  a  second  phase  (metal  oxide  deficient).* 

Then  in  the  metal  oxide  deficient  region,  due  to  the  decreased  overlap, 
and  the  presence  of  other  species  presenting  a  wrong  energy  level 


*An  alternative  model  would  consist  of  two  interconnected  phases,  each 
having  a  slightly  different  metal  oxide  content,  such  that  the  easiest 
path  takes  the  polaron  across  the  different  phases,  rather  than 
restricting  it  to  one  phase. 
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to  the  electron,  the  wave  function  will  decrease  sharply;  while  in 
the  cells  of  the  disordered  metal  oxide  rich  region,  the  wave  function 
will  decrease  only  gradually.  The  exponent  variation  would  thus 
arise  essentially  from  the  distribution  of  metal  oxide  rich  or  ordered 
region  cell  sizes. 

In  order  to  be  compatable  with  x-ray  data,  the  size  of  the 

O 

ordered  regions  can  be  no  larger  than  '15  A;  otherwise  crystalline 
diffraction  lines  will  appear.  There  are  no  restrictions  on  the  size 
of  the  phase  separated  regions  however,  for  if  the  metal  deficient 
phase  assumes  a  very  small  total  volume  of  the  glass  (being  only  one 
or  two  oxide  complexes  thick  about  the  metal  oxide  deficient  phase) , 
then  low  angle  x-ray  scattering  will  not  detect  the  phase  separation. 

The  phase  separation  explanation  appears  to  be  more  desirable 
since  it  allows  the  polaron  to  travel  large  distances ,  over  which 
the  effect  of  the  applied  field  will  be  appreciable.  It  should  be 
noted  that  Pollack  and  Geballe 1 6  found  that  in  impurity  conduction 
in  silicon  at  low  concentrations  the  distance  traveled  in  a  single 

O 

hop  was  the  order  of  100-300  A.  However,  in  the  absence  of  specific 
calculations  which  our  present  model  does  not  allow,  it  is  questionable 
as  to  which  effect  is  the  cause  of  the  change  in  the  exponential  of  to. 

In  any  event,  the  electron  transport  mechanism  at  low  frequency 
would  be  similar  to  that  considered  by  Pollack  and  Geballe,  but  as 

the  frequency  increases,  we  approach  a  range  where  all  i  are  either 

11  2 
>•>  —  or  <<  — .  The  former  contribute  as  u>  as  previously  shown,  the 

0)  'Ji 

latter  do  not  contribute  to  the  frequency  dependent  conductivity. 

The  temperature  dependence  of  the  AC  conductivity  is  an 
extremely  complex  situation.  In  the  first  place,  we  note  that  the 
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transition  time  is  strongly  dependent  on  temperature  through  the  term 
2s(t)  /AE  l 

e  '  '  tanhl-^pl  The  conductivity  expression  aiso  contains  a  term 


!  /AE 
2kT 


.  Furthermore  there  are  multiple  phonon  process  which 


enter  in  as  T  increases;  these  terms  become  increasingly  important 

as  the  electron-lattice  coupling  increases.  Finally,  the  y  in  the 

exponential  will  be  temperature  dependent,  increasing  as  temperature 

increases.  This  last  point  can  be  seen  by  considering  Anderson's 

Green's  function  method.  The  amplitude  on  site  j  for  a  wave  function 

localized  on  site  k  will  be  given  by  equation  (77).  Each  of  the  matrix 

elements  (pN|v|p'N)  contains  the  term  exp  [— S (T) ] .  Therefore,  the 

sum  over  £,  which  gives  the  wave  function  at  site  j  contains  the 

term  exp  [-s(T)].  Furthermore,  because  of  (2|w|h),  the  wave  function 

_n  ^ 

must  decrease  approximately  as  e  '  for  the  n  site  from  the 

site  j.  Hence,  the  wave  function  will  become  more  localized  as  the 
temperature  increases.  This  is  the  localized  state  analogue  of  the 
decrease  in  mean  free  path  with  increasing  temperature  for  the  crystal¬ 
line  problem,  and  is  the  reason  for  the  appearance  of  the  term 
exp  [S (T) ]  in  the  criterion  for  localized  states. 

Before  proceeding,  it  will  be  of  interest  to  examine  the  term 
exp  f-s(T)].  At  absolute  zero,  this  term  will  assume  a  value  e  S^. 

As  the  temperature  increases,  this  term  will  decrease  slowly,  since 


the  coth 


will  be  large  only  for  a  few  modes.  At  some  temper¬ 


ature  (of  the  order  of  ®j)e^ye) »  a  large  number  of  modes  will  be 
activated,  and  the  exponent  will  decrease  more  rapidly.  The  larger 
the  electron-lattice  coupling,  the  more  rapid  will  be  decrease.  The 


qualitative  behavior  of  exp  [-S(T)]  is  plotted  in  Figure  4. 


The  AC  conductivity  is  given  by  equation  (75) 


Re(oAC)  =  -J2  N  c^1-c)  fli  4"r  T 


2  2 

UJ  ’ 


1  +  U) 


2  2 


P(r)  f ( AE)  f'(y) 
i  2  ME  t 

c°sh  tei 


dr  d(AE)  dy 


where : 


h2  Mu3  (-few  )2 

t  =  - ~2 - ^Y~  exp  exP  ^vR  tanh('2kfl 

256  J  Z  c 

We  shall  consider  the  effect  first  on  the  transition  time 
At  sufficiently  low  temperatures,  kT  <<•  all  AE.  Under  these  circum¬ 
stances  tanh  ("Hy)  '  exP  [2S(T)]  ;  constant,  and  y  is  not  highly 
dependent  on  temperature  (provided  that  the  electron  lattice  inter¬ 
action  is  not  too  strong).  We  conclude  then,  that  at  low  temperatures, 
the  transition  time,  ,  is  practically  independent  of  temperature. 

As  the  temperature  increases,  for  those  transitions  for 

which  kT  >>  AE,  tanh  *  0.  At  the  same  time  exp  [ 2S (T )  J  will 

increase.  Thus  some  transition  times  become  longer,  while  others 

become  shorter  (AE  and  R  small).  At  sufficiently  high  temperatures, 

AE 


tanh  exP  [  2S  (T)  J  -*■  *  for  all  transitions;  thus  all  transition 

times  increase.  The  increase  is  even  more  pronounced  due  to  the 
temperature  caused  increase  of  y.  Also,  as  the  temperature  is 
increased  multiphonon  transitions  will  come  into  play,  providing  short 
transition  times  again. 

As  previously  stated,  the  low  temperature  AC  conductivity 
supports  the  view  that  the  density  of  transition  times  come  in  clumps 
(Figure  5).  As  the  temperature  increases  therefore,  the  frequency 
at  which  the  exponent  of  v  becomes  equal  to  2  should  decrease. 
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Plot  of  density  of  transition  times. 


Furthermore,  at  frequencies  above  .10  Hz,  the  exponent  should  also 

decrease;  the  frequency  at  which  the  decrease  in  the  exponent  is  observed 

should  also  decrease  with  temperature. 

Looking  now  at  the  term  - — — — ,  at  low  temperatures, 

kT  costi 


kT  <<  all  AE ,  cosh^  |-|^rj  "*  00 •  such  that  kT  cosh^  (’ffx)"*  **  ^hat  *s» 
the  AC  conductivity  should  decrease  with  decreasing  temperature  at 
very  low  temperatures.  For  kT  AE,  kt  cosh^  "  kT,  That  is, 

at  sufficiently  high  temperatures,  the  AC  conductivity  should  decrease 
as  1/T.  At  temperatures  in  between  these  two  extremes,  only  those 
transitions  for  which  AE  «■  kT  will  contribute  appreciably  to  the 
conductivity.  As  the  temperature  increases,  more  such  states  will 
be  brought  into  play,  but  offsetting  this  increase  in  total  conductivity 
will  be  the  factor  1/T.  This  paragraph  is  depicted  graphically  in 
Figure  6. 

It  should  be  noted  that  in  this  discussion  we  have  assumed 
that  none  of  these  transitions  will  contribute  to  the  DC  conductivity. 

As  the  temperature  increases,  many  of  the  transition  with  small  AE 
will  contribute  to  the  DC  conductivity;  these  must  be  subtracted  from 
consideration  in  our  AC  conductivity. 


DC  Conductivity 

The  low  temperature  conductivity  will  be  determined  essentially 
by  the  most  difficult  step  of  the  easiest  path  through  the  glass.14 
Thus,  most  of  the  steps  through  the  glass  are  between  localized  levels 
of  nearly  the  same  energy.  It  is  probable,  however,  that  the  most 
difficult  step  will  be  of  sufficiently  high  energy  that  the  transition 
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AE  + 


Figure  6.  Transition  contribution  as  a  function  of 


will  be  via  a  multiphonon  hop.  Therefore  we  shall  resort  to  using 
the  high  temperature  approximation  which  takes  into  consideration 
these  multiphonon  hops.  A  detailed  calculation  of  the  conductivity 
must  be  based  on  averaging  procedures  similar  to  those  of  Miller  and 
Abrahams.14  In  lieu  of  the  uncertainty  of  the  exact  nature  of  the 
localized  states  at  the  present  time,  as  was  discussed  in  the  preceeding 
section,  we  shall  not  go  into  such  a  derivation  here.  Instead,  we 
shall  consider  only  the  most  difficult  step,  which  should  largely 
determine  the  conductivity.  The  conductivity  contribution  of  this 
step  will  be  from  equation  (72): 

2_2 

o  -  n  e  v  «  —  W(r-r') 


2_2t2_  1/2 

ne  R  J  Zir  ,  , 

■yj~  exp  (-4yR)  exp 


2E1/2(kT) 

a 


ae  +  M_ 

a  2  16E 


a 


kT 


(79) 


Over  a  narrow  temperature  range,  at  sufficiently  low  temper¬ 
atures  where  we  can  ignore  the  temperature  effect  on  y,  the  activation 
energy  will  be  given  by: 


-  r  .AE  ,  AE 

Q  ■  Ea  +  r  +  i6T 


(80) 


where  AE  represents  the  energy  difference  of  the  most  difficult  step. 

As  the  temperature  increases,  AE  will  decrease,  as  has  been 
noted  by  Miller  and  Abrahams,14  and  Mott.1  The  decrease  is  a  result 
of  the  activation  of  paths  other  than  the  easiest.  However,  partially 
offsetting  the  decrease  will  be  an  increase  caused  by  the  temperature 
dependence  of  y.  This  effect  will  be  more  pronounced  the  larger  the 
electron-lattice  interaction. 
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Effect  of  Heat  Treatment 


Experimentally,  the  effect  of  heat  treatment  on  a  glass  has 
been  shown  to  be  (1)  the  nucleation  and  growth  of  crystallites  in 
the  glassy  matrix,  and/or  (2)  the  separation  of  the  glass  into  two 
or  more  phases  of  differing  composition. 19 

We  shall  consider  the  effect  of  nucleation  and  growth  first, 
which  is  an  ordering  of  the  glassy  matrix.  This  ordering  must  lead 
to  a  reduction  in  the  deviation  of  the  distribution  of  the  energy 
levels,  f(AE).  The  decrease  in  disorder  must  lead  to  an  increase 
in  the  magnitude  of  the  lov  temperature  DC  conduction  due  to  the 
increased  overlap  between  more  distant  localized  states,  and  a  decreased 
transition  time.  The  activation  energy  may  also  change  to  reflect 
the  dielectric  properties  o 2  the  ordered  or  partially  ordered  lattice, 
and  the  decreased  AE  term.  As  the  disorder  continues  to  decrease, 
the  criteria  for  localized  states  will  not  longer  be  met  at  low  temper¬ 
atures  and  the  polaron  motion  will  take  place  largely  by  band  diffusion 
through  large  regions  of  the  glass.  Further  reduction  in  the  disorder 
will  result  in  the  introduction  of  more  paths  through  which  band 
diffusion  is  possible  and  an  increase  in  the  bandwidth.  At  the  same 
time  the  number  of  localized  states  would  be  reduced. 

The  effect  of  nucleation  on  the  AC  conductivity  will  be  highly 
dependent  on  the  starting  glassy  model  and  on  the  exact  way  in  which 
nucleation  takes  place.  Concrete  experimental  evidence  is  lacking 
on  either  of  these  factors.  Furthermore,  there  will  be  a  marked 
dependency  on  the  exact  nature  of  the  localized  wave  functions,  a 
factor  which  our  model  in  its  present  development  has  not  determined. 
These  three  factors  determine  the  overlap  of  the  wave  functions  as 
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a  function  of  heat  treatment.  As  the  localized  levels  disappear, 
however,  the  AC  effect  will  disappear. 

Some  comment  is  in  order  on  the  two  models  postulated  in  the 
AC  conductivity  section;  that  is,  the  order-disorder  cell  model,  and 
the  two-phase  cell  model.  In  the  order-disorder  model,  the  partial 
ordering  of  the  regions  which  are  highly  disordered  will  lead  to  a 
decrease  in  the  reduction  of  the  wave  function  amplitude  as  it  traverses 
the  disordered  region  (Figure  7).  On  the  average,  the  effect  on  the 
wave  function  is  equivalent  to  a  decrease  in  the  disorder  parameter  y. 
Alternatively,  in  the  two-phase  cellular  model  (or  in  a  model  where 
the  phase  interconnect  as  discussed  in  the  AC  conductivity  section)  , 
order  in  the  metal  oxide  deficient  phase  will  have  the  same  effect 
as  in  the  order-disorder  model. 

A  phase  separation  from  a  homogeneous  phase  into  a  metal  oxide 
rich  and  a  metal  oxide  deficient  region  will  also  ip.ad  to  an  increase 
in  the  DC  conductivity  due  to  an  increased  overlap  of  the  electronic 
wave  functions.  A  continuous  phase  will  be  formed  provided  that* 

(1)  the  separation  occurs  by  spinodal  decomposition  or  (2)  separation 
occurs  by  nucleation  and  growth,  and  that  the  metal  oxide  rich  phase 
assumes  approximately  40  percent  or  more  of  the  total  volume  of  glass. 

In  the  latter  case,  the  continuous  phase  (as  well  as  the  second  phase) 
is  nucleated  as  spherical  particles  in  the  glassy  matrix.  These  spher¬ 
ical  particles  grow,  eventually  intersecting  to  form  a  continuous 
phase.  In  spinodal  decomposition,  compositional  fluctuations  arise.. 

As  the  heat  treatment  progresses,  the  compositional  gradient  across 
the  phase  boundary  sharpens,  until  complete  separation  has  occurred. 


*A  review  of  phase  separation  in  a  glass  is  given  by  Hench.19 
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It  is  probable  that  the  new  phases  will  be  more  ordered  than 
the  parent  phase.  The  reasons  for  this  are:  (1)  The  parent  phase 
reflects  the  temperature  from  which  the  glass  is  cast  and  the  quenching 
rate  (as  well  as  composition) ,  while  the  disorder  of  the  separated 
phases  will  strongly  reflect  the  annealing  temperature  (as  well  as  the 
composition  of  the  phases),  which  is  lower  than  the  casting  temperature. 
(2)  The  "driving  force"  for  the  phase  separation  is  a  reduction  in 
the  free  energy  caused  by  satisfying  the  preferred  coordination  numbers 
of  the  species  participating. 

The  result  of  such  a  phase  separation  will  be  an  increase 
in  DC  conductivity.  The  region  of  rapidly  increasing  conductivity 
will  be  the  region  where  nulceation  and  growth  is  taking  place,  or 
where  the  spinodal  structure  is  coarsening.  The  activation  energy 
will  change  (either  increase  or  decrease)  to  reflect  the  dielectric 
properties  of  the  high  metal  oxide  phases,  and  the  new  disorder 
parameter  and  energy  differences. 

The  effect  of  phase  separation  on  the  AC  conductivity  again 
will  be  dependent  on  the  model  treated.  If  the  compositional  gradient 
across  the  phase  boundary  is  sharp  enough,  then  conduction  in  the 
high  metal  oxide  phase  may  be  via  band  diffusion,  which  the  metal 
oxide  deficient  phase  will  act  as  an  insulator.  Thus  after  sufficient 
heat  treatment,  all  conduction  would  be  essentially  confined  to  the 
metal  oxide  rich  phase,  via  band  diffusion,  and  the  AC  conductivity 
would  disappear. 

The  above  discussion  was  based  on  the  assumption  that  the 
initial  composition  was  homogeneous,  and  is  therefore  applicable  only 
to  the  order-disorder  cell  model.  Turning  our  attention  now  to  the 


two-phase  cellular  model  and  the  alternative  interconnected  phase 
model,  we  must  take  into  consideration  the  fact  that  phase  separation 
has  already  begun.  Thus  the  two-phase  cellular  model  would  have  to 
stem  from  a  rather  complicated  peretectoid  or  eutectic  reaction. 

We  shall  therefore  not  consider  it  further.  If,  however,  decomposition 
has  already  begun  spinodally,  then  the  coarsening  of  the  phase  compo¬ 
sition  can  lead  to  the  formation  of  band  states  in  the  metal  oxide 
rich  phase  and  a  loss  in  the  AC  conductivity. 

Turning  now  to  the  experimental  data  (Figure  8) ,  we  see  that 
the  DC  conductivity  rapidly  increases  with  heat  treatment.  The  leveling 
off  of  the  DC  conductivity  after  a  short  time  indicates  that  (1)  a 
path  of  continuous  high  metal  oxide  content  has  been  formed  or  (2)  that 
a  path(s)  has  been  established  which  permits  conduction  essentially 
by  band  diffusion. 

The  AC  conductivity  during  the  early  stages  of  heat  treatment 
supports  either  of  these  views.  The  increase  in  the  magnitude  of 
the  AC  conductivity  indicates  that  more  localized  states  can  contribute 
to  the  conduction.  This  can  be  a  result  of  either  an  expansion  in 
the  spatial  extent  of  the  localized  wave  functions,  allowing  hops 
to  be  made  between  more  distant  localized  states  for  which  the  field 
will  cause  a  greater  variation  in  energy,  and/or  a  decrease  in  the 
energy  levels.  The  loss  of  the  AC  effect  with  sufficiently  long  heat 
treatment  indicates  that  the  localized  levels  have  spread  out  to  such 
an  extent  that  they  may  be  considered  as  band  states.  Any  of  the 
models  discussed  so  far  supports  these  observations. 

Additional  insight  into  the  distribution  of  transition  times 
is  gained  in  considering  the  effect  of  temperature  on  the  AC  conductivity 
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of  a  glass  heac  created  for  a  short  time.  The  AC  effect  Is  found 
to  be  present  at  low  temperatures,  but  decreases  in  magnitude  as  the 
temperature  increases  (Figure  9).  At  a  temperature  of  56°C,  the  AC 
conductivity  becomes  less  than  the  DC  conductivity.  Since  AC  conduction 
Is  via  the  bulk  of  the  glass,  the  only  way  in  which  the  AC  conductivity 
can  be  less  than  the  DC  conductvitiy  is  if  the  DC  conductivity  at 
temperatures  above  the  inversion  temperature  is  also  a  bulk  effect, 
and  not  confined  to  a  limited  number  of  paths. 

The  depression  in  conductivity  will  arise  from  the  presence 

in  the  DC  conduction  expression  of  transition  times  r  <  — .  That  is, 

0) 

both  in  AC  and  DC  conduction,  hops  are  made  mainly  by  thermal  acti¬ 
vation,  but  as  the  frequency  of  the  applied  field  increases,  the  long 
transition  times  will  not  contribute  to  the  conductivity.  This  is 
a  different  phenomenon  from  the  low  temperature  AC  enhancement  effect; 
in  that  case,  the  localized  states  which  caused  an  increase  in  the 
conductivity  with  frequency  did  not  contribute  appreciably  to  the 
DC  conductivity.  Instead,  they  were  activated  by  a  variation  in  the 
energy  levels  caused  by  the  applied  AC  field.  Both  effects  may  be 


present;  we  require  only  that  the  suppression  effect  outweigh  the 
enhancement  effect  at  high  temperatures. 

Referring  to  equation  (76)  ,  we  see  that  long  transition  times 
will  arise  from  (1)  the  effect  of  the  exponential  term  exp  [2S(T)] 


which  increases  with  temperature;  (2)  the  exponential  term  exp  (+4yR) , 

AE 

since  y  increases  with  temperature;  and  (3)  the  term  (tanh  -jj^)  which 


increases  with  increasing  AE. 


When  a  DC  field  is  applied,  the  metal  oxide  complexes  assume 
a  fixed  set  of  energies.  Assuming  bulk  conduction,  hops  between  all 


;ure  9.  AC  conduct 


states  will  contribute  to  the  total  conductivity.  We  will  also  assume 
that  band  diffusion  contributes  little  at  the  high  temperatures.  When 
an  AC  field  is  applied,  the  transitions  for  which  t  <  will  contribute 

little  (i.e.,  proportional  to - —  ^  .  The  AC  enhancement  effect 

1  +  (j  r 

may  still  be  present;  we  require  only  that  the  suppression  effect 
outweighs  the  enhancement  effect.  The  long  transition  times  arise 
primarily  from  distant  states  with  large  energy  difference. 

As  the  heat  treatment  continues,  the  localized  states  become 
less  and  less  localized  until  band  diffusion  appears.  This  will  mean 
that  the  temperature  at  which  the  inversion  in  the  AC  conductivity 
takes  place  will  increase,  primarily  since  v  decreases.  Similarly 
at  early  stages  during  the  heat  treatment,  the  enhancement  effect 
outweighs  the  depression  effect  at  all  but  the  highest  temperatures, 
when  all  states  participate  in  the  hopping.  Once  again  we  see  that 
either  the  order-disorder  cellular  or  the  two-phase  models  are  qual¬ 
itatively  correct  in  dealing  with  the  depression  phenonemon. 

Conclusions 

1.  The  localized  polaron  in  the  Nearly  Crystalline  Ionic  Lattice 
model  is  found  to  qualitatively  describe  the  observed  conduction 
properties  of  a  potassium  vanadium  phosphate  glass. 

2.  The  conduction  phenomenon  of  a  glass  is  found  to  be  highly 
structure  dependent,  and  with  further  numerical  development 
should  yield  useful  quantitative  information  on  the  disorder 
and  phase  separation  in  the  glass. 
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Mott's  proposal1  that  the  localized  wave  functions  vary  as 
exp  (-y|xn-xj)  is  found  to  bo  true  only  on  the  average.  In 
actuality,  the  wave  functions  appear  to  be  highly  sensitive 
to  the  exact  disorder. 

The  AC  conductivity  of  a  potassium  vanadium  phosphate  glass 
can  be  explained  in  terms  of  an  inhomogeneous  structure;  the 
inhomogeneities  being  a  result  of  either  a  phase  separation 
or  varying  degress  of  order  and  disorder. 

Future  experimental  work,  should  be  centered  around  proving 
or  disproving  the  existence  of  a  phase  separation.  This  work 
should  include  small  angle  x-ray  scattering,  transmission 
electron  microscopy,  and  electrical  measurements  of  splat- 
cooled  V2O5  (where  a  phase  separation  is  impossible).  Future 
theoretical  work  hould  center  on  numerical  solutions  to  the 
wave  functions,  based  on  the  proposed  models.  The  theory  should 
also  be  extended  to  glasses  containing  other  transition 
metal  oxides. 
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B.  THE  EARLY  STAGES  OF  CRYSTALLIZATION  OF  GLASS  (L.  L.  Hench,  D.  L.  Kinser 
and  S.  W.  Freiman) 

1.  INTRODUCTION 

It  has  been  well  known  for  many  years  that  when  glass  is  heated  at  a 
temperature  in  the  softening  range,  physical  and  chemical  properties  change 
irreversibly  and  eventually  lead  to  complete  crystallization.  Tamman"*" 
established  that  the  bulk  crystallization  sequence  in  glasses  consists  of 
two  distinct  processes,  nucleation  and  crystal  growth.  These  processes  may 
occur  separately  or  concurrently  during  the  crystallization  of  glasses.  Thermo¬ 
dynamic  theories  have  since  been  developed  that  describe  the  kinetics  of  nucle- 

2  3 

ation  and  growth  in  terms  of  activated  complexes  *  .  The  effect  of  hetero¬ 
geneities  on  nucleation  and  growth  kinetics  have  also  been  described  by 

thermodynamic  arguments.  Application  of  these  concepts  has  been  the  subject 

4  5  6 

of  considerable  research  ’  ’  and  has  led  to  the  development  of  a  new  class 

7,8 

of  glasu-ceramic  materials 

However,  in  spite  of  the  above  advances,  there  is  still  considerable 
question  as  to  the  structural  mechanisms  involved  in  the  nucleation  of  crystal 
phases  in  a  glass.  In  particular,  little  is  known  about  the  relationship  be¬ 
tween  the  mechanisms  of  nucleation  and  growth  and  the  rates  of  bulk  crystal¬ 
lization  of  glasses.  There  are  at  least  several  reasons  responsible  for  this 
state  of  affairs,  including  the  added  complexity  of  including  phase  separation 
in  structural  theories,  the  difficulty  of  interpreting  changes  in  multicomponent 
commercial  glass  systems,  and  perhaps  most  important,  the  scale  of  observation 
required  to  observe  nucleation  processes. 

It  is  the  objective  of  the  present  paper  to  establish  a  model  for  the 
nucleation  process  iu  binary  L^O-SiC^  an{*  Na20-SiC>2  glasses  based  on  a 
variety  of  analytical  results.  A  second  objective  is  to  relate  the  nucleation 
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process  to  bulk  crystallization  kinetics  of  the  glass  and  the  microstructure 
developed  during  devitrification. 

2.  EXPERIMENTAL  PROCEDURE 

Sample  Preparation:  The  glass  was  prepared  for  these  experiments  by 
blending  99. 9%  silica  (5P  Min-u-Sil)  and  99.3%  lithium  carbonate  in  a  jar 
mill  and  melting  for  18-70  hours  in  air  in  a  covered  platinum  crucible.  The 
glass  compositions  and  melting  temperatures  employed  are  given  in  Table  I. 
Samples  2.5  cm  long  and  1.3  cm  in  diameter  of  1.8  cm  diameter  and  0.5  cm 
thick  were  cast  in  a  steel  mold,  placed  in  a  furnace  at  a  temperature  range 
of  230°  C  -  300°  C  and  heated  for  1-4  hours  to  relieve  stresses  to  a  level 
sufficient  to  prevent  sample  breakage.  The  cast  specimens  were  subsequently 
given  nucleation  and  crystallization  heat  treatments  in  nichrome  wound  tube 
furnaces  equipped  with  aluminum  blocks  to  eliminate  temperature  gradients. 
Temperature  control  was  maintained  at  ±  1#C. 

Dielectric  Relaxation  and  D.  C.  Conductivity:  Dielectric  absorption 
spectra  were  measured  on  polished  glass  samples  with  vapor  deposited  gold 
electrodes  over  a  frequency  range  of  20  Hz  to  5  MHz  using  Wayne  Kerr  bridges. 
The  source  and  detector  system,  sample  holder  configuration  and  environmental 

9 

control  system  have  been  discussed  in  a  previous  publication  .  Dielectric 
loss  peaks  were  measured  as  a  function  of  temperature  over  a  range  of  -80°  to 
200®  C  with  a  control  of  ±  1°  C.  The  oscilloscope,  VTVM,  and  picoammeter 

system  employed  to  make  D.  C.  conductivity  measurements  within  a  range  of 

5-4  9 

times  of  10  sec  to  10  sec  have  also  been  described  previously  . 

X-Ray  Analysis:  Bulk  crystallization  kinetic  data  were  obtained  on  -200 
mesh  powders  according  to  the  technique  of  Munch  and  Pierron  (10)  using  a 
Norelco  x-ray  diffractometer  with  a  fine  focus  Cu  tube  and  curved  crystal 
monochromator.  Quantitative  analysis  of  the  x-ray  data  has  been  discussed  in 
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a  recent  publication 

X-ray  identification  of  the  crystals  produced  in  the  early  stages  of 

the  reaction  sequence  was  obtained  through  the  use  of  a  vacuum  Guinier- 

12 

DeWolff  camera  using  CuKa  radiation 

X-ray  small  angle  scattering  (SAS)  measurements  were  made  on  samples 

thinned  to  0.1  to  0.2  mm  thickness  after  thermal  treatments  in  the  form  of 

13 

1.3  cm  diameter  by  2.5  cm  long  cylinders  .  A  Siemans  Kratky  camera  using 

CuKa  radiation  and  entrance  slits  of  150u  and  exit  slits  of  50p  was  used  to 

measure  the  scattering  intensity  over  the  range  of  .01°  -  .2°  20. 

Microscopy:  Microscopy  of  the  reaction  sequence  occurring  in  the  early 

stages  of  crystallization  was  obtained  using  hot  stage  transmission  electron 
14 

microscopy  .  Samples  were  reduced  to  transmission  thickness  using  a  special 
14 

etching  procedure  .  A  Philips  EM  200  was  used  for  bright  field,  dark  field, 
and  diffraction  photographs. 

Analysis  of  the  microstructures  of  the  intermediate  stages  of  crystalli¬ 
zation  was  made  using  carbon  replicas  of  Pt  preshadowed  etched  fracture  sur¬ 
faces  . 

A  Bausch  and  Lomb  Dynoptic  polarizing  microscope  and  a  Leitz  research 
metallograph  were  used  for  quantitative  analysis  of  the  microstructures  of 
the  final  stages  of  crystallization. 

3.  EVIDENCE  FOR  A  METASTABLE  REACTION  SEQUENCE 
The  compositions  of  the  glasses  studied  in  the  Li20-Si02  system  are 
shown  on  the  phase  diagram  of  the  system  in  Figure  1.  Only  the  disilicate 
composition  in  the  Na20-Si02  system  has  been  investigated  to  date.  The  as- 
cast  glasses  prepared  from  the  compositions  outside  the  disilicate  contain 

microheterogeneous  regions,  as  seen  in  the  replica  electron  micrograph  of 

12 

Figure  2.  It  has  been  concluded  in  this  laboratory  and  by  other  investi¬ 
gators'*^’^  that  the  heterogeneous  microstructure  consists  of  Si02  droplets 
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Figure  2  Replica  electron  micrograph  of  30  mole  %  lithia- 
silica  glass  as  cast- 
-70- 


in  a  disilicate  matrix.  It  has  also  been  concluded  *  that  a  secondary 
phase  separation  of  alkali-rich  regions  occurs  within  the  disilicate  matrix. 

Of  concern  is  what  happens  structurally  when  the  microheterogeneous 
glass  is  heat  treated  at  temperatures  in  the  nucleation  range.  Quantitative 
dielectric  absorption  has  been  an  important  tool  in  following  such  structural 
changes  because  the  alkali-rich  regions  of  the  glass  gradually  lower  their 
free  energy  by  increasing  in  modifier  content,  forming  isolated  islands  of 
near  alkali  metasilicate  composition.  When  an  A.C.  field  is  applied,  the 
presence  of  the  high  conductivity,  modifier-rich  regions  in  a  low  conductivity 
Si02-rich  matrix  gives  rise  to  interfacial  polarization.  The  interfacial 

polarization  results  in  large  Maxwell-Wagner-Sillars  (MWS)  heterogeneous 
12 

dielectric  losses 

It  can  be  observed  in  Figure  3  that  the  magnitude  of  the  loss  peaks  devel¬ 
oped  in  the  initially  loss-free  30  mole  percent  Li20-Si02  glass  decreases 
as  the  duration  of  thermal  treatment  is  increased.  The  decrease  in  absorption 
indicates  that  either  the  size  of  the  highly  conductive  regions  is  decreasing 
or  that  the  compositional  difference  between  the  modifier-rich  regions  and 
the  Si02-rich  phases  is  decreasing  with  thermal  treatment. 

A  similar  sequence  of  appearance  and  disappearance  of  loss  peaks  has  also 

18 

been  observed  in  other  alkali-silicate  glasses  .  The  summary  of  these  results 
in  Figure  A  also  shows  that  the  changes  in  loss  peak  height  vary  for  the  dif¬ 
ferent  glasses.  Consequently,  it  appears  that  the  kinetics  of  precipitation 
and  resorption  of  the  modifier-rich  phase  varies  with  the  modifier  content  of 
the  parent  glass. 

Small  angle  x-ray  scattering  is  also  sensitive  to  heterogeneous  structures 
19 

in  glasses  and  can  be  used  to  determine  the  size  of  the  heterogeneous  regions, 

20 

the  composition  of  the  regions,  and  a  number  of  other  structural  features 
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A  parameter  (M)  directly  proportional  to  the  total  integrated  small 
angle  x-ray  scattering  intensity  is  plotted  as  a  function  of  thermal  treat¬ 
ment  of  the  lithia  silicate  glasses  in  Figure  5.  M  is  defined  as  M  =  Q  /K' 

M  =  Q0/K'  and  K'=  4itt q/A  where  tq  is  a  factor  involving  the  distance  between  the 
sample  and  the  detector  slit,  and  A  is  the  x-ray  wavelength.  Figure  5  shows 
that  the  scattering  intensity  increases  to  a  maximum  with  thermal  treatment 
followed  by  a  decline  in  scattering.  The  magnitude  of  scattering  intensity 
is  a  function  of  the  total  volume  of  the  scattering  particles  and  the  electron 

22  7 

density  difference  between  the  particles  and  the  matrix  ;  Qq  =  V  (p  -  pQ)  . 
Therefore,  it  appears  that  either  the  total  volume  of  the  scattering  particles 
is  changing  with  thermal  history  or  the  compositional  difference  between  the 
particles  and  the  matrix  is  changing.  The  total  volume  of  particles  can  in¬ 
crease  by  two  alternative  methods.  The  particles  can  grow  during  thermal 
treatment,  causing  the  average  radii  to  increase.  Or,  new  particles  may  be 
grown  during  thermal  treatment  which  would  cause  the  size  distribution  of 
particles  to  change. 

Table  I  represents  a  summary  of  the  Guinier  radius  and  the  Porod  radius 

measured  from  the  scattering  curves  for  various  heat  treatments.  Two  features 

of  the  data  can  be  noted.  First,  since  both  sets  of  measured  radii  are  close 

to  the  same  value,  the  size  distribution  of  the  scattering  particles  is  very 
23 

narrow  and  does  not  change  with  time.  Secondly,  the  size  of  the  scattering 
particles  does  not  change  appreciably  with  thermal  treatment. 

Consequently,  the  stability  of  the  particle  size  and  the  size  distribution 
indicates  that  the  change  in  scattering  intensity  is  due  to  the  compositional 
difference  between  the  scattering  particles  and  the  matrix  becoming  greater 
during  the  first  part  of  the  thermal  treatment,  followed  by  a  reduction  in 
compositional  differences  in  the  later  stages.  Thus,  the  x-ray  small  angle 
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Figure  5a  Transmission  electron  micrograph  of  30  mole  % 

lithia-silica  glass  at  480°C  after  35  minutes  at 
480°C  in  the  microscope.  Bright  field 


Figure  5b  Transmission  electron  micrograph  of  30  mole  X 

lithia-silica  glass  at  480°C  after  35  minutes  at 
480°C  in  the  microscope.  Dark  field. 


scattering  results  are  consistent  with  the  interpretation  that  metastable 

metasilicate  particles  precipitate  in  the  glass  and  then  resorb  as  the 

equilibrium  crystal  phase  appears. 

Confirmation  that  it  is  the  alkali-metasilicate  that  appears  as  the 

isolated  metastable  phase  in  the  above  experiments  was  obtained  from  both 

x-ray  and  transmission  electron  microscopic  analysis.  Lithium  metasilicate 

crystals  were  detected  in  the  30  mole  percent  Li20-Si02  glass  using  a  Guinier- 

12 

DeWolff  vacuum  x-ray  diffraction  camera  .  A  concentration  of  approximately 
0.1%  was  observed  at  heat  treatment  times  corresponding  to  the  maximum  scat¬ 
tering  intensity  and  maximum  dielectric  losses.  Additional  diffraction  lines 
in  the  metasilicate  pattern  indicate  that  a  transition  phase  of  lower  symmetry 
or  extensive  solid  solution  may  be  present.  With  longer  thermal  treatments 
the  metasilicate  phase  disappears,  within  the  detection  limits  of  the  diffrac¬ 
tion  system,  leaving  only  diffuse  glass  peaks.  For  longer  thermal  treatments 
of  up  to  50  hours  at  500°C,  the  equilibrium  lithium  disilicate  phase  appears. 

Figures  6  and  7  show  the  appearance  of  the  metasilicate  crystals  in  a  30 
mole  percent  Li20-Si02  glass  after  a  thermal  treatment  of  the  electron  trans¬ 
mission  specimen  for  five  minutes  at  480°C.  A  bright  field  image  and  the  dif¬ 
fraction  pattern  of  the  metasilicate  crystals  are  shown.  Indexing  the  diffrac¬ 
tion  pattern  shows  that  some  disilicate  crystals  are  *lso  present  (Table  III). 

A  complete  sequence  of  the  precipitation  and  resorption  process  will  be  reported 

in  another  publication  including  a  discussion  of  the  increased  reaction  rate 

14 

due  to  beam  heating  and  sample  thickness. 

The  simultaneous  appearance  of  both  lithium  metasilicate  and  lithium  di¬ 
silicate  in  selected  area  diffraction  patterns  also  provides  evidence  that  the 
resorption  of  the  metasilicate  is  not  completed  before  the  equilibrium  phase 
begins  to  grow.  This  observation  lends  support  to  the  latter  argument  that 
the  disilicate  growth  proceeds  at  the  decomposing  metasilicate  interface. 
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Figure  7  The  effect  of  thermal  treatment  time  at  500°C  on 
a  parameter  related  to  the  total  integrated  small 
angle  intensity  for  two  compositions  in  the 
Li  O-SiO  system 


As  summarized  in  Table  IV,  all  alkali-silicate  compositions  investigated 
appear  to  follow  a  metastable  reaction  sequence  similar  to  that  observed  for 
the  30  mole  percent  Li20-Si02  glass.  The  major  distinction  between  the  various 
compositions  is  that  glasses  of  the  disilicate  composition  do  not  undergo  a 
liquid-liquid  phase  separation  prior  to  the  metastable  precipitation. 

By  using  available  thermodynamic  data,  it  has  betn  possible  to  construct 

12  18 

a  free  energy-composition  diagram  for  the  various  glasses  studied  *  .  One 

of  the  diagrams  reproduced  in  Figure  8  shows  that  the  metastable  reaction 
sequence  observed  experimentally  is  consistent  with  a  sequential  lowering  of 
the  free  energy  of  the  system.  Crystallization  of  the  alkali-rich  regions 
such  as  the  metasilicate  is  favored  because  of  its  structural  simplicity 
relative  to  the  glass  structure.  The  average  free  energy  of  the  system  is 
decreased  by  the  amount  G3  -  Gg  through  the  metasilicate  crystallization. 

Reaction  of  the  metasilicate  with  the  Si02~rich  matrix  produces  the 
structurally  more  complex  disilicate  crystals  and  further  lowers  the  free 
energy  by  the  amount  Gg  -  Gg.  Thus,  it  is  the  kinetic  barriers  to  nucleation 
of  the  structurally  complex  equilibrium  phases  that  is  responsible  for  the 
intermediate  reaction  sequence  observed. 

4.  RELATION  OF  THE  METASTABLE  REACTION  TO  CRYSTALLIZATION  KINETICS 

A  recent  publication  by  two  of  the  authors  demonstrated  that  the  crystal¬ 
lization  kinetics  of  a  33  mole  percent  Li20-Si02  glass  which  has  had  a  three 
hour  nucleation  treatment  at  475°C,  insufficient  to  form  a  significant  concentra¬ 
tion  of  the  metasilicate  phase,  is  appreciably  slower  than  a  glass  of  the  same 
composition  which  has  had  a  24  hour  nucleation  treatment  at  475°C^.  The  24 
hour  treatment  produces  a  large  concentration  of  metastable  metasilicate 
crystals.  The  accelerating  effect  of  the  metasilicate  phase  on  the  crystal¬ 
lization  rate  of  the  glass  is  shown  in  Figure  9.  The  rate  of  crystallization 
is  increased  by  a  factor  of  1.5.  However,  it  can  also  be  noted  that  the  incu- 
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Figure  8  Free  'anergy-composition  diagram  for  lithia-silica 
system  at  500°C. 
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Figure  9  Crystallization  of  L^Si^O^  from  a  Li20-2Si02 
glass  for  two  nucleation  treatments 


bation  period  is  constant  for  the  two  samples  at  all  temperatures.  There¬ 
fore,  the  activation  energy  of  the  incubation  process  for  both  samples  is 
independent  of  the  metastable  reaction  process. 

The  activation  energy  for  bulk  crystallization  is  strongly  dependent  on 

prior  nucleation  treatment.  The  activation  energy  for  samples  with  a  24 

hour  treatment  corresponds  to  that  of  crystal  growth  in  these  systems 

(52  kcal/mole) ,  while  a  sample  with  a  three  hour  nucleation  treatment  exhibits 

a  much  higher  activation  energy  (92  kcal/mole).  Application  of  the  Johnson- 
24  25 

Mehl-Avrami  [JMA]  *  analysis  to  the  crystallization  data  reveals  that  the 
growth  of  the  disilicate  crystals  is  rod-like^,  as  is  observed  in  the  spheru- 
lites  grown  in  the  glass  studied.  The  JMA  analysis  also  showed  that  it  is 
the  activation  energy  of  nucleation  that  is  decreased  when  the  sample  has  been 
given  a  heat  treatment  that  produces  metasilicate  as  a  nucleation  step. 

The  question  to  be  answered  is  what  is  the  role  of  the  metasilicate  particles 
in  controlling  the  kinetics  of  nucleation  of  the  equilibrium  disilicate  crystals. 
Figure  10  depicts  a  structural  model  that  is  proposed  to  describe  the  relation 
between  the  metastable  reaction  sequence  and  the  equilibrium  crystallization.  The 
structural  model  shows  the  structural  simplicity  of  the  metasilicate  lattice  and 
the  structural  similarity  of  the  disilicate  lattice  when  alternate  L^O  groups 
are  removed.  Thus,  it  is  proposed  that  metasilicate  crystals  form  in  the  glass 
by  diffusion  of  L^O  into  L^O-rich  regions  because  of  preferred  Li  coordinations 
and  the  relative  structural  similarity  of  the  metasilicate  to  the  initial 
glass  structure.  With  additional  heating,  L^O  groups  diffuse  out  of  the  meta¬ 
silicate,  probably  as  ionic  species,  producing  a  disilicate  structure  from 
the  SiC^-rich  matrix  at  the  metasilicate  interface.  It  is  this  diffusional 
step  which  is  responsible  for  the  incubation  period  measured.  When  a  critical 
size  of  disilicate  region  is  reached,  the  growth  proceeds  rapidly  towards 
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Figure  10  Structural  model  of  the  lithium  metasilicate 

nucleation  of  lithium  disilicate  crystals.  Upper 
chain  structure  is  a  planar  representation  of  the 
Li20-Si02  lattice.  The  lower  sheet-like  structure 

is  a  planar  projection  of  the  Li20*2Si02  lattice 


completion.  Because  of  the  difficulty  of  L^O  groups  diffusing  through  the 
disilicate  structure,  it  is  energetically  easier  for  the  diffusion  to  occur 
at  the  ends  of  the  dendrite  arms,  thus  producing  the  spherulitic  structures 
occurring  in  these  glasses. 

When  a  glass  has  been  given  a  nucleation  treatment  which  develops  a  large 
concentration  of  the  highly  ordered  metasilicate  regions,  rapid  crystal  growth 
proceeds  from  many  sites.  The  final  spherulite  size  is  small,  the  number  per 
unit  volume  is  large,  and  the  size  distribution  is  small,  as  observed.  As 
the  spherulites  envelop  the  ordered  regions  developed  around  each  metasilicate 
site,  the  nucleation  barrier  in  the  growth  of  new  dendrite  arms  is  lowered, 
as  observed,  and  the  overall  bulk  crystallization  rate  is  accelerated,  also 
as  observed.  This  is  the  process  the  authors  refer  to  as  internally  activated 
nucleation. 

The  basis  for  internally  activated  nucleation  may  be  seen  thermodynamically 

26 

if  one  looks  at  the  equation  given  for  heterogeneous  nucleation  rate  . 


j=max  Nj 

=  T  tt  K. 

3-0  v  J 


AF*  +  AF 
exp  -  J-  a 


kT 


I  is  the  nucleation  rate;  Nj/V  is  the  number  of  nucleation  sites  per  unit 
volume;  is  the  frequency  factor;  AF*  is  the  energy  necessary  to  form  a 
stable  nucleus  at  a  particular  site  j  ;  AF  is  the  activation  energy  for  a  mole- 
cule  to  detach  itself  from  the  glass  and  attach  itself  to  the  growing  nucleus. 
At  large  degrees  of  supercooling,  such  as  occurs  in  the  crystallization  of 
glass,  AF  is  the  predominant  term  and,  therefore,  governs  the  rate  of  crystal- 

SL 


lization.  It  is  essentially  this  energy  that  is  measured  as  the  activation 
energy  for  nucleation.  It  is  easy  to  see  that  ordering  around  the  metasilicate 


regions  should  reduce  the  energy  barrier  which  must  be  overcome  in  order  to 
transfer  a  molecule  from  the  glass  to  the  crystalline  state,  thereby  reducing 
the  activation  energy  for  nucleation. 

When  the  nucleation  treatment  has  been  insufficient  to  develop  a  large 
concentration  of  ordered  metasilicate  regions,  the  incubation  time  for  crystal 
growth  remains  the  same  as  the  internally  activated  process  because  the 
diffusional  step  must  still  occur.  The  incubation  activation  energy  remains 
the  same,  100  kcal/mole  as  observed.  Much  fewer  sites  reach  critical  size, 
however,  and  the  bulk  crystallization  rate  is  therefore  slow.  The  spherulites 
do  not  envelop  ordered  regions  during  growth  and  the  nucleation  activation 
energy  is  high. 

5.  SUMMARY 


The  appearance  of  metastable  phases  as  a  step  in  a  transformation  se¬ 
quence  is  not  without  precedent.  Ordered  regions  in  super-saturated  metallic 

solid  solutions  called  Guinier-Preston  zones  appear  and  eventually  disappear 

27  28 

prior  to  equilibrium  precipitation  ’  .  Roy  proposed  metastable  reaction 


processes  in  the  thermal  treatment  of  glass  several  years  ago 


29 


McMillan 


and  Partridge  suggested  that  the  presence  of  a  metastable  lithium  metasilicate 

30 

phase  influenced  the  properties  of  Li20-Zn0-Si02  glass-ceramics  .  Recently, 

several  investigators  have  reported  variations  in  properties  with  thermal 

treatments  which  cannot  be  explained  by  liquid-liquid  immiscibility  effects 

31,32 

and  imply  that  metastable  phases  are  present 

Consequently,  it  appears  that  the  metastable  reaction  sequence  investi¬ 
gated  in  the  Li20-Si02  and  Na20-Si02  glasses  may  be  generally  observed  as  a 
precursor  to  equilibrium  crystallization.  It  also  appears  likely  that  the 
influence  of  ternary  or  multicomponent  additions  in  the  crystallization  behavior 
of  alkali-silicate  glasses  could  be  one  of  either  enhancing  or  retarding  the 
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development  of  the  metastable  phase.  For  example,  the  retardation  of 
crystallization  ol  L^O-SiC^  glasses  by  3%  TiC^  additions  is  apparently  due 
to  the  metasilicate  failing  to  appear  because  of  the  titanium  ion  compe¬ 
tition  for  oxygen  ions'*''*'.  When  the  metastable  phase  does  precipitate,  it 
is  a  very  efficient  nucleating  site;  consequently,  it  is  possible  to  control 

particle  sizes  through  the  use  of  internally  activated  nucleation  without 

33 

adding  heterogeneous  agents  to  the  glass 

Another  consequence  of  the  metastable  reaction  sequence  is  the  possibility 
of  drastic  structural  and  property  changes  in  a  glass  without  crystals  appear¬ 
ing  in  a  quantity  detectable  by  routine  diffraction  methods.  This  may  mean 
that  glasses  which  have  been  given  an  extended  annealing  at  high  temperatures 
or  have  been  heated  for  long  durations  at  elevated  temperatures,  such  as 
diffusion  experiments  or  many  property  measurements,  can  yield  data  character¬ 
istic  of  a  heterogeneous  glass-crystal  morphology  and  not  that  intrinsic  to 
the  glass  itself.  Detailed  documentation  of  the  thermal  history  of  all  glass 
experiments  must  therefore  be  deemed  essential  in  light  of  these  results. 
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III.  Measurement  Techniques  and  Analysis  (R.  W.  Gould,  J.  J.  Hren,  L.  L. 

Hench,  and  J.  K.  Watson 

A.  X-RAY  ANALYSIS  (R.  W.  Gould  and  J.  J.  Hren) 

The  x-ray  research  conducted  during  this  period  comprises  two  related 
efforts:  1)  supporting  analytical  analysis  and  2)  basic  x-ray  oriented 
research.  Both  will  be  discussed  in  this  report. 

The  basic  x-ray  research  has  been  largely  devoted  to  particle  size  dis¬ 
tribution  analysis  as  determined  by  small  angle  x-ray  scattering.  This  work, 
which  has  been  under  development  for  approximately  one  year,  has  been  stimu¬ 
lated  by  the  need  for  such  information  in  connection  with  the  degenerate 
semiconductor  research  program. ^  Small  angle  x-ray  scattering  research  has 
been  augmented  by  the  optical  diffraction  analog  approach,  which  will  also 
be  discussed  in  this  report. 

A  new  laboratory  for  x-ray  topographic  analysis  (Lang  technique)  has 
been  set  up  during  the  summer  of  1968.  This  laboratory  will  be  utilized  for 
basic  and  support  research  in  the  degenerate  semiconductor  program. 

A  significant  portion  of  the  research  effort  of  the  x-ray  group  has  been 
devoted  to  the  solution  of  analytical  problems  arising  from  various  THEMIS 
research  groups.  The  most  notable  problem  has  been  the  quantitative  micro¬ 
probe  analysis  of  Permalloy  thin  films  (See  Section  III-B) .  In  addition,  crystal 
orientation  analysis  has  been  performed  for  various  groups. 

1.  PARTICLE  SIZE  ANALYSIS  BY  SMALL  ANGLE  X-RAY  SCATTERING 
Small  angle  x-ray  scattering  has  been  used  to  find  the  average  dimensions 
of  particles  in  dispersed  suspensions.  It  is  especially  convenient  in  studies 
of  small  particles  with  diameters  of  the  order  of  30  to  5,000  angstroms.  A 
knowledge  of  the  average  particle  dimension  is  often  insufficient  to  character¬ 
ize  a  particle  system,  and  information  must  be  obtained  on  the  distribution 
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of  particles  if  meaningful  correlation  with  properties  is  to  be  obtained. 

X-ray  methods  for  the  determination  of  particle  distribution  analysis 

2  3 

have  been  attempted  for  some  20  years.  The  early  work  of  Schull  ,  Roess  » 

4 

and  Risemann  independently  showed  that,  in  principle,  small  angle  x-ray 
scattering  could  be  used  to  determine  the  distribution  of  particle  diameters 
in  dilute  suspensions  of  spherical  particles.  It  must  be  emphasized  that 
the  systems  in  which  particle  distribution  analysis  has  been  successful  have 
been  dilute  systems  containing,  in  general,  less  than  2  vol.  %  of  dispersed 
phase.  This  early  work  of  Roess  and  Risemann  was  never  applied  to  experi¬ 
mental  scattering  data  since  at  the  time  of  their  work  (1940's  and  early 
1950's)  small  angle  scattering  theory  and  techniques  were  not  sufficiently 

developed  to  make  the  method  practical. 

5,6 

Recent  investigations  by  Schmidt  and  his  students  suggest  that,  by 
modifying  the  earlier  work  of  Roess  and  Risemann,  diameter  distributions  can 
be  found  from  small  angle  x-ray  scattering  provided  that  the  systems  are 
dilute  and  the  particles  are  equiaxed . 

The  restriction  for  dilute  particle  densities  is  based  upon  the  need  for 
independent  particle  scattering.  When  particle  densities  exceed  a  certain 
lower  limit,  interferences  occur  between  scattering  originating  at  the  indi¬ 
vidual  particles.  The  interparticle  interference  causes  great  difficulty  with 
interpretation  of  the  small  angle  scattering  data,  especially  with  regard  to 
the  determination  of  particle  size  distribution  by  the  integral  transform 
method . 

Schmidt  and  his  students  have  approached  the  matter  of  particle  size 
distribution  from  an  elegant  point  of  view  that  was  originally  taken  by  Rcess. 
This  school  obtains  the  particle  size  distribution  directly  from  the  integral 
transform  of  the  small  angle  scattered  intensity.  This  method  has  the  advan- 


-87- 


tage  that  it  requires  no  assumptions  concerning  the  type  of  particle  distri¬ 
bution  found  in  the  real  systems.  The  integral  transform  is  conveniently 
evaluated  by  numerical  analysis  techniques,  and  Schmidt  has  tested  the  prac¬ 
ticality  of  his  method  on  theoretical  scattering  curves  and  some  experimental 
scattering  curves  of  Ludox  (small  spherical  particles  of  silica  of  known  size 
distribution) • 

Small  angle  x-ray  scattering  has  been  applied  to  pre-precipitation  metal¬ 
lurgical  systems  for  the  past  ten  years.  These  systems  contain  5  to  10%  by 
weight  of  scattering  particles  and  thus  cannot  be  considered  dilute.  In 
general,  interparticle  interference  effects  are  clearly  evidenced  in  the 
small  angle  scattering  curves  of  metallic  alloys.  Thus  the  methods  described 
by  Roess,  Risemann  and  Schmidt  are  clearly  out  of  their  theoretical  bounds  when 
alloy  systems  are  being  considered. 

For  this  reason,  the  work  in  this  laboratory  has  been  largely  devoted  to 
determining  a  system  which  could  provide  the  particle  diameter  distributions 
in  systems  where  interparticle  interference  did  occur.  Recntly,  Harkness, 

Gould  and  Hren  have  addressed  this  problem  and  derived  a  method  for  obtaining 
particle  diameter  distributions  from  systems  such  as  aluminum-zinc  and  aluminum- 
silver,  where  the  particles  are  spherical  but  the  system  is  not  dilute.  This 
method  has  been  described  in  a  recent  paper  and  will  be  briefly  outlined  below. 

The  basis  of  this  system  is  a  presupposition  concerning  the  form  of  the 

particle  distribution  function.  This  particle  distribution  function  is  then 

related  to  parameters  which  can  be  obtained  from  the  small  angle  scattering 

7 

curve.  Harkness,  Gould,  and  Hren  showed  that  the  size  distribution  function, 
in  this  case  log  normal,  could  be  written  in  terms  of  simple  small  angle  scat¬ 
tering  parameters.  Equation  (1)  is  the  log  normal  size  distribution  function. 
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The  geometric  mean  and  the  variance  o  may  be  expressed  as  shown  in  Equation  (2) 
in  terms  of  a  general  moment  equation  of  the  log  normal  distribution  function. 
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Baur  and  Gerold  showed  (for  a  line  shaped  x-ray  beam)  that  the  Guinier 

k 

radius,  R  »  could  be  expressed  by  Equation  (3)  and  the  Porod  radius,  R  , 
G 

by  Equation  (4) . 
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These  equations  are  simply  ratios  of  moments  of  the  particle  size  distribution 

function.  Thus,  if  the  size  distribution  function  is  assumed,  these  ratios 

* 

may  be  easily  calculated.  When  thj.s  is  done  for  the  log  normal  distribution, 


Equations  (5)  and  (6)  are  obtained. 


* 

R  =  exp (In 

U  + 

2.5  In2  o) 

(5) 

Rq  =  exp (In  \ 

j  + 

6  ln2o) 

(6) 

when  Equations 

(5)  and  (6)  are  solved  for  V 

and 

0,  one  obtains 

the  following 

expressions : 

Rr 

In  =  In  R  - 

G 

1.71 

ln  “I 

R 

(7) 

*  While  this  is  true  for  the  log  normal  distributions,  calculations  of  this 
nature  for  the  normal  distributions  prove  to  be  extremely  difficult. 
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This  then  gives  the  form  of  the  size  distribution  curve  in  terms  of  the  small 

angle  scattering  parameters  alone,  for  if  M  and  o  are  substituted  into 

Equation  (1),  f(x)  may  be  obtained  for  any  x  desired. 

9 

During  the  summer  of  1968  Gould  and  Bates  have  applied  this  method  of 

analysis  to  the  determination  of  larger  particle  sizes,  such  as  those  found 

in  dispersion  hardened  nickel  alloys.  The  details  of  this  work  appear  in 

Section  III-A-2.  This  effort  demonstrated  the  applicability  of  this  method 

to  larger  particle  sizes  than  are  normally  found  in  the  GP  zone  metal  systems. 

10 

In  addition,  Harkness  and  Gould  have  also  utilized  this  method  to  follow 
the  growth  of  Guinier-Preston  zones  by  following  the  change  in  the  particle 
size  distributions  with  time.  This  is  a  major  advance  over  previous  methods 
in  which  only  one  particle  diameter  could  be  monitored  during  the  growth  of 
a  dispersed  phase. 

5 

Utilizing  the  data  obtained  by  Letcher  and  Schmidt  ,  Gould  in  some  unpub- 

6 

lished  results  has  applied  the  HGH  analysis  to  known  particle  distribution 

systems,  such  as  Ludox*  spheres.  Figure  1  shows  the  small  angle  scattering 

intensity  obtained  by  Schmidt's  students  for  three  Ludox  samples.  Figures 

4 

2  and  3  show  the  particle  distribution  results  from  the  Letcher-Schmidt 
analysis  (integral  transform  method)  and  Figure  4  is  the  corresponding  result 
obtained  by  the  simpler  HGH  analysis. 

In  Ludox  I  (Figure  2)  the  most  probable  particle  diameter  as  obtained  from 
the  manufacturer  was  approximately  70  angstroms.  For  Ludox  III  (Figure  3)  the 
average  diameter  obtained  by  Schmidt  was  55  angstroms,  while  the  manufacturer's 
average  diameter  was  140  angstroms.  Note  that  the  present  analysis  (Figure  4) 

*  Ludox  samples  were  colloidal  silica  obtained  from  Dupont  Chemical  Corporation 
whose  particle  size  distribution  has  been  determined  by  electron  microscopy. 
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Figure  1 

Slit  Corrected  x-ray  Scattering  Curves  for  Ludox  Samples. 
Curves  A,  B,  and  C  represent  Ludox  III,  Ludox  II  and  Ludox  I 
respectively  (after  Letcher  and  Schmidt  .(^) 


Figure  2 

Particle  Size  Distribution  for  Ludox  I 
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Figure  3 

Particle  Size  Distribution  for  Ludox  III 


Particle  Diameter  (A  ) 


shows  an  average  particle  diameter  of  approximately  40  X,  which  is  also  lower 
than  the  manufacturer's  data. 

Note  also  that  the  curves  given  by  Schmidt  do  not  appear  to  have  a  real¬ 
istic  size  distribution,  particularly  at  the  large  particle  end  of  the  size 

distribution.  This,  Schmidt  explains,  is  due  to  a  series  termination  effect 

6 

and  has  been  corrected  somewhat  in  a  later  paper  .  In  order  to  carry  out 
the  analysis  of  the  Schmidt  data,  it  was  necessary  to  rederive  the  equations 
given  in  the  HGH  paper  since  the  Schmidt  data  had  been  "desmeared"  to  obtain 
collimat-on  corrected  data.  These  equations  for  the  mean  and  the  variance 
of  the  log  normal  distribution  for  "desmeared"  small  angle  scattering  inten¬ 
sity  are  given  below. 
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The  Harkness-Gould-Hren  method  has  also  been  applied  to  some  additional 
Ludox  scattering  data  obtained  from  a  paper  by  Rothwell^ .  Figure  5  shows 

the  small  angle  scattering  curves  for  various  concentrations  of  type  AS  Ludox, 

20 

where  h  =  21T(— ).  Note  the  interparticle  interference  effects  appearing  in 

the  0.16  vol.  %  sample,  as  evidenced  by  the  small  interference  peak  at 
h  =  2.6  x  10~2. 

Figure  6  shows  the  effect  that  concentration  has  upon  the  Guinier  radius 
which  is  often  used  to  characterize  a  particle  system.  Note  that  in  this 
figure  it  is  given  as  a  diameter  where  D  =  2x(Guinier  radius).  Note  also 
that  the  diameter  changes  markedly  with  concentration  and,  in  this  case,  the 
diameter  decreases  as  the  concentration  of  particles  increases. 

Figure  7  shows  the  small  angle  scattering  curves  for  glassy  carbon  con¬ 
taining  a  rather  high  concentration  of  small  spherical  voids.  The  curve 
marked  with  the  small  circles  represents  a  sample  that  was  heat  treated  at 

3000°  C,  producing  voids  having  a  Guinier  radius  of  approximately  25  X  and 

-18 

containing  approximately  5  x  10  particles. 

One  of  these  curves  was  analyzed  with  the  HGH  method  described  in  the 
previous  paragraphs  and  the  size  distribution  obtained  is  shown  in  Figure  8. 
The  analysis  of  all  of  these  results  has  not  been  completed,  but  indications 
are  that  the  method  described  by  Harkness,  Gould,  and  Hren  can  be  applied  to 
particle  size  distribution  analysis  in  non-dilute  systems.  It  appears  to 
suffer  somewhat  from  interparticle  interference  effects,  probably  due  to  the 
effect  of  interparticle  interference  on  the  determination  of  the  Guinier 
radius . 

Figures  5A  and  5B  show  the  result  of  the  HGH  analysis  when  applied  to  the 
Ludox  scattering  data  given  by  Rothwell  (10). 
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Figure  5 

Small  Angle  X-Ray  Scattering  Curves  for  Various 
Concentration  Solutions  of  Type  AS  Ludox,  h  =  2tt  (20 /X 

( Pof  m\  3  '  ' 


Figure  6 

Effect  of  Concentration  of  Ludox  on  Particle 
Diameter  Measurement.  Diameter  =  2x  (Guinier  radius^ 
Ref.  10)  ’ 
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Figure  7 

Small  Angle  Scattering  data  for  Glassy  Carbon  after 
Rothwell  (10) 


Particle  Diameter  (A) 


a.  OPTICAL  DIFFRACTION  MODELING  IN  PARTICLE  SIZE  DISTRIBUTION  ANALYSIS 

When  applied  to  x-ray  iiffraction,  Fourier  transform  theory  states  that 
the  observed  diffraction  pattern  is  directly  related  to  the  Fourier  transform 
of  the  real  crystal.  This  Fourier  transform  of  the  real  crystal  may  be  shown 
to  consist  of  a  product  of  Fourier  transforms  of  various  factors:  a  peak 
function  representing  the  crystal  lattice,  a  shape  function  to  represent  the 
size  and  shape  of  the  crystal,  and  a  function  representing  the  contents  of 
the  single  repeat  unit  in  real  space.  When  the  Fourier  transforms  of  these 
functions  are  convoluted  with  the  Fourier  transform  of  the  several  instrument 
factors,  the  result  is  directly  proportional  to  the  observed  diffraction  pat¬ 
tern  . 

It  may  be  shown  that  this  approach  is  not  limited  to  crystalline  materials, 
but  may  easily  be  applied  to  other  states  of  matter.  Bragg  in  1939  constructed 
a  diffraction  grating  consisting  of  holes  drilled  in  an  opaque  plate  at  the 
relative  positions  of  reciprocal  lattice  points.  When  light  was  transmitted 
by  these  holes,  the  diffraction  pattern  was  formed.  This  diffraction  pattern 
is  remarkably  similar  to  the  x-ray  diffraction  pattern  obtained  form  the  real 
crystal  whose  reciprocal  lattice  net  was  used  as  a  model  for  the  optical  mask. 

This  correspondence  is  an  exposition  of  the  fundamental  of  image  formation  ex¬ 
pounded  by  Abbe  in  1873. 

This  analogy  between  the  scattering  of  light  by  a  series  of  regularly  spaced 
holes  and  the  diffraction  of  x-rays  or  electrons  by  regularly  spaced  atoms  on 
a  periodic  lattice  has  been  utilized  to  understand  diffraction  theory  more 
fully.  The  Fourier  transform  method  requires  that  the  position  of  the  various 
scattering  centers  (atoms,  molecules,  clusters,  etc.)  be  specified.  Equation  (16) 
shown  below  then  gives  the  scattering  wave  amplitude,  G(§),  from  n  scattering 
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units  whose  positions  arc*  given  by  a  vector  (r^)  . 

N 

G(S)  =  fn(S)  2tt  i?n-  t  (16) 

n=l 

The  scattering  amplitude  may  be  evaluated  in  a  computer  and  the  resulting 
function  G(S)  may  be  obtained  over  a  portion  of  reciprocal  space. 

A  far  easier  and  more  convenient  way  of  carrying  out  this  Fourier  trans¬ 
form  is  by  means  of  optical  diffraction  analogs.  Utilizing  this  method, 
models  may  be  made  by  simple  drafting  techniques,  photographically  reduced, 
and  a  diffraction  pattern  obtained  in  a  matter  of  minutes.  Thus  rapid  testing 
of  various  models  is  facilitated.  An  optical  diffraction  laboratory  has  been 
set  up  in  the  Department  of  Metallurgical  and  Materials  Engineering  by  Dr. 

R.  W.  Gould.  At  the  heart  of  the  laboratory  is  a  Rank-Pullin  optical  diffract¬ 
ometer  which  has  been  modified  to  be  used  with  a  gas  laser  light  source.  This 
device  is  being  used  to  model  particle  distribution  scattering  theory.  Scatter 
ing  from  such  lattice  heterogeneities,  including  dislocations,  vacancy  clusters 
and  GP  zones  is  being  compared  with  theoretical  scattering  equations  and  with 
experimental  results  from  small  angle  x-ray  scattering. 

12 

Such  optical  diffraction  modeling  is  not  new.  Hosemann  ,  Taylor  and 
13 

Lipson  ,  among  others,  have  been  active  in  this  field  although  it  has  not 
obtained  the  attention  that  it  justly  deserves.  Recently,  the  optical  diffrac¬ 
tion  method  has  been  applied  to  the  interpretation  of  low  energy  electron 

14 

diffraction  patterns  ,  providing  a  new  and  powerful  method  for  the  unraveling 

of  these  sometimes  complex  surface  diffraction  effects.  Figure  9  shows  some 

optical  diffraction  models  of  systems  containing  spherical  particles  of  varying 

12 

concentration  and  poiydispersity .  These  figures  were  obtained  by  Hosemann 
and  the  corresponding  optical  diffraction  patterns  are  given  below  each  model. 
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Figure  9 

Optical  Diffraction  Models  and  Pattern:  for  Spherical 
Particle  Systems  of  Varying  Concentration  (after  Hosemann  11) 


IP 
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Microdensitometer  curves  through  these  optical  diffraction  patterns  (after 
suitable  correction  for  intensity)  are  shown  in  Figure  10.  These  curves  are 
analogous  to  small  angle  x-ray  scattering  curves.  Thus,  in  this  case,  the 
particle  system  is  completely  under  control  and  one  can  analyze  scattering 
effects  from  a  predetermined  distribution  of  sizes  with  a  far  greater  freedom 
than  is  available  in  natural  particle  systems. 

b.  APPLICATION  TO  DEGENERATE  SEMICONDUC  ’  ORS 
The  optical  diffraction  method  will  be  used  to  model  the  distribution  of 
scattering  particles  in  degenerate  silicon  and  germanium.  Observed  scattering 
curves  will  be  interpreted  with  the  aid  of  optical  models,  such  as  the  ones 
shown  in  Figure  9.  Note  that  a  large  number  of  scattering  centers  are  required 
to  simulate  natural  conditions.  The  models  shown  in  Figure  9  contain  an 
average  of  2,000  or  more  spherical  particles.  In  order  to  eliminate  diffrac¬ 
tion  "ghosts"  it  is  necessary  to  rotate  these  masks  in  the  optical  diffracto¬ 
meter  during  exposure. 
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Figure  10 

Microdensitometer  Trace  of  Middle  Optical  Diffrac¬ 
tion  Pattern  of  Figure  9-  Log  intensity  plotted  versus 
diffraction  angle. 
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2.  DETERMINATION  OF  PARTICLE  SIZE  DISTRIBUTION  IN  A  DISPERSION  HARDENED 
NICKEL  ALLOY  USING  SMALL  ANGLE  X-RAY  SCATTERING  (R.W.  Gould  and  S.R. 
Bates) 

a.  INTRODUCTION 

An  effective  method  for  strengthening  a  ductile  matrix  is  the  inclusion 
of  a  finely  dispersed  second  phase.  Physical  properties  of  such  dispersion 
strengthened  materials  are  strongly  influenced  by  the  particulate  nature  of 
the  dispersant.  Particle  parameters  such  as  size,  shape,  number  fraction 
and  distribution  of  particle  sizes  will  influence  the  physical  properties  of 
the  final  products  made  from  such  alloys.  Many  of  these  particle  paremeters 
are  difficult  to  measure  owing  to  the  small  size  of  most  dispersants.  A 
typical  example  may  be  found  in  dispersion  strengthened  nickel  base  alloys 
such  as  the  "TD  nickels."  In  this  example,  finely  divided  thoria  (Th09)  is 
incorporated  into  pure  nickel  or  a  nickel  alloy.  The  thoria  particles  are 
generally  less  than  1000  angstroms  in  diameter,  making  it  necessary  to  use 
transmission  electron  microscopy  (TEM)  to  observe  any  direct  particle-matrix 
interactions . 

Particle  size  distribution  measurements  by  TEM  are  time  consuming  and 
tedious;  thus,  in  most  cases  only  a  few  thousand  particles  will  be  measured. 
This  paper  will  describe  another  technique  for  obtaining  particle  size  distri¬ 
butions  from  dispersion  strengthened  alloys  such  as  TD  nickels.  The  technique 
to  be  described  in  this  paper  utilizes  small  angle  scattering  to  determine 
the  particle  parameters. 

b.  THEORETICAL  BACKGROUND 

General  small  angle  x-ray  scattering,  hereafter  SAXS,  theory  has  been 

1,2 

discussed  thoroughly  in  the  literature  .  In  this  section  we  will  only 

present  those  SAXS  equations  that  have  particular  bearing  on  the  problem  at 

3 

hand.  Letcher  and  Schmidt  have  used  an  integral  transform  of  the  scattered 
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intensity  to  obtain  particle  size  distribution  of  polydisperse  spherical 
particles.  Their  method  does  not  require  assumptions  about  the  form  of  the 

distribution  function  but  is  more  difficult  to  apply  than  the  method  described 

4 

in  this  paper.  Harkness  et  al  have  recently  shown  that  in  certain  particu¬ 
late  systems  it  is  possible  to  obtain  the  particle  size  distribution  in  a 
very  simple  manner  using  SAXS  if  one  apriori  assumes  a  size  distribution 
function.  Naturally  the  particle  system  must  meet  all  the  requirements  for 
SAXS.  Furthermot e ,  the  analysis  has  been  carried  out  only  for  spherical 
particles  having  either  a  log  normal  or  a  normal  distribution  of  particle 
diameters.  In  this  section  we  will  briefly  discuss  the  theoretical  back¬ 
ground  of  this  type  of  analysis. 

Size  Distribution  and  Small  Angle  X-ray  Scattering 
4 

Harkness  has  shown  that  a  size  distribution  function,  in  this  case 
log  normal,  could  be  written  in  terms  of  simple  SAXS  parameters.  Equation  (1) 
is  the  log  normal  size  distribution  function. 


F(x)  = 


(/2tt)  (x)  (In  o) 


exp 


-  —  In  P  -  In  x 

2 


In  o 


(1) 


The  geometric  mean  and  the  variance  may  be  expressed  as  shown  in 
Equation  (2)  in  terms  of  a  general  moment  equation  for  the  log  normal  distri¬ 
bution  function. 


2 

rn  =  exp  n  In  P  +  In  o 


(2) 


Baur  and  Gerold^  showed  that  the  Guinier  radius,  R^,  could  be  expressed  by 
Equation  (3)  and  the  Porod  radius,  R  »  by  Equation  (4). 


<R^> 

<r5> 


(3) 
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(4) 


These  equations  are  simply  ratios  of  moments  of  the  particle  size  distri¬ 
bution  function.  Thus,  if  the  size  distribution  function  is  assumed,  these 
ratios  may  be  easily  calculated.  When  this  is  done  for  the  log  normal 
distribution,  Equations  (5)  and  (6)  are  obtained. 


R  =  exp (In  p  +  2.5  ln^  o) 


Rg  =  exp (In  p  +  6  lnz  a  ) 


When  Equations  (5)  and  (6)  are  solved  for  and  one  obtains  the  following 


expressions . 


In  p  =  In  R„  -  1.71  In  — 
G  R* 


In  RG 


In2  o 


This  then  gives  the  form  of  the  size  distribution  curve  in  terms  of  the  SAXS 
parameters  alone,  for  if  P  and  d  are  substituted  into  Equation  (1),  F(x) 
may  be  obtained  for  any  x  desired.  This  analysis  has  also  been  carried  out 
for  a  normal  distribution. 

Small  Angle  X-ray  Scattering  Parameters 

It  is  evident  then  that  only  two  SAXS  parameters  are  necessary  to  charac¬ 
terize  the  size  distribution  and  these  parameters  may  be  obtained  easily  from 
line  collimated  SAXS  data  provided  that  the  small  angle  camera  is  properly 
aligned.  R^  is  obtained  using  the  familiar  Guinier  approximation  for  spherical 
particles,  Equation  (9),  where  e  is  the  scattering  angle  in  radians. 


I(£)  =  n  exp 


4*2  rg  "2 
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Hr 


7  * 

Gerold  has  shown  that  R  .  the  so-called  Porod  radius,  can  be  obtained 
from  line  collimated  SAXS  data  if  the  "infinite"  beam  conditions  are  met. 
This  relationship  is  given  in  Equation  (10)  where  s  =  e/>  is  the  modulus  of 
s,  the  diffraction  vector,  and  sQ  is  a  value  of  s  lying  in  the  Porod  region 
of  the  SAXS  curve,  i.e.,  where  the  slope  of  In  I(s)  versus  In  of  s  has  a 
value  three  for  oo 


* 

R 


fo 


8  71  (1  -  f) 


s  I(s)  ds 


I(s0) 


(10) 


In  Equation  (10)  f  is  the  volume  fraction  of  scattering  particles.  The 

integral  in  Equation  (10)  is  evaluated  using  a  finite  sum  with  an  appropriate 

7  8 

termination  approximation  ’  . 

c .  EXPERIMENTAL  RESULTS 

Alloy  Preparation 

A  sample  of  TD  NiC  (Ni-2.13  wt.  %  Cr-2.6  wt .  %  Th02>  was  kindly  provided 
by  Dr.  B.  A.  Wilcox  of  the  Metals  Science  Group  at  Battelle  Memorial  Insti- 

9 

tute,  Columbus,  Ohio.  Dr.  Wilcox  has  studied  the  creep  behavior  of  this 
alloy,  and  in  the  process,  has  measured  the  particle  size  distribution  using 
TEM  (see  Figure  1).  Other  investigators^  have  measured  the  particle  size 
distributions  in  basic  TD  nickel  alloys  containing  no  chromium.  Some  typical 
results  on  alloys  of  this  latter  type  are  shown  in  Figure  2.  Unfortunately, 
alloys  such  as  these  were  not  available  for  this  research. 

The  alloy  used  in  this  research  was  cold  rolled  to  a  thickness  of  0.10 
millimeters  and  stress  relieved  at  1000°  C  in  vacuum  for  two  hours.  Prior 
to  x-ray  examination  the  sample  was  chemically  thinned  to  a  uniform  thickness 
of  0.027  ±  0.003  millimeters. 
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Figure  1.  Transmission  electron  micrograph  of  Ni  »  20  Cr  -  2  Th02 

used  in  this  research,  after  Wilcox  et  al9. 
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Figure  2.  Particle  size  distribution  of  ThC>2  particles  in  a 
dispersion  hardened  nickel  alloy  containing  no 
chromium,  after  Mee  and  Sinclair1^. 
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Experimental  Results 

Table  I  presents  the  small  angle  scattering  results  of  this  investigation. 

* 

The  Guinier  radius,  R_,  and  the  Porod  radius,  R  ,  were  obtained  using  Equations 

G 

(9)  and  (10),  respectively.  Ny  was  calculated  from  Equation  (11)  below,  using 
Equation  (2)  to  obtain  a  value  for  <R^>.  The  volume  fraction,  f,  of  TtU^ 
particles  was  2.2%  in  the  alloy  used  in  this  research. 


(ID 


The  parameters  given  in  Table  I  have  been  used  to  obtain  the  quantities 
In  M  and  ln^  o  by  means  of  Equations  (7)  and  (8).  These  values  were  then 
substituted  into  Equation  (1)  from  which  the  particle  size  distribution  curve 
(Figure  3)  was  established. 

d.  DISCUSSION 

Wilcox^  obtained  a  particle  size  distribution  for  the  same  alloy  that 
was  used  in  this  research.  The  particle  diameters  that  he  measured  by  trans- 
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TABLE  I 


Figure  3 


Small  Angle  X-ray  Scattering  Results  from  the 
TD  Ni-C  Alloy  Used  in  This  Rasesrch 


Rg  -  110  +  10  A 

R*  -  60  +  5  A 

16  3 

Nv  ■  2.  A  x  10  particles  per  cm 


Size  distribution  of  Th02  particles  in  TD  NiC  measured 
by  small  sngle  x-ray  scattering. 


mission  microscopy  ranged  from  25  angstroms  to  460  angstroms,  with  the 
average  diameter  being  145  angstroms.  These  data  are  not  inconsistent  with 
the  results  obtained  in  this  study.  As  can  be  seen  from  Figure  3,  the  SAXS 
results  give  an  average  particle  diameter  of  approximately  80  angstroms,  with 
particle  diameters  ranging  from  20  angstroms  to  250  angstroms.  These  diffe¬ 
rences  may  have  been  the  result  of  the  resolution  of  the  SAXS  arrangement  being 
too  low  to  measure  particles  near  the  large  end  of  the  size  distribution.  In¬ 
vestigations  are  currently  underway  using  a  high  resolution  entrance  slit 
arrangement  on  the  Kratky  camera.  It  is  alsp  evident  that  the  position  of 
the  maximum  in  the  size  distribution  curve  determined  by  SAXS  is  very  sensi¬ 
tive  to  Rq  and  to  the  ratio  of  Rq/R*.  as  can  be  seen  in  Equation  (7).  In  the 
present  research,  we  did  not  obtain  a  long  straight  line  portion  in  the 
In  1(e)  versus  e  plot  used  to  obtain  the  Guinier  radius  [see  Equation  (9)]. 
This  occurs  as  a  result  of  the  particle  size  distribution  and  interparticle 
interference  effects  and  will  cause  some  uncertainty  in  the  determination  of 
R„  and,  therefore,  in  the  mean  (u)  of  the  log  normal  distribution. 

U 

e.  SUMMARY 

The  results  of  this  research  indicate  that  the  small  angle  x-ray  scat¬ 
tering  method  can  be  used  to  estimate  the  particle  size  distribution  of  thoria 
particles  in  dispersion  strengthened  nickel  alloys  if  a  particle  size  distri¬ 
bution  function  can  be  assumed.  The  method  described  is  rapid  and  may  be 
easily  applied  to  other  alloys  of  this  type.  Current  investigations  are 
underway  on  nickels  containing  no  chromium. 
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B.  ELECTRON  MICROPROBE  ANALYSIS  OF  PERMALLOY  THIN  FILMS  (R.  W.  Gould  and 
J.  K.  Watson) 

The  fundamental  problem  in  the  quantitative  analysis  of  Permalloy  thin 
films  resides  in  the  thickness  of  the  samples:  they  are  too  thin  to  utilize 
the  normal  quantitative  microprobe  technique  applicable  to  an  "infinitely" 
thick  specimen. 

There  are  few  techniques  described  in  the  literature  which  are  capable 
of  analyzing  alloy  thin  films.  Thick  films  (greater  than  3,000  angstroms) 
will  often  be  analyzed  using  the  infinitely  thick  approach  if  the  alloy 
linear  coefficient  is  large  enough.  Ultra-thin  films,  less  than  500  angstroms, 
may  be  analyzed  by  the  so-called  infinitely  thin  sample  technique  (which 
neglects  absorption)  if^  the  take-off  angle  of  the  emergent  x-rays  is  suffi¬ 
ciently  large.  Unfortunately,  the  microprobe  in  the  Department  of  Metallurgy 
has  an  18°  take-off  angle  which  necessitates  films  less  than  200  to  300 
angstroms  if  the  infinitely  thin  sample  method  is  to  be  used  successfully. 

A  general  "infinitely  thin"  technique  for  binary  alloys  of  A  and  B 
atoms  has  been  developed  and  is  described  by  Equation  (1)  given  below. 

XA  *  Va(1  +  V(A6S)Pz  (1) 


where, 


W, 


the  proportionality  constant  which  contains  factors  such  as 
crystal  and  detector  efficiency,  etc.  Unless  element  A  and 
element  B  are  near  neighbors  on  the  periodic  chart,  will 
vary  markedly  with  different  atoms. 

the  weight  fraction  of  element  A 

the  enhancement  coefficient 

the  atomic  number  correction  factor  where  z  equals  the 
atomic  number 


ABS  *  the  absorption  factor  for  A  radiation  by  the  sample  [exp(-uD)] 
where  u  equals  a  linear  absorption  coefficient  of  the  sample 
for  A  radiation  and  D  is  the  absorption  path  length. 
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I 


For  element  B  a  similar  expression  is  obtained 

IB  =  kbwb  (1  +  Kf)(ABS)  Pz 


(2) 


The  absorption  factor  is  generally  the  roost  troublesome  when  the 

infinitely  thick  approximation  cannot  be  used.  In  order  to  calculate  ABS, 

a  composition  and  path  length  must  be  assumed.  From  Equations  (1),  (2)  and 

(3)  W  and  W  could  Le  calculated  and  these  values  inserted  in  the  ABS  factor. 
*  A  B 


W.  +  Wn 

A  B 


(3) 


Iteration  of  this  process  should  yield  progressively  better  values  of  and 

Wg.  Unfortunately,  the  emergent  path  length  of  the  x-rays  cannot  be  known 

accurately  unless  the  film  thickness  is  determined  with  great  accuracy.  The 

approach  being  attempted  by  the  x-ray  group  assumes  that  the  ABS  factor  is 

nearly  unity;  i.e.,  the  product  uD  is  very  small.  Another  assumption  which 

is  justified  for  elements  of  small  Sz,  such  as  iron  and  nickel  is  that  P  is 

z 

the  3ame  for  both  1^  and  I  and  that  K^  =  Kg. 

Since  NiK  radiation  enhances  FeK  ,  K„  may  be  calculated  for  the  assumed 

r 

Permalloy  composition  as  shown  in  Equation  (4)  below. 


.Nix 


yFeNi  (V  -  Vc  ) 

*  <0'6><W  V  -HI 


^Alloy 


(V  -  V*e) 


Fe\2 


(4) 


U  =  mass  absorption  coefficient  of  Fe  for  NiKa  (400  cm  /g) 

FeNi 

Ni 

u  =  the  mass  absorption  coefficient  of  alloy  for 

°y  NiKa =  radiation  (128  cm^/g). 


Ni 

c 


=  critical  excitation  potential  for  nickel  radiation 
in  kilovolts 


EfeNi  =  excitation  efficiency  of  iron  by  NiKa  radiation  =  0.213. 
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For  the  excitation  voltage  of  25  kv,  F.  is 

2 

Kpl  -  (0.6)  (0.213)  (0.82)  [3.20]  i25  -  8-'3).  -  0.292  (5) 

*5  kv  (25  -  7.1)2 

Combining  Equations  (1) ,  (2)  and  (5)  and  imposing  the  simplifying  assumptions 
discussed  above,  one  obtains 


m  Tp-  (1.292)  (6) 

hi  WNi 

which  can  be  combined  with  Equation  (3)  to  obtain  the  weight  fraction  of 
nickel. 

Ipe 

1  +  (f^)  0.774  (7) 

xNi 

Table  I  presents  some  calculated  nickel  compositions  using  Equation  (7) 
for  a  Permalloy  thin  film  whose  thickness  was  not  known  accurately  but  believed 
to  be  approximately  1,000  angstroms. 


TABLE  I 


^Nickel 

IIron 

... BBStS  n  ! 

651 

108 

0.166 

88.6 

642 

110 

0.171 

88.3 

647 

109 

0.168 

88.4 

658 

110 

0.167 

88.6 

662 

111 

0.168 

88.4 

The  values  of  Table  I  are  considerably  higher  than  the  expected  value  of  82Z 
nickel;  thus  the  assumption  of  the  infinite  thickness  is  probably  not  justified 
for  this  film.  Work  is  currently  in  progress  utilizing  very  thin  (less  than 
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300  angstroms)  films  which  have  been  removed  from  the  glass  substrate  and 
placed  on  a  wire  mesh  holder.  This  latter  step  removes  the  uncertainty 
caused  by  substrate  electron  back  scatter  with  its  resulting  increased 
x-ray  excitation. 
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C.  FACILITIES  FOR  PERMALLOY  THIN  FILM  FABRICATION  AND  EVALUATION 
(J.  K.  Watson) 

During  the  past  six  months  we  have  established  facilities  for  the 

routine  vacuum  deposition  of  planar  permalloy  thin  films  and  for  the 

measurement  with  fair  accuracy  of  thickness,  coercive  force  H  ,  and  aniso- 

c 

tropy  field  A  previously-developed  pulsed-field  system  has  been  adapted 

to  the  measurement  of  dispersion.  Magnetic  domains  can  now  be  observed 
routinely  using  the  Kerr  magnet-optical  effect. 

The  vacuum  deposition  process,  dispersion  measurements,  and  preliminary 
observations  of  magnetic  domains  are  each  discussed  briefly  below.  Metal¬ 
lurgical  measurements  of  alloys  composition  are  discussed  by  Dr.  Gould  else¬ 
where  in  this  report. 

1.  VACUUM  DEPOSITION  PROCESS 

Film  fabrication  is  by  vapor  deposition  in  vacuum  from  an  alloy  which 
is  heated  by  induction  from  a  5-kw  Lepel  RF  generator.  A  Sloan  OMNI  crystal 
oscillator  provides  both  thickness  and  deposition  rate  information,  with 
capability  for  automatic  control  of  these  parameters  in  the  future.  We 
usually  deposit  films  at  the  rate  of  25  angstroms  per  second,  with  our  present 
source  configuration,  Figure  1,  and  with  the  RF  generator  at  about  2/3  of 
maximum  plate  current.  Thus,  there  is  a  reasonable  margin  for  varying  the 
deposition  rate.  The  source  is  nominally  5  grams  and  the  source-substrate 
distance  is  13  Inches.  The  1-cm  diameter  films  are  deposited  through  a  mask 
onto  5/8  inch  square  substrates  of  Corning  7059  pyrex  of  .030  inch  thickness. 
Substrate  temperature  is  controlled  (by  a  Barber  Colman  471-P)  using  a 
standard  500  watt  quartzline  lamp  as  the  heating  element  for  the  0FHC 
copper  heat  sink  which  backs  the  substrates. 
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Molten  permalloy 


Neoted  crucibles 
(99.  7%  A1203) 

-  Crucible  holder 

[Machined  from  Grade  A  Aluminum 

Silicate  (Lava)  ] 

-  Water-cooled 

OFHC  copper  work -coil 


Figure  //I 

Source  Configuration  for  Vacuum  Deposition 


Several  problems  with  this  system  still  prevent  its  operation  as  origi¬ 
nally  intended.  First,  the  background  pressure  in  the  system  rises  by  2  to  3 
orders  of  magnitude  during  deposition,  such  that  the  purity  of  the  films  is 
degraded  by  capture  of  an  estimated  1  to  2  atomic  percent  of  residual  gas 
during  deposition.  Although  this  is  not  an  unreasonable  level  of  contamination, 
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we  seek  an  order  of  magnitude  improvement,  by  the  addition  of  a  Meissner 
trap  and  by  Improving  outgassing  and  deposition  techniques.  Two  tougher 
problems  are  concerned  with  the  melt  confinement:  first,  the  crucible 
usually  breaks  after  only  one  deposition;  and  second,  a  thin  layer  of 
contaminant  sometimes  forms  over  the  surface  of  the  melt,  causing  incon¬ 
sistent  evaporation.  Breakage  has  been  reduced  by  changing  the  crucible 
support  to  reduce  the  thermal  gradients  during  deposition,  and  seems  to  be 
further  reduced  by  a  slow  heat-up  to  the  evaporation  temperature.  Since 
the  surface  contaminant  does  not  always  form  over  the  melt,  it  is  thought  to 
be  affected  by  the  smoothness  of  the  crucible  inner  surface. 

The  surface  contaminant  now  seems  to  be  under  control,  by  prior 
outgassing  of  crucible  at  1500°C,  followed  by  the  addition  of  0.1  percent 
gadolinium  to  the  permalloy  melt. 

Further  plans  call  for  the  continuing  evaluation  of  other  crucible 
materials,  and  for  the  exploration  of  the  feasibility  of  using  a  levitation 
work  coil.  These  developments  are  expected  to  move  gradually,  as  we  are  able 
to  make  more  sophisticated  characterizations  of  film  performance  parameters. 

2.  DISPERSION 

Ue  have  identified  dispersion  (of  the  easy  axis  of  magnetization)  as  a 
key  parameter  for  systematic  investigation.  Together  with  skew,  film  dispersion 
is  the  parameter  which  determines  the  minimum  required  reset  bit-current  in  a 
film  memory  application.  At  high  values,  gross  dispersion  has  been  found  to  be 
affected  by  annealing  procedures,  alloy  composition,  and  film  thickness.  It 
is  difficult  to  achieve  consistent  low  values  of  dispersion,  thus  it  is  a 
sensitive  parameter  of  the  process  variables.  We  have  made  an  arbitrary 
definition  of  six  categories  of  dispersion  in  a  geometric  progression  ranging 
from  "very  low"  (less  than  0.6°)  to  ’very  high"  (greater  than  9.2°)  plus  two 
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additional  categories  for  peculiarities. 

Some  investigators  have  found  a  thickness  dependence  of  dispersion, 

attributed  to  a  magnetization  buckling  caused  by  the  demagnetizing  field.  It 

is  too  early  for  our  investigation  to  be  conclusive  but  our  approach  is: 

In  the  ill-uncerstood  range  of  moderate  to  very  low  dispersion,  we 
attribute  dispersion  to  several  effects  which  can  be  identified  initially 
by  "thickness  dependence"  and  "other".  By  hypothesis,  for  any  thickness 
there  is  a  minimum  value  of  dispersion,  as  suggested  in  Figure  2  we  are 
attempting  to  experimentally  identify  the  components  of  dispersion,  and 
to  establish  their  meaning  theoretically. 

3.  DOMAIN  OBSERVATION  BY  KERR  EFFECT 

Earlier  investigations  have  suggested  that  a  simple  dispersion  measurement 
is  inadequate  as  a  statistical  film  description.  We  have  initiated  a  companion 
study  of  magnetic  domain  configurations  as  a  way  of  gaining  additional  insight 
into  the  influence  of  micromagnetic  structure.  It  is  found  that  residual  domain 
configurations  depend  not  only  on  the  film  properties  but  also  on  the  procedure 
for  demagnetizing  the  film.  For  a  given  domain  pattern,  we  seek  an  identifi¬ 
cation  of  its  key  attributes,  and  to  identify  those  attributes  with  demagne¬ 
tizing  procedures  and/or  appropriate  parameters  of  film  properties. 

Demagnetization  of  the  film  by  a  slowly-decaying  a-c  magnetic  field  has 
been  found  to  cause  strikingly  different  domain  patterns,  as  a  function  of  the 
angle  of  the  applied  field  relative  to  the  film  easy-axis.  Typical  patterns 
are  shown  in  Figure  3  for  easy-axis  and  hard-axis  demagnetization,  where  far 
fewer  domains  are  found  for  the  easy-axis  case.  Preliminary  observations  indi¬ 
cate  an  increase  in  domain  density  with  film  thickness,  with  surprisingly  little 
correlation  with  film  dispersion.  For  high  dispersion  films,  the  optimum 
polarizer/analyzer  orientation  for  best  domain  contrast  is  found  to  vary  over 
the  surface  of  the  film.  For  an  occasional  film,  the  domain  edges  are  not 
straight  as  in  Figure  3  but  are  curved,  due  to  local  skew. 
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Film  Thickness  in  Angstroms 

Figure.  2  Preliminary  Measurements  of  Magnetic  Film  Dispersion  vs  Thickness 


Dispersion  Categories,  (Arbitrary) 

ver^_|^i0w^J_^10rmai^_| - ►moderate  - 1— —  high 

low  * 


a)  AFTER  EASY-AXIS  DEMAGNETIZATION 


b)  AFTER  HARD-AXIS  DEMAGNETIZATION 


FIG,  3  TYPICAL  DOMAINS  IN  A  MAGNETIC  FILM  AFTER  A-C 

DEMAGNETIZATION 
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Demagnetization  of  the  films  by  a  hard-axis  pulse  requires  precise  film 
alignment  and  excellent  homogeneity  of  the  field,  yet  causes  domain  configu¬ 
rations  which  are  not  only  difficult  to  Interpret  but  are  apparently  dependent 
on  the  derivatives  of  the  pulsed  field.  We  have  not  yet  seriously  considered 
these  problems. 

In  summary,  we  have  established  a  facility  for  creating  and  observing 
domain  configurations  in  magnetic  films.  A  simple  procedure  of  a-c  demagneti¬ 
zation  has  helped  identify  some  inhomogeneous  films  with  easy-axis  skew.  We 
are  attempting  to  establish  a  more  quantitative  interpretation  of  observations 
on  typical  films. 
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D.  AN  ANALYSIS  OF  THICKNESS  ERROR  IN  QUANTITATIVE  TRANSMISSION  MICROSCOPY 

OF  CERAMICS  (C.V.  Gokularathnam,  S.W.  Freiman,  R.T.  DeHoff  and  L.L.  Hench) 

1 .  INTRODUCTION 

Optical  microscopy  has  been  a  fundamental  tool  in  the  identification 
of  crystals  for  over  a  century,  and  electron  microscopy  has  been  essential 
for  analysis  of  structural  defects  for  nearly  two  decades.  However,  only 
since  the  development  of  quantitative  stereology  in  recent  years  has  ii  been 
possible  to  quantitatively  analyze  structural  features  such  as  volume  fraction, 
surface  area  of  phases,  number  of  particles  per  unit  volume,  etc.  from  micro¬ 
scopic  measurements.  It  is  well  recognized  that  accurate  quantitative 
analysis  of  defects  and  phase  assemblages  is  essential  if  correlations  between 
structure  and  properties  are  to  be  developed.  It  has  been  shown  that  the 
counting  measurements  and  mathematical  analysis  necessary  for  quantitative 
structural  determinations  are  relatively  simple  and  straightforward  as  long 
as  the  microscopy  is  conducted  on  a  sample  illuminated  in  reflection. 1  The 
surface  of  the  sample  in  such  illumination  is  a  planar  section  of  the  micro¬ 
structure  and  corresponds  to  the  model  employed  in  deriving  the  fundamental 
stereological  parameters,  Pp,  NL,  and  NA<  (The  definitions  of  the  stereology 
terms  and  their  methods  of  measurement  are  reviewed  in  the  appendix.) 

However,  in  many  microscopic  analyses,  it  is  necessary  to  employ 
transmitted  illumination.  For  example,  electron  microscopy  is  limited  to 

O 

transmission  measurements  when  particles  or  defects  are  less  than  100  A 
in  size.  The  optical  transparency  of  many  ceramic  crystals,  glasses  and 
multiphase  bodies  makes  transmission  microscopy  of  thin  sections  necessary. 

The  problem  that  arises  in  transmission  microscopy  is  that  the  thickness 
of  the  sample  must  be  taken  into  account  in  order  for  quantitative  evaluations 
of  structure  to  be  valid.  When  the  sample  thickness  to  particle  diameter 
ratio  (t/b)  is  greater  than  about  2,  a  considerable  error  can  be  introduced 
if  the  effect  of  thickness  is  neglected.  It  is  the  objective  of  the  present 
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paper  to  examine  the  magnitude  of  error  resulting  from  thickness  effects  when 
the  t/b  ratio  is  in  the  range  of  2.  In  many  electron  or  optical  microscopy 
studies  it  is  possible  to  achieve  a  t/b  ratio  of  20  or  more  and  in  these  cases, 
the  potential  thickness  error  is  even  greater  than  that  herein  reported. 

An  additional  objective  of  the  present  paper  is  to  demonstrate  that 
complete  stereological  treatment  of  the  transmission  microscopic  analysis  does 
yield  calculated  values  comparable  to  the  experimental  values  obtained  from 
reflection  measurements  on  the  same  samples. 

The  system  employed  in  this  study  was  a  partially  devitrified  Li20-2Si02 
glass  which  could  be  examined  in  both  reflection  and  transmission  over  a  range 
of  volume  fractions  of  crystals.  These  petrographic  studies  could  also  be 
compared  with  the  detailed  x-ray  analysis  previously  reported  on  the  same 
samples.2 


2 .  THEORY 

In  reflection  microscopy,  there  exist  three  fundamental  stereological 
relationships  that  permit  the  estimation  of  geometric  properties  of  a  three- 
dimensional  structure  from  simple  counting  measurements  made  upon  a  random 
plane  section:  1 


The  quantities  on  the  left  of  these  equations  are  the  counting  measurements, 
which  are  performed  under  the  microscope,  or  on  photomicrographs;  they  are 
fully  explained  in  the  Appendix.  The  quantities  on  the  right  are  properties 
of  the  structure  which  may  be  estimated  from  these  measurements.  is  the 
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volume  fraction  of  the  phase  of  Interest  In  the  sample.  Is  the  total  area 
of  Interface  between  the  matrix  and  the  phase  of  Interest,  In  unit  volume  of 
sample.  Is  the  number  of  particles  in  unit  volume,  and  their  average 
size  (caliper  diameter).  These  relationships  are  valid  for  structures  of 
arbitrary  geometry,  i.e.,  for  aggregates  of  particles  of  arbitrary  shape  and 
size  distribution. 

When  a  sample  is  observed  in  transmission  microscopy,  the  image  contains 
a  projection  of  all  of  the  structure  contained  within  the  volume  of  the  thin 
section  under  study.  Particles  that  intersect  the  top  and  bottom  surface  of 
the  section  contribute  to  the  projected  image  in  proportions  that  are  essentially 
equal  to  those  that  would  be  observed  in  reflection  microscopy.  Two  additional 
kinds  of  contributions  to  the  projected  image  must  also  be  accounted  for;  these 
alter  and  complicate  the  relationships  between  the  counting  measurements  and 
the  properties  of  the  three-dimensional  structure  in  the  thin  section. 

Particles  whose  centers  lie  within  the  volume  of  the  section  contribute 
projected  images  which  must  be  added  to  those  that  result  from  particles  that 
intersect  the  surfaces  of  the  section.  The  equations  in  set  (I)  become:1’3 


(II) 


P  '  -  V„  +  N  tA 

P  V  V  proj 

N/  -  7  S„  +  -  N„tL 
L  2  V  7T  V  proj 


Na'  -  N„D  +  N„t 
A  V  V  V 


(4) 

(5) 

(6) 


where  Pp' ,  N^',  and  N^'  are  the  values  of  the  counting  measurements  described 
in  the  Appendix  measured  in  transmission  microscopy,  A  is  the  average 
projected  area  of  particles  in  the  structure,  L  tlie  average  perimeter 

of  the  projections  of  all  particles  in  the  structure,  and  t  is  the  section  thick¬ 
ness.  Thus,  if  the  equations  in  set  (I)  are  used  to  estimate  properties  of  the 
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structure,  the  values  of  these  properties  will  tend  to  be  overestimated  by 
amounts  corresponding  to  the  second  term  in  each  of  these  equations. 

If  the  section  thickness  is  large  compared  to  the  particle  size,  the 
projections  of  some  particles  will  be  partially  or  wholly  obscured  by  the  projec 
tlons  of  other  particles  in  the  system.  This  effect  of  particle  overlap, 
called  the  Holmes  Effect,3-4  will  tend  to  yield  an  underestimate  of  the 
geometric  parameters  of  the  structure.  The  Holmes  effect  increases  as  the 
section  size  Increases,  and  as  the  volume  fraction  of  the  second  phase  Increases 
The  Holmes  effect  has  been  neglected  in  the  present  study,  because  both  the 
volume  fraction  and  section  thickness  are  small. 

If  it  is  assumed  that  all  of  the  particles  in  the  structure  are  the  same 
size  and  shape,  the  geometric  properties  of  the  structure  may  be  written  in 
terms  of  the  particle  size,  b,  and  shape  factors: 


Dv  - 

kib 

(7) 

sv  - 

NyS  -  Nyk2b2 

(8) 

■ 

> 

> 

V  -  Nvk3b3 

(9) 

Lproj 

-  k4b 

(10) 

*proj 

-  k5b2 

(ID 

where  the  constants,  kj,  k2,  k3,  k4 ,  and  ks  are  shape  factors  which  must  be 
evaluated  for  given  particle  shape  and  choice  of  size  parameter.  Further,  it 
can  be  shown,  that  for  convex  particle  shapes. 


k-4 

-  TTkj 

(12) 

*5 

- 

(13) 
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Substitution  of  these  results  into  equations  (II)  gives: 


V 

-  Nvb2  [k3b  +  ±  k2t] 

(14) 

(IV) 

V 

-  Nyb  [|  k2b  +  2k i t ] 

(15) 

V 

-  Ny  [k jb  +  t] 

(16) 

The  experimental  determination  of  the  three  counting  measurements  on  the 
projected  image,  coupled  with  the  evaluation  of  the  three  shape  factors,  thus 
yields  three  simultaneous  equations  which  may  be  solved  for  the  three  unknowns, 
Ny,  b,  and  t.  The  insertion  of  the  number  of  particles  and  their  size  into 
equations  (8)  and  (9)  then  yields  a  complete  description  of  the  structure  in 
terms  of  its  extensive  geometric  properties:  volume,  surface  area,  number, 
and  size  of  particles.  In  practice,  it  is  necessary  to  obtain  an  independent 
measurement  of  sample  thickness  (t)  in  order  to  solve  the  equations  in  Set  IV 
with  accuracy. 


3.  EXPERIMENTAL  PROCEDURE 

The  partially  crystalline  Li20-2Si02  glass  used  in  this  experiment  was 
prepared  according  to  the  procedure  outlined  by  Freiman  and  Hench.2  The  heat 
treatment  shcedule  consisted  of  a  3  hr  nucleation  treatment  at  475°C  followed 
by  a  crystallization  treatment  of  1  minute  at  595°C  for  sample  1,  13  minutes 
at  592°C  for  sample  2,  and  140  minutes  at  333°C  for  sample  3.  A  petrographic 
thin  section,  20+3  microns  thick,  was  prepared  from  each  sample.  The 
remainder  of  each  sample  was  retained  for  use  in  the  reflection  studies. 

The  thin  sections  were  studied  in  transmission  in  a  Bausch  and  Lomb 
model  LI  dynoptic  polarizing  microscope.  Due  to  the  birefringence  of  the 
devitrified  phase,  the  Li2Si205  crystals  appeared  light  against  a  dark  back¬ 
ground  when  the  section  was  viewed  under  crossed  nicols.  A  typical  micro¬ 
structure  is  shown  in  Figure  1. 
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The  quantitative  stereological  parameters,  PD',  N  *  and  N  '  ,  were 
measured  in  transmission  by  placing  an  8  x  8  grid  in  the  eyepiece  of  the  micro¬ 
scope.  Twenty-five  measurements  of  each  parameter  were  made  on  each  section, 
and  their  average  was  calculated.  At  any  given  position  of  the  microscope 
stage,  a  constant  fraction  of  L^S^Os  crystals  was  invisible  under  crossed 
nicols  because  of  partial  extinction.  The  quantitatve  measurements  were 
corrected  for  this  effect  by  multiplying  the  calculated  average  values  of  Pp'  , 

N^'  and  N^'  by  a  factor  which  accounted  for  the  invisible  fraction  of  the 
particles.  By  measuring  the  angular  extinction  of  100  particles  this  factor 
was  determined  to  be  1.2 

Specimens  were  prepared  for  measurements  in  reflection  by  polishing 
through  0.3  micron  alumina  grit,  and  etching  for  10  minutes  in  2.5%  HF. 5  A 
typical  microstructure  after  such  a  treatment  is  shown  in  Figure  2.  An  average 
of  the  twenty-five  measurements  of  Pp,  and  made  in  a  Bausch  and  Lomb 
Research  Metallograph  was  calculated.  The  thickness  of  each  thin  section  was 
also  measured  in  the  metallograph.  A  typical  cross-sectional  view  of  a  thin 
section  is  shown  in  Figure  3. 

The  evaluation  of  the  average  particle  size  and  number  of  particles  in 
the  system  requires  a  determination  of  the  three  shape  factors,  kj,  k£ ,  and 
k2-  This  requires  an  assumption  or  determination  of  the  shape  of  the  particles 
involved.  Examination  of  the  microstructure  of  the  system  under  study  suggested 
that  the  particles  could  be  modeled  as  ellipsoids  of  revolution.  The  shape 
of  particles  is  specified  in  terms  of  the  axial  ratio,  q,  of  the  generating 
ellipse,  where  q  is  the  ratio  of  the  minor  to  major  axis  (b/a).  Values  of 
the  three  shape  factors  have  been  calculated  as  a  function  of  q,6  and  are 
plotted  in  Figure  4. 

It  has  been  shown  that  the  average  axial  ratio  of  sections  through 
ellipsoids,  F,  is  related  to  the  value  of  q  for  the  structure;7  this  relationship 
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is  plotted  in  Figure  5.  The  value  of  F  was  determined  for  each  sample  by 
measuring  the  axial  ratio  of  100  particle  sections  chosen  at  random.  The  axial 
ratio  of  the  three-dimensional  particles,  q,  was  then  determined  from  Figure  5, 
and  the  corresponding  shape  factors  read  from  Figure  A. 

A .  RESULTS 

The  results  of  the  present  investigation  are  summarized  in  Tables  I  and 
II.  In  Table  I  the  values  of  the  parameters  measured  in  transmission  are  given 
in  the  first  column  for  each  sample,  the  calculated  reflection  values  in  the 
second  column,  and  th _>  values  measured  in  reflection  in  the  third  column.  The 
fourth  column  lists  the  "thickness  error",  the  error  which  would  be  encountered 
if  the  values  measured  in  transmission  were  not  corrected  for  the  thickness 
of  the  thin  section.  The  last  column  gives  the  error  after  correction,  the 
difference  between  the  calculated  values  and  the  values  measured  in  reflection. 

The  "correction  error"  can  be  attributed  to  a  number  of  causes.  There 
was  approximately  A%  error  in  determining  Pp,  and  N^,  resulting  from  human 
subjectivity  in  counting  and  from  reproducing  counting  measurements.  There 
was  probably  some  error  due  to  the  Holmes  effect,  especially  for  the  sample 
of  higher  volume  fraction.  Rather  than  the  particle  size  being  constant,  as 
assumed  in  the  theoretical  analysis,  there  was  approximately  a  3%  dispersion 
in  particle  size,  producing  an  error  in  the  transmission  correction  factor. 
Finally,  there  appeared  to  be  a  small  variation  in  section  thickness  which 
will  influence  transmission  values. 

Table  II  summarizes  the  complete  microscopic  analysis  of  the  devitrified 
glass  samples  measured  in  transmission  and  corrected  for  thickness  effects. 

The  volume  fraction  of  lithium  disilicate  spherulites,  the  number  of 
spherulites  per  unit  volume,  spherulite  size,  the  surface  area  of  the  spherulites. 
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TABLE  I 


SUMMARY  OF  THICKNESS  ERROR  ANALYSIS 


Thickness  Error 


Transmission  -  (Ref lection) ca^ 
(Reflection) ca^ 


Correction  Error 


(Reflection)  ,  -  (Reflection) 

cal _ me  as 

(Reflection) 

meas 


Sample 

Number 

Trans¬ 

mission 

Reflection 

Calculated 

Reflection 

Measured 

Thickness 
Error  % 

Correction 
Error  % 

1 

0.090 

0.023 

PP 

0.020 

291 

15.0 

2 

0.106 

0.025 

0.029 

324 

14.0 

3 

0.213 

0.045 

0.053 

373 

15.1 

1 

177 

103 

nl 

107 

71 

3 

2 

154 

86 

85 

79 

1 

3 

232 

135 

106 

72 

27 

1 

9.3  x  104 

5.0  x  104 

na 

4.8  x  104 

66 

4 

2 

7.0  x  10  4 

4.1  x  104 

4.1  x  104 

71 

0 

3 

8.2  x  104 

3.4  x  104 

4.6  x  104 

141 

26 
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and  Che  average  distance  between  the  particles  is  presented  for  three  samples. 

The  average  distance  between  particles  is  obtained  in  the  relation  A  =  2[1-P,J/N 

i  L> 

It  can  be  seen  that  although  the  same  volume  fraction  of  crystals  can  always 
be  obtained  from  different  heat  treatment  schedules,  the  other  structural 
parameters  may  vary  depending  on  nucleation  treatments  employed,  differences 
in  glass  composition,  experimental  scatter  or  other  factors. 

5.  DISCUSSION  OF  RESULTS 

It  is  clear  from  Table  I  that  the  quantitative  stereological  parameters 
measured  from  transmission  specimens  differ  considerably  from  the  true 
parameters.  This  is  demonstrated  by  the  magnitude  of  the  "thickness  error". 

The  measured  and  calculated  reflection  values  agree  very  closely,  however, 
justifying  the  assumptions  made  in  the  analysis  regarding  particle  size  and 
shape.  Comparison  of  the  values  of  determined  in  this  study  with  those 
expected  from  the  crystallization  kinetics  determined  by  Freiman  and  Hench2 
from  x-ray  diffraction  data  also  confirms  the  accuracy  of  the  present  analysis. 

Examination  of  Table  II  reveals  that  more  insight  into  the  behavior  of 
the  devitrified  glass  samples  is  possible  when  all  the  structural  parameters 
are  known  than  when  volume  fraction  alone  is  measured.  For  example,  Hasselman 
and  Fulrath8  have  shown  the  strength  of  glasses  containing  second  phase 
oxide  particles  is  dependent  on  the  mean  free  path  between  particles.  Two 
of  the  authors  have  recently  demonstrated  the  same  dependence  for  partially 
devitrified  glasses.9  The  porosity  dependence  of  strength  of  polycrys calline 
ceramics  is  also  probably  related  to  the  mean  free  path  between  pores.  As 
shown  in  Table  II,  measurement  of  average  mean  free  path  between  particles 
can  be  determined  analytically  ra’.her  than  assumed  from  volume  fraction 
measurements . 
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Surface  area  per  volume  Is  another  parameter  which  can  be  determined 
analytically  and  is  often  overlooked  in  the  evaluation  of  the  behavior  or 
processing  of  ceramics.  DeHoff  and  Rhines  have  recently  shown  that  the  sequence 
of  geometric  changes  occurring  in  the  densif ication  of  a  sinterbody  are 
directly  related  to  the  change  in  surface  area  of  the  porosity.10  Surface 
area  analysis  is  also  important  in  the  analysis  of  catlysis,  permeability, 
and  sorption-desorption  processes. 

Quantitative  measurements  of  particle  size  and  number  of  particles  per 
unit  volume  are  particularly  useful  in  evaluating  process  variables.  Pore  and 
particle  size  distributions  that  result  from  density  gradients  in  forming, 
or  temperature  gradients  in  firing,  can  be  determined  from  fired  ware  in 
this  manner. 


6.  CONCLUSIONS 

The  results  of  this  investigation  show  that  quantitative  stereological 
parameters  can  be  overestimated  by  as  much  as  300Z  when  measured  in  trans¬ 
mission  If  corrections  for  thickness  effects  are  not  made.  It  has  also  been 
shown  that  structural  parameters  corrected  for  thickness  effects  correspond 
with  values  obtained  from  reflection  measurements  and  x-ray  analysis  within 
about  10%.  Consequently,  it  is  possible  to  obtain  quantitative  correlations 
of  physical  properties  and  process  variables  with  ceramic  microstructures 
from  either  reflection  or  transmission  microscopy  if  the  aoove  analysis  is 
followed. 
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APPENDIX 


NOMENCLATURE  AND  DEFINITIONS  OF  QUANTITIES  MEASURED  IN  THE  MICRO  STRUCTURE 

1.  The  point  fraction,  Pp  (Figure  Al) .  A  grid  is  superimposed  on  the 
structure  and  the  number  of  points  lying  within  the  phase  of  interest 
is  counted.  Pp  is  this  number  divided  by  the  total  number  of  points 
in  the  grid. 

2.  The  number  of  intersections  per  unit  length  of  test  line,  N^  (Figure  A2) . 
A  test  line  is  superimposed  on  the  structure.  The  number  of  times  this 
.'■.ine  intersects  particle  outline,  divided  by  the  total  length  of  test 
line  on  which  measurements  are  made  is  the  intercept  count,  N  . 

l_i 

3.  The  particle  count,  N^  (Figure  A3).  This  parameter  is  defined  as  the 
total  number  of  particles  observed  in  the  structure,  divided  by  the 
area  of  observation.  Particles  intersecting  the  edge  of  the  area  of 
observation  are  counted  as  1/2  particle. 
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Figures  Al,  A2  and  A3.  Illustration  of  the  method  of  measurem* 


IV.  Semiconductors  and  Semiconductor  Devices  (J .  A.  Lindholm,  S.  S.  Li, 
A.  J.  Brodersen,  E.  R.  Chenette,  J.  Salz,  R.  W.  Gould,  J.  J.  Hren, 
L.  L.  Hench,  C.  T.  Sah) 


A.  THE  DENSITY  OF  STATES  IN  DEGENERATE  N-TYPE  GERMANIUM  AND  SILICON 
(S.  S.  Li  and  J.  Agraz-G) 

In  a  nondegenerate  semiconductor,  the  conventional  expression  for  the 


density  of  states  in  the  conduction  band  is  given  by: 

*  .  3/2 


8o(E) 


Att(2M  d) 
h^ 


(E  -  E  )' 
c 


Equation  [1]  predicts  a  zero  density  of  states  at  the  edge  of  conduction  band 
(i.e.,  at  E  ■  Ec)«  However,  as  a  result  of  the  interaction  of  conduction  elec¬ 
trons  with  donor  impurities,  the  density  of  states  in  a  heavily-doped  semi¬ 


conductor  does  not  become  zero  at  the  bottom  of  the  conduction  band:  the 


region  of  permitted  energies  is  shifted  downward.  Sharp  band  edges  disappear 
and  a  tail  in  the  density  of  states,  extending  directly  into  the  forbidden 


gap,  appears. 

Recent  experiments  on  tunneling  (1),  optical  absorption  (2),  and  lumi¬ 
nescence  (3)  have  yielded  experimental  evidence  for  a  "tail"  in  the  density 
of  states  associated  with  impurities  in  degenerate  semiconductors.  In  addition, 
a  considerable  amount  of  theoretical  work  has  been  done  in  determining  the 
effect  of  a  high  density  of  impurities  on  the  density  of  states  (A).  None  of 
these  theories  has  been  used  extensively  to  the  study  of  the  effects  of  doping 
on  the  density  of  states  in  either  germanium  or  silicon. 

The  present  work  deals  with  theoretical  estimates  of  the  density  of  states 

as  a  function  of  impurity  concentrations  for  degenerate  n-type  silicon  and 

(A) 

germanium.  The  study  of  Halperin  and  Lax  (H-L)  gives  the  density  of  states 
deep  in  the  forbidden  gap  and  the  results  of  Bonch-Bruyevich  (B-B) ^  gives 
the  density  of  states  continuously  for  energies  slightly  below  and  above  the 
conduction  band  edge  (i.e.  E  >  Ec)  and  approaches  asymptotically  the  parabolic 
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approximation  (see  Eqn.  [1])  for  energies  exceeding  the  Fermi  level. 


The  H-L  (4)  approach  predicts  the  density  of  states  deep  inside  the 


forbidden  gap  as  given  below: 


where 


g(E)  =  [ (Q  Eg)3/?2]  a (v)  exp  [-E2  b(v)/2^]  [2] 


r  =  —  —  NT  •  —  is  the  screening  length,  N,  is 

^  0  2  I  ’  0  6  ’  ! 


the  impurity  concentration.  e 


2M 


denotes  dielectric  constant. 

M*  =  3(i-  +  —  )-^  is  the  elec¬ 
tron  effective  mass. 

Q2  =  [4Md  e2/e1l  ]  (3^/*)  , 

is  the  density  of  states  effec¬ 


tive  mass. 

v  =  (E^  -  E)/Eq  is  the  energy  below  mean  poten¬ 

tial  E  in  units  of 

°  4ire2  Nj 

EQ  >  Eo  “  "  “ 

a(v)  and  b(v)  are  given  in  reference  (4). 

We  have  made  use  of  Equation  [2]  to  estimate  the  density  of  states  deep  in 
the  forbidden  gap  for  degenerate  N-type  Germanium  and  Silicon  for  different 
doping  concentrations.  The  results  are  displayed  in  Fig.  1  and  Fig.  2. 

For  energies  slightly  below  and  above  the  conduction  band  edge,  the  density 
of  states  g(E)  is  given  by  B  -  B  (5): 

(2M*  E)3^2 

8(E)  "  ~372  exp("x)  [I1/4<X>  +  I-l/4<x>  +  1 3/ 4 (x)  +  I_3/4<x)]  [3] 

E2 

where  _  „04  „  _  — *• 


a2  =  •  Nj,  and  x  =  ^ 

€  Q 

The  Ir(x)  are  Bessel  functions  with  imaginary  a-.guments. 


In  the  energy  region  E  ■  Ep  (Fermi  level),  Equation  [3]  reduces  to: 

g(E)  =  g  (E)  (1  -  -4  ) 

E2 


14] 


Near  the  bottom  of  the  conduction  band  (i.e.  I  E 


0) ,  the  density  of 


states  changes  radically.  In  this  region,  Equation  [3]  reduces  to: 


■»/ it  a  1  .  , 

g(E)  “  g0<E>  i  E  r t  1 


4r(|) 


where  T (r)  denotes  the  Gamma  function. 

Equations  [2],  [4]  and  [5]  provide  a  direct  means  for  displaying  the 
density  of  states  g(E)  as  a  function  of  energy  E  for  arbitrary  impurity  concen¬ 
trations.  Fig.  1  and  Fig.  2  demonstrate  g(E)  versus  E  for  N-type  germanium 

19  21  -3 

and  silicon  with  donor  concentrations  varying  from  10  to  10  cm  .  Note  that 

in  the  intermediate  region  where  the  two  approaches  (i.e.  H-L  and  B-B  approaches) 

fail,  certain  approximations  have  been  made  to  smooth  out  the  g(E)  curves. 

Fig.  1  shows  that  in  an  N-type  germanium  with  lO^cm  ^  donor  impurities, 

the  band  edge  shifted  downward  by  0.022  ev  from  the  unperturbed  conduction  band 

edge;  the  density  of  states  in  the  tail  decreases  rapidly  with  depth  into  the 

forbidden  bands  (faster  than  exponentially).  As  a  result  of  screening,  the 

number  of  states  in  the  forbidden  band  turn  out  to  be  small  compared  with  the 

number  of  states  inside  the  conduction  band.  Shepherd  et  al  (1)  have  observed 

19  i 

an  exponential  band  edge  in  degenerate  N-type  germanium  doping  of  1.2  x  10  cm  J  Sb. 
This  gives  reasonable  agreement  with  the  present  results. 

Deep  in  the  conduction  band  (in  particular,  near  the  Fermi-level  of  N-type 
uncompensated  semiconductors),  the  density  of  states  in  degenerate  semiconductors 
differs  only  slightly  from  the  conventional  parabolic  approximations  (See  Fig.  1 
and  Fig.  2)  . 
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The  above  approaches  yield  single-band  approximations  and  are  valid 
only  for  the  uncompensated  or  weakly-compensated  degenerate  materials. 
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FIGURE  CAPTIONS 


Density  of  states  g(E)  vs  energy  E(ev)  for  N-type  germanium 
for  Nq  =  1019,  5  x  1019,  lO^^cm  ^  respectively. 

♦  denotes  the  lower  limit  of  B-B  approach  (5) 

x  denotes  the  upper  limit  of  H-L  approach  (4) 

Density  of  States  g(E)  vs  energy  E(ev)  for  N-type  silicon 

1 Q  on  oi  -3 

for  Nq  =  5  x  10  ,  5  x  10  ,  10  cm  respectively. 

4  denotes  the  lower  limit  of  B-B  approach  (5) 

x  denotes  the  upper  limit  of  H-L  approach  (4) 
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B.  EXTENSIONS  OF  THE  THEORY  OF  PN  JUNCTIONS:  PART  I  (D.  D.  Kleppinger, 
F.  A.  Lindholm) 


I.  INTRODUCTION 


This  interim  report  extends  the  conventional  equations  for  non-degen¬ 
erate  pn  step  junctions  by  removing  the  Boltzmann  approximation  for  carrier 
densities.  The  major  approximation  underlying  the  present  analysis  is  that 
the  density  of  states  is  not  a  function  of  impurity  concentration.  The 
validity  of  this  approximation  is  discussed  in  Section  VII. 

The  problem  of  describing  junctions  in  material  for  which  Fermi-Dirac 
statistics  must  be  used  is  important  since  most  practical  devices  are  made 
with  such  material.  The  emitter-base  diode  of  a  diffused,  planar  transistor 
is  an  example. 

The  works  of  Bonch-Bruyevich  [l]  and  Ivanchik  [2]  are  examples  of  the 
theoretical  treatments  of  degenerate  material  and  junctions  that  have  appeared 
in  the  literature.  The  results  that  they  report  are  unsuited  for  direct  use 
in  engineering  design. 

The  brief  review  of  some  aspects  of  semiconductor  theory  that  follows  is 
intended  to  provide  a  basis  for  the  remainder  of  the  paper  and  to  introduce 
some  of  the  approximations  to  be  used. 

The  assumption  of  non-degeneracy,  which  is  employed  in  conventional  junc¬ 
tion  theory,  permits  the  use  of  Maxwel 1 -Bol tzmann  statistics, 


f(E)  =  exp  ), 

instead  of  the  Fermi-Dirac  statistics, 

f(E)  *  [1  +  exp  (Lil)]"1 


(1) 

(2) 


This  replacement  [3,4,5]  makes  the  resulting  equations  tractable  in  junction 
design  and  is  appropriate  when  the  Fermi-level  Ep  is  several  kT  units  of  energy 
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away  from  the  energy  of  interest,  E.  For  intrinsic  germanium  and  silicon, 
the  Fermi-level  lies  near  the  middle  of  the  forbidden  energy  gap  and  thereby 
reflects  the  non-degeneracy  of  the  material.  As  impurities  are  added  to  the 
material  the  Fermi-level  approaches  one  of  the  band  edges  (k]  and,  at  great 
enough  impurity  concentrations,  enters  the  allowed  band.  Since  in  junction 
analysis  the  energy  of  interest  which  determines  the  use  or  non-use  of  Boltzmann's 
statistics  is  the  band  edge,  the  Fermi-Dirac  statistics  must  be  used  at  high 
impurity  concentrations. 

The  equation  for  the  conduction  electron  density,  for  example  is  given  by 

(6)  r" 

v  *  5  3/2  1/2 

\  4i»(2mc/h  J  (E-Ec)  dE  (3) 

°  JZC  I  +  exp[ (E  Er)/kT] 

where  (I)  and  the  density  of  states  function 

*  2  3/2  1/2 

g(E)dE  =  l*jr(2mc/h  )  (E-Eq)  (M 

have  been  used.  Two  assumptions  are  made  in  the  derivation  of  ( 3 )  *  The  first 
is  that  the  upper  limit  of  integration  can  be  justified  by  arguing  that  f ( E ) 
decreases  more  rapidly  than  g(E)  increases.  The  dependence  of  the  density  of 
states  Function  on  the  square  root  of  energy  is  also  an  assumption  that  is  valid 
near  the  conduction  band  edge.  The  basis  for  this  assumption  also  allows  us  to 
use  the  scalar  effective  mass.  Some  further  comments  on  this  assumption  are 
made  in  Section  VII.  An  equation  analogous  to  (3), 

jv  *  2  3/2  i/2 

V  W2m  /h  )  (EyE)  dE  , 

P0  -  \  - * - - ^ -  (5) 

I  ♦  exp  [ (Ep-E) /kT] 

holds  for  holes  and  is  subject  to  similar  assumptions.  Equation  (3)  has  the 


-  NC^'2l£f  -  Ei> 
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where 


N 


C 


*  3/2 

2(2*1^  kT) 

h3 


(7) 


and 


^a  -1/2  \  -I 

*  <!*-!*)  -  2n  J  (E  *-E^)  (1  ♦  exp(E '  -Ep)  ]  d(E'-E’)  (8) 

o 


is  a  Fermi-Dirac  integral  [6].  We  will  assume  (6)  correctly  describes  the 
electron  density  but  we  will  discuss,  in  Section  VII,  the  approximation  for 
g(E)  and  its  implications. 

Quasi-Fermi  levels  are  used  in  the  development  to  describe  the  densities 
of  carriers  in  non-equilibrium  [4,7] •  Sah  has  shown  [8]  that  the  quasi-Fermi 
levels  in  space  charge  regions,  for  the  special  case  of  a  symmetric  step 
junction  ,are  constant  to  a  good  approximation  in  material  doped  at  high  impurity 
concentrations.  In  our  analysis,  we  will  assume  that  the  quasi-Fermi  levels 
are  constant  in  junction  regions. 

The  remainder  of  this  report  proceeds  in  the  following  order.  In  Section  II 
we  introduce  an  approximation  for  the  Fermi-Oirac  integral  and  use  it  to  calculate 
the  np  product  in  equilibrium  and  non-equilibrium.  Section  III  describes  the 
development  of  equations  for  the  contact  potential  of  a  junction  and  the  Fermi- 
level-  An  equation  for  the  excess  minority  carrier  density  at  the  edge  of  the 
space  charge  region  is  developed  in  Section  IV.  Section  V  presents  c  quantita¬ 
tive  definition  of  the  voltage  limit  for  low  injection  in  terms  of  impurity 
concentration.  In  Section  VI  we  derive  an  equation  for  the  current  density  and 
show,  through  the  use  of  an  extension  of  the  Einstein  relationship,  that  it  is 
identical  to  the  non-degenerate  equation.  Section  VII  discusses  the  results  and 
the  assumptions  made  in  obtaining  them.  Each  of  the  new  equations  reduces  to 
or  is  Identical  with  its  non-degenerate  counterpart  in  the  limit  of  low  impurity 
concentration . 

Illustrative  calculations  are  made  for  silicon.  The  values  of  the  material 
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parameters  and  various  constants  used  in  the  calculations  are  listed  in 
Appendix  B.  A  ..able  of  symbols  used  throughout  the  paper  is  also  included, 


as  Appendix  A. 


2.  THE  NP  PRODUCT 


The  Fermi  integral  (8)  cannot  be  integrated  analytically  but  it  can  be 


approximated  by  [6] 


3  (Ep  -  Ef) 

i  n  F  b  ' 


exp (Ep  ‘  Eq) 

1  +  0.27  exp (E^  -  E£) 


which  is  a  tractable  function  for  engineering  design  of  junctions.  This  approxi 
mation  is  accurate  to  within  +  3%  for  (Ep  -  E^)  —  1.3  and  is  applicable  to 
impurity  concentrations  of  interest  in  this  paper. 

The  equilibrium  density  of  conduction  electrons  is,  from  (6)  and  (9), 


no  “  NC 


The  hole  dens i ty  is 


p  *  N 

o  V 


exp(Ep  -  Ej!) 

I  +  0.27  exp  (Ep  -  E^) 


exp (Ey  -  E,i) 

I  +  0.27  exp  (Ey  -  Ep) 


in  thermal  equilibrium.  For  non-equilibrium,  the  carrier  densities  can  be 
expressed  in  terms  of  quasi-Fermi  levels  ( 4 , 7 j 


n  =  N 


\  exp(E^e  ~  Ec) 
C  I  +  0.27  exp 


P  =  N. 


The  product  of  (12)  and  (13)  is 


e*p(E^  •  Eph) 

I  +  0.27  exp(Ey  - 


NcNy  exp(-E^)  exp (qV 1 ) 


I  +  1 0 . 27  [exp  (Ep^  -  V)  ♦  exp  ( E^  -  E^)’  +  0.0729  expf-E^)}  exp  (qV’J 


,  (14) 
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The  energy  equivalent  of  the  applied  voltage,  qV 1  ,  is  equal  to  the  difference 
between  the  quasi-Fermi  levels  because  of  our  assumption  that  the  quasi-Fermi 
levels  are  constant.  Equation  (14)  can  be  rewritten  as 


where 


np  = 


2 

n.  exn(qV') 

~~\  ♦  K  exp (qty' } 


(15) 


Nc\  e’<P<‘EG, 


(16) 


The  factor  K  can  be  evaluated  at  the  n-  and  p-type  edges  of  the  space  charge 
region  as  fol lows : 

Kn  -  0.27  [exp (Ep  -  EjL  -  qV')n  ♦  exp(E^  -  E,i)p]  ♦  0.0729  exp(-E^)  (17) 

and 

Kp  -  0.27  [exp  (Ep  -  E^)p  ♦  exp(E^  -  E'  -  qV)p]  +  0.0729  exp(-E(l)  .  (18) 

In  (17)  and  ( 1 8) ,  E^  is  to  be  interpreted  as  the  Fermi-level,  normalized  with 
respect  to  the  thermal  energy  kT,  which  accounts  for  the  equi I ibr ium  carrier 
densities  in  the  n-type  material  for  (17)  and  the  p-type  material  for  (18). 

In  other  words,  it  is  the  quasi-Fermi  level  far  enough  from  the  junction  that 
the  carrier  densities  have  attained  equilibrium  values.  E'  appears  only  with 
respect  to  a  band  edge  in  the  equations  and  each  of  the  terms  are  subscripted 
to  remind  us  to  calculate  the  energy  differences  on  particular  sides  of  the 
junct ion. 

The  equilibrium  (qV‘  =  0)  form  of  (15)  , 


noPo 


n . 1 

i 


I  ♦  0.27  [exp(Ep  -  E£)  ♦  exp(E^  -  Ep)J  ♦  0.0729  exp(-E£) 


(19) 


reduces  to  the  conventional  non-degenerate  form 

2 

n  p  =  n . 
oKo  i 


(20) 
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in  the  limit  of  low  doping.  Equation  (19)  is  valid  in  bulk  material  as  well 
as  at  space  charge  regions. 

Figure  I  shows  the  behavior  of  the  np  product  in  silicon  predicted  by  (15). 
The  product  decreases  from  that  expected  from  the  conventional  expression  as 
the  impurity  concentration  increases.  This  decrease  follows  from  consideration 
of  the  statistics  which  show  that  the  minority  carrier  density  decreases  faster 
with  increasing  doping  than  the  majority  carrier  density  increases.  The  denomi¬ 
nator  of  (15)  is  a  function  of  voltage  as  well  as  impurity  concentration  but 
the  voltage  dependence  does  not  become  important  until  the  applied  voltage  is 
within  a  few  percent  of  the  low  injection  voltage  limit  to  be  defined  in  Section 
V.  Figure  I  is  therefore  applicable  to  low  injection  non-equilibrium  in  addition 
to  equ i I i br i urn . 

Equation  (14)  is  actually  more  general  than  we  have  thus  far  implied  since 
it  can  be  applied  to  any  situation  where  the  differences  (E^  -  E^)  and  (Ey  -  Epfi) 
can  be  evaluated.  In  particular  these  differences  can  be  determined  from  the 
equations  which  describe  the  spatial  variation  of  the  band  edges  throughout  a 
junction  region  and  the  resulting  np  product  can  be  used  to  study  recombination 
in  the  junction  region  as  a  function  of  position  and  applied  voltage. 

3.  THERMAL  tUUlUBRIUM:  THE  FERMI  LEVEL  AND  THE  CONTACT  POTENTIAL 


Equations  for  the  Fermi-levels  may  be  obtained  by  using  charge  neutrality 

along  with  (10)  for  n-type  material  or  (II)  for  p-type  material.  Let  us  consider 

uncompensated  n-type  material  which  is  sufficiently  extrinsic  that  the  number 

of  holes  is  negligible.  Space  charge  neutrality  then  requires  that 

+ 


n  = 


o 


(21) 


where 


+ 


I  +  2  exp ( Ep  -  E£) 


(22) 


The  factor  of  2  in  (22)  accounts  for  the  degeneracy  of  the  impurity  level  [9]. 


-  ' 
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Therefore , 


"C  exp(E'  -  £•) 


1  +  0.27  exp(E '  -  E  ' ) 
F  C 


+  2  exp(Ep  -  Ep) 


(23) 


When  the  substitutions 


a  -  exp(E^),  6C  ■  exp(-E(l),  and  ■  exp(-E^) 


(2k) 


are  made,  (23)  becomes 


NC  aBC 


I  +  0.27  a8, 


I  +  2 


(25) 


The  solut ion  for  a  is 


Nr  -  0.27  N_ 

a  - \L  + 

H  Bo 


<"c  -  "D> 

^  \  V2 


1/2 


N 


2  N  6  B 
CCD 


(26) 


where  the  positive  square  root  is  chosen  so  that  Ep  is  real 


Us  i  rn  the  i dent i ty 


-I 


ln(x  +  ^x2  +  a2)  ■  In  a  +  sinh  *.  , 

a 

we  can  obtain  the  Fermi-level  as 

N  -I  ”1/2  e ’  -  El 

E^  -  y(  E 1  *  E ' )  +I  -  sinh"  ( (Nc  -  0.27  ^>(8^10  exp(-~ — -)  ]  .(27) 

C  D 


Similarly,  in  uncompensated  p-type  material  under  the  assumption  of  negligible 
electron  density,  the  Fermi-level  is 

N,  -1/2 

'NV 


EF  "  7<Ev  +  EA}  •  I  'n(Wu)  +  5inh''  ((NV  -  0.27  NA)  (8NVNA)  exp(LL_iA)  j  (28) 
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If  the  impurity  concentration  is  such  that  0.27  Nq  may  be  neglected  with 
respect  to  N^.  in  (27),  the  equation  reduces  to  the  non-degenerate  expression  [9], 

Ep  -  j  (E £  +  E’)  +  I  ln(^)  -  sinh  (nC/8Nd)  exP<  °  C)]  *  (29: 


When  Np  has  increased  to  a  value  such  that  the  reduction  of  (27)  to  (29)  cannot 
be  allowed  in  good  conscience,  we  may  define  that  level  of  doping  as  the  border¬ 
line  separating  non-degeneracy  and  degeneracy.  This  is  an  arbitrary  definition 
but  (27)  does  permit  a  comparison  directly  between  the  impurity  concentration 
and  the  effective  density  of  states  for  the  particular  material.  The  arbitrari¬ 
ness  of  the  definition  arises  from  the  need  to  choose  how  large  N^,  will  be 
allowed  to  get  with  respect  to  N^.  before  (29)  is  not  a  valid  reduction  of  (27). 

Figure  2  illustrates  the  behavior  predicted  by  (27)  and  (28)  for  phosphorous 
or  boron  doped  silicon.  The  Fermi-level  monotonical ly  approaches  the  band  edge 
as  the  impurity  concentration  increases.  At  equal  impurity  concentrations,  the 
Fermi-levels  are  not  in  the  same  position  relative  to  the  nearest  band  edoe  in 
n-  and  p-type  material  because  of  the  difference  between  the  effective  density 
of  state  values.  A  difference  between  activation  energies  for  the  donors  and 
acceptors  would  also  cause  a  difference  between  the  relative  Fermi-level  posi¬ 
tions.  (See  Appendix  B  for  values  used  in  the  illustrative  calculations). 

Equations  10  and  11  may  be  solved  for  the  Fermi-level: 


in  n-type  material  and 


__  F;_,nfNC-°-27  nno  * 

F  C  1  nno 

fNw  -  0.27 

:*  =  E '  +  In  \J. _ _££_ 

‘F  V  L  pnn 


in  p-type  material.  Equations  (30)  and  (31)  can  be  used  to  calculate  the  Fermi- 
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level  position  as  a  function  of  carrier  density  while  (27)  and  (28)  show  the 
dependence  of  the  Fermi-level  on  the  concentration  and  the  characteristics  of 


the  impurities  through  the  activation  energy.  The  two  sets  of  equations  there¬ 
fore  provide  complementary  information  and  additional  Insight. 

Equations  (10)  and  (II)  can  be  used  together  to  develop  an  equation  for 
the  contact  potential  of  a  pn  junction.  The  definition  for  Ny  is  analogous  to 


(7)  for  Nj.  > 


1/2 

2(2irmy  kT) 


Then  Nc  is  related  to  Ny  by 


*c  3/2 


so  that  (16)  may  be  rewritten 


mr  3A  Er 

4-  n.  «xp<-§-> 
mv  1 


Er 

nj  exp(-|  ) 


The  density  of  electrons  on  the  n-type  side  of  the  junction  is  by  (10) 

"ife  exp(Ei  •  vt ♦  4-  > 


I  +  0.27  exp (Ep  -  E') 


and  the  density  of  holes  on  the  p-type  side  of  the  junction  is  by  (II) 


K  |3A  n 
"i  4  w(E;  •  E^r  > 

'  mC  1 

I  +  0.27  exp (E^  -  Ep) 


Choosing  the  reference  energy  as  Ey  *  0  (therefore  E£  *  Eq)  and  taking  the 
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logarithm  of  these  equations,  we  obtain 
n  E' 

In  -  -  (E  •  •  A  )  -  I  n[  I  ♦  0.27  explt;  -  El  )) 

K  3/<.  F  2  F  c 

n .  — 
i  m“ 
v> 


Ppo  JL 


■E '  +  JL  -  I n [ I  +  0.27  exp (E '  -  E')] 
r  ?  F  C 


The  contact  potential  can  be  obtained  as  the  potential  equivalent  of  the 

energy  difference  between  any  two  identically  defined  energies  on  the  two  sides 

of  the  junction.  In  the  present  case,  the  contact  potential  can  be  obtained 

from  the  difference  in  energy  between  the  mid-gap  energies  ( E ' / 2 )  shown  in 

G 

the  above  equations.  If  (38)  and  (39)  are  rearranged  and  subtracted,  the 
result  is  .  *  *  ,/i. 


E  ' 

.  G 

=  1  n 

2 

P 

n 

no 

*  3/t» 

r  c\ 


ppo  4 

i  m« 


from  which 


+  I n [ I  +  0.27  exp(E '  -  E')]  +  I n ( I  +  0.27  exp(-E')] 

F  G  F 


n  p 
norpo 


n2  +  I  n ( I  +  0.27  exp(Ep  -  E^)n) 


+  I n ( I  +0.27  exp(E '  -  E')  ]  . 

V  F  p 


In  (1*1)  the  reference  energy  has  been  changed  so  that  the  exponentials  explicitly 
show  dependence  on  the  position  of  the  Fermi-level  with  respect  to  the  band  edges. 
The  subscripts  on  the  arguments  of  the  exponentials  are  to  remind  us  exactly 
what  energy  difference  is  required.  Equation  (41)  correctly  predicts  that  the 


potential  of  the  n-type  region  Is  positive  with  respect  to  the  p-type  region. 
(Recall  that  qV^  is  a  normalized  energy  equivalent  of  the  contact  potential  Vg.) 

The  form  of  (1*1)  is  interesting.  The  conventional  non-degenerate  contact 
potential  equation  Is 


I  n 


nnoPPQ 

n? 


(42) 


or  more  usual ly , 

.  ,  N°N* 

where  1 00%  ionization  is  assumed  in  (**3)  •  The  first  term  of  (41)  is  the  same 
as  (1*2)  hence  the  last  two  terms  of  (1*1)  can  be  thought  of  as  correction  terms 
to  the  non-degenerate  equation.  Equation  (1*1)  reduces  to  (**2)  when  the  correction 
terms  are  negligible  and  this  occurs  for  low  doping  concentrations  on  both  sides 
of  the  junction.  The  form  of  (1*1)  makes  it  clear  that  there  will  be  a  correction 
term  even  if  only  one  side  of  the  junction  is  heavily  doped.  The  correction 
terms  are  always  positive  so  (42)  underestimates  the  contact  potential.  Equation 
(1*3)  always  overestimates  the  contact  potential  for  doping  densities  greater  than 
I0*^cm  The  estimates  made  by  using  (42)  or  (43)  become  less  accurate  as 
doping  increases  with  (43)  leading  to  very  considerable  errors,  as  is  discussed 
be  I ow . 

Figure  3  illustrates  the  behavior  predicted  by  ( 4 1 )  for  junctions  in 
silicon  doped  with  phosphorous  and  boron.  The  contact  potential  increases  as 
either  or  both  sides  of  the  junction  become  more  heavily  doped. 

Figure  4  indicates  the  magnitude  of  the  correction  terms  in  (41)  and 
therefore  the  error  which  would  arise  from  using  (42).  The  error  involved  in 
using  (43)  is  also  shown  for  the  single  case  of  equal  doping  on  both  sides  of 
the  junction.  The  errors  do  increase  with  increasing  doping  as  mentioned  above. 
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Figure  5  was  generated  in  the  course  of  the  investigation  and  is 


included  for  reference.  The  curves  clearly  show  that  1 00%  ionization  is  a 
poor  approximation  over  most  of  the  range  of  useful  impurity  densities.  The 
difference  between  the  curves  tor  phosphorous  and  boron  is  due  to  the  difference 
in  the  carrier  effective  masses  and  hence  the  effective  density  of  states,  (7) 
and  (32).  The  ionization  energy  of  the  impurity  will  also  determine  the  position 
of  the  curve.  The  values  shown  in  Figure  5  were  calculated  using  1 )  and  (10) 
or  (II). 

An  alternative  equation  for  the  contact  potential  is 

("“i >  ■  e;  *  (ef  -  E;>„  *  <EJ  ■  efV  m 

Equation  (bb)  can  be  obtained  by  combining  (1*1),  (9),  (10),  and  (16)  or  by 
inspection  from  an  energy  band  diagram  of  a  junction  in  equilibrium.  Figure  (10) 
in  Appendix  A  is  such  a  diagram.  Equations  (bb)  and  (1*1)  are  valid  in  general; 
however,  (1*1)  shows  explicitly  the  correction  terms  while  (kb)  provides  ease  of 
computation  since  both  equations  require  the  same  energy  differences  to  be  known. 

4.  THE  EXCESS  MINORITY  CARRIER  DENSITY  AT  THE  SPACE  CHARGE  REGION  EOGE 

An  equation  which  gives  the  excess  hole  concentration  at  the  edge  of 
the  space  charge  region  at  the  n-type  side  of  the  junction  can  be  obtained  by 
starting  with  (13).  For  a  forward  biased  junction,  the  excess  minority  hole 
density  in  the  n-type  quas i -neu t ra  I  region  achieves  its  largest  value  at  the 
edge  of  the  space  charge  region.  If  the  total  density  of  holes,  equilibrium 
plus  excess,  at  the  edge  of  the  space  charge  region  is  small  compared  to  the 
density  of  majority  electrons  present,  we  define  this  condition  as  low  injection. 
Equation  (13)  can  be  rewritten  as 

N V  exp*EV  "  EF^  cxp(EF  *  EFh)  (1*5) 

Pn  "  I  +  0.27  exp(Ey  -  E<)  exp(E^  -  E£h) 

The  majority  carrier  density  is  assumed  to  deviate  negligibly  from  its  equilibrium 
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I 


IMPURITY  DENSITY  CM 


value. 


This  assumption  implies 


E'  =  E1 
F  Fc 


and  (^5)  may  be  re-expressed  as 


(46) 


Nv  exp  (Ey  -  Ep  exp  (q  V') 

Pn  =  I  7  0.27  exp (Ey  -  E^)  expTqV 1 ) ‘ 


(47) 


The  assumption  of  constant  quasi-Fermi  levels  was  used  as  in  Section  II  to 
relate  the  applied  voltage  to  the  difference  between  the  quasi-Fermi  levels. 
The  tota I  density  of  minority  holes  at  the  space  charge  region  edge  is 


(I  +  0.27  exp(Ey  -  Ep ) ]  [exp(qV')] 

Pf1  "  Pn°  1  +  0.27  exp(E*  -Ep)  exp(qV) 

by  rewriting  (f»7) ;  therefore  the  excess  minor  i  ty  hole  density  is 

[ exn  (qV  1 )  -  I] 

P'  ‘  Pn  ’  Pno  =  I  ♦  0.27  “exp(E(,  -  E  p)  n  exp(q  V '  T  ’ 


(48) 


(49) 


Equations  (48)  and  (^9)  are  based  on  the  assumption  of  low  injection.  The 
result  of  the  development  to  be  made  in  Section  5  shows  that  the  denominator  of 
(f<9)  may  become  different  from  unity  for  certain  combinations  of  impurity  con¬ 
centration  and  applied  voltage  without  violating  the  assumption  of  low  injection 
Thus,  for  these  combinations,  the  excess  minority  carrier  density  will  be  less 
than  the  value  predicted  by  the  non-degenerate  equation 

P|  =  Pno  (exp(qV')  -  I)  (50) 

Notice  further  that  the  value  of  p  in  (49)  should  be  calculated  by  using 

no 

2  2 

(II)  or  (19)  and  not  by  using  pn0  ■  n./n  =  ni/N0  as  wou,d  be  done  for  non- 
degenerate  cases. 

The  excess  density  of  electrons  at  the  edge  of  the  space  charge  region  in 
p-type  materia 

by  analogy  with  (1»9)  . 


I  i  s 


-  n 


npo  (exp(qV')  -  I] 


po 


I  +  0.27  exp(Ei  -  E')  exp(qV) 


(50 


*»**•'' 
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5.  A  QUANTITATIVE  DEFINITION  OF  LOW  INJECTION 


The  qualitative  definition  of  low  injection  given  In  Section  IV  can  be 
made  quantitative.  Low  injection  was  described  as  being  the  condition  existing 
when  the  total  density  of  minority  carriers  is  small  with  respect  to  the  density 
of  majority  carriers  at  the  space  charge  region  edge.  Recall  that  the  minority 
carrier  density  in  the  quas i -neutral  region  will  have  its  greatest  value  at  the 
space  charge  region  edge.  The  definition  can  be  written  as 

?n  ^  F  nno '  (52) 

at  the  edge  of  the  space  charge  region  on  the  n-type  side,  where  F  is  some 
number  less  than  unity.  The  voltage  at  which 

Pn  "  F  nno  (53) 


is  now  defined  as  the  voltage  limit  for  low  injection  and  its  value  may  be  found 


by  substituting  (9)  and  (k7)  into  (53), 

Nv  exp(Ey  -  E')nexp(qV^)  fH  exp(E'  -  E!)n 

.  *  v  r  L  (5 q) 

I  +  0.27  exp (Ey  -  exp(qV^)  |  +  0-27  exp(E,L  -  E£)n 

The  solution  to  (5I1)  is  * 

•>  f mA  \  /m”\3/2 

qV'  »  In  F  ♦  y  In  |CJ  +  (Ep  -  E^  +  E'  -  In  <exp[-(E£  -  E^J*  0.27[I-F  4  ]} 

Vmv  y  n  vV 

( 

where  the  relationships  (Ey  -  Ep)^  -  E£  -  (Ep  -  E^)p  and  (33)  have  been  used. 

For  a  given  material  (which  determines  m"/^y)  the  maximum  voltage  for  which  low 
injection  exists  (as  defined  by  F)  can  be  calculated  as  a  function  of  impurity 
concentration  by  using  (55). 

Figure  6  indicates  the  predictions  of  (55)  as  a  graph  of  qV^  versus  (E^  -  E^)^ 
with  F  as  a  parameter.  For  a  given  F,  the  voltage  required  to  violate  the 

definition  of  low  injection,  (52)  increases  as  the  impurity  concentration  increases. 
This  is  reasonable  since  a  larger  number  of  injected  carriers  will  be  required  to 
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attain  the  given  fraction  of  the  increasing  number  of  majority  carriers.  Also 
reasonable  Is  the  observation  that,  at  a  given  impurity  concentration,  the 
maximum  voltage  allowed  decreases  as  the  fraction,  F,  decreases  since  less 
voltage  is  required  to  cause  the  injection  of  the  smaller  density  of  carriers. 
The  abscissa  is  plotted  as  (E^  -  E^)  to  make  the  graph  applicable  to  silicon  in 
qcneral.  The  graph  may  lx-  made  specific  for  phosphorous  doped  silicon  by  using 
Figure  2  to  obtain  (E^  -  Eji)  for  .1  qiven  impurity  concentration. 

The  equation  analogous  to  (551  for  injection  into  p-type  material  is 


In  (cxpHE'-E')  ]  +  0.27 

V  r  p 


[l-F 


mv  3/2 


(56) 


The  voltage  limit  for  low  injection  can  be  compared  to  the  contact  poten¬ 
tial  by  using  (Al)  and  (55).  Figure  7  illustrates  this  comparison  for  the  F  =  0.1 

by  displaying  the  ratio  qV/qV'  vermis  impurity  concentration  in  phosphorous 

L  D 

doped  silicon.  If  ihe  low  injection  limit  voltage  exceeds  the  contact  potential, 
low  injection  cannot  be  violated  since  a  voltage  greater  than  the  contact  poten¬ 
tial  cannot  be  applied  to  the  junction.  The  curves  show  that  the  low  injection 
criteria  can  be  violated  on  the  more  lightly  doped  side  of  the  junction  but  cannot 
be  violated  on  the  more  heavily  doped  side.  The  qualitative  results  arc  not 
unexpected  but  the  quantitative  results  are  interesting. 

Equation  (55)  may  be  used  to  study  the  behavior  of  the  denominator  of  (f'°)  , 
the  equation  for  the  excess  minority  carr ior (hole)  density  at  the  edge  mr  the 
space  charge  region.  The  questions  to  be  answered  are:  (I)  Docs  the  denominator 
differ  from  unity  before  the  low  injection  assumption  is  violated?  and  (2)  Poes 
the  denominator  differ  from  unity  because  of  its  dependency  on  impurity  concen- 
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tration?  The  answers  developed  below,  to  these  questions  show  that  (^9)  does 

not  reduce  to  the  non-degenerate  form  (50)  in  general  because  the  denominator 

does  differ  from  unity  under  certain  conditions. 

If  the  applied  voltage  is  cqV 1  ,  where  c  5  I,  the  denominator  of  (1*9)  becomes 

L 

Denom.  •  I  ♦  0.27  exp(E‘  -  E  * ) _ - 

V  F  n 

where  exp(cqV')  *  [exp(qV')]  and  (55)  have  been  used. 

L  L 

Figure  8  presents  the  predictions  of  (57)  for  the  behavior  of  the  denomi¬ 
nator  at  several  values  of  c  with  F  =  0.1  and  for  c  *  I  with  F  =  0.05.  When 
the  number  calculated  from  (57)  is  greater  than  unity,  the  excess  density  of 
minority  carriers  (holes)  is  less  than  the  density  predicted  by  the  non-degenerate 
equation  (50).  This  deviation  does  occur  for  higher  concentrations  at  applied 
voltages  near  the  low  injection  limit  of  voltage.  Equation  (^9)  does  not  reduce 
to  (50)  in  all  cases  but  does  in  the  proper  limits. 

The  discussion  of  the  voltage  dependence  of  the  denominator  of  (15),  which 
we  postponed  in  Section  II,  proceeds  similar  to  the  above  discussion  for  (i»9)  . 

Our  purpose  is  to  show  that  the  denominator  does  not  depend  on  voltage  for  volt¬ 
ages  below  approximately  97?  of  the  low  injection  limit  and,  therefore,  that 
Figure  I  is  valid  for  voltages  (including  zero)  less  than  this.  Figure  9  illus¬ 
trates  the  departure  of  the  denominator  from  its  equilibrium  value  for  voltages 
equal  to  0.97  qV 1  and  1.0  qV^.  The  curves  substar1 iafe  the  remark  in  Section  II 
that  Figure  I  is  valid  for  voltages  less  than  approximately  97?  of  the  low 
injection  voltage  limit.  The  calculations  for  Figure  9  were  made  with  F  ■  0.1 
which  is  probably  the  greatest  amount  one  would  define  as  low  injection.  Smaller 

values  of  F  would  lead  to  smaller  departures  from  the  curve  of  Figure  1  for  any 
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given  voltage  including  q\T.  At  •  3  x  10  cm  for  example,  with  F  ■  0.1 
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p 

4)  exp(EfEi)" 

1 

/ 

exp(E'-E')  ♦  0.27 
v.  V  F  n 

mc  1 

1  -  F  T, 
m..  1 

3/2' 

1 

i 

cxp(-E^) j 

(57) 


and  qV*  -  qV^  ,  Figure  9  shows  that  the  ratio  is  0.943  whereas  a  calculate 

would  show  that  the  ratio  is  0.991  If  F  -  0.01  and  qV 1  =  qV' . 

l' 

6.  THE  CURRENT  DENSITY  EQUATION 


The  electrical  current  density  due  to  electrons  is  qi 


given  by 


Jn  =  n  lx  +  9  D  22 


in  one  dimension.  The  n  and  —  fartnrc  I  c\  . 

n  ana  dx  Tactors  in  (5)  may  be  evaluated  by  using  the 


non-equilibrium  form  of  {(■>), 


"  ■  "r  <t;«  -  t;»  . 

0 

d. 

where  the  Fermi-Dirac  integral  is  analogous  to  (8). 
The  derivative  in  (57)  may  be  written 


dn  I  3_n_ 

dx  I  ,  E ' 


dEc  in 

57  +  se: 


ond  it  is  noted  that  [4] 


dE'  qc 

C  _  __x 

dx  ~  kT 


The  two  partial  derivat 


ves  are  evaluated  using  Liebnitz1  rule  [10]  to  obtai 


Nc  -  i  <E;e  -  £;> 
2 


a 

■.  1 


V  1/2 

E')  df(E') 
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The  current  density  equation  becomes 


Jn  *  qiJnncx  +  q  0 


-  ['  N‘  \ 


«;.-*£)  5T-  k  Nc  \  (*•-«;)  /2^r>  dE- 


C  dE1 


f  3_ 

Jn  *  qUnnex  -  ,0  i  I»c  *l/2  <E^-E£)  ]  U 


%r  'x  ♦ 


i  N  \ 

/T  c  J 


df(E') 


If  Li  and  Lindholm's  generalized  Einstein  relationship  [ll] 


^-l/2(Ep  -  E^)  kT 

"il/2(E*  -  E‘)  "  q  °n 


is  assumed  to  be  valid  in  non-equilibrium,  equation  S65)  can  be  rewritten 


n  dx 


-  la  N  0  \  ( E  •  -  E ' 

x  rr  c  n  3  c 


1/2  df(E') 


dE',  (67) 


where  (59)  has  been  used  also.  If  (67)  ' s  integrated  by  parts  and  combined  with 
(59)  and  (66),  the  result  is 


J  =  u  n 

n  n 


This  equation  is  identical  to  that  derived  for  the  non-degenerate  case  [4] 
but  it  has  been  derived  by  considering  Fermi-Dirac  statistics  and  the  variation 
of  the  Einstein  relation  with  impurity  concen  rat  ion ;  therefore,  (68)  has  been 
shown  more  generally  to  be  valid  at  higher  impurity  concentrations  than  previously 
demonstrated  with  the  conventional  non-degenerate  Einstein  relation. 


The  corresponding  equation  for  the  current  due  to  holes  is  given  by 


The  thought  may  occur  to  the  reader  that  one  of  the  basic  assumptions  of  this 
work  has  been  the  constancy  of  the  quasi-Fermi  level,  but  (68)  includes  the 
spatial  derivative  of  the  quasi-Fermi  level.  The  rationale  for  this  paradox 
lies  in  the  meaning  of  the  assumption.  The  assumption  that  was  really  made  was 
that  the  change  in  the  quasi-Fermi  level  across  the  junction  was  so  small  that 
it  was  negligible  for  the  prupose  of  computing  the  voltage  across  the  junction. 
However,  the  small  change  that  is  present  is  represented  by  an  equally  small 
slope  which  is  sufficient,  when  multiplied,  cf.(68),  by  the  enormous  number 
of  carriers  present,  to  account  for  current. 

7.  DISCUSSION  AND  CONCLUSIONS 

We  set  out  to  develop,  under  certain  assumptions,  a  set  of  equations  for 
pn  junction  behavior  which  would  be  valid  at  impurity  concentrations  beyond 
the  conventional  non-degenerate  range.  The  approach  we  took  was  to  remove 
the  assumption  of  non-degeneracy  and  to  make  an  approximation  (9)  for  the 
Fermi-Dirac  statistics.  This  approximation  is  tractable  for  engineering  calcu¬ 
lations  and  is  accurate  to  +  1%  if  the  Fermi  level  is  not  more  than  1.3  kT  into 
the  allowed  band.  The  resulting  equations  for  step  junctions  reduce  to,  or 
are  identical  with,  the  conventional  equations  in  the  limit  of  low  impurity 
concentration. 

The  major  assumption  which  served  as  the  basis  for  the  work  was  that  the 
semiconductor  band  edges  and  the  impurity  energy  levels  were  unaffected  by 
increasing  impurity  concentrations.  Experimental  and  theoretical  results  show 
that  this  assumption  is  not  valid  [12,13*1^*3-  The  single  energy  level  of  the 
impurity  atoms  at  low  concentration  widens  to  a  band  of  energies  due  to  the 
impurities  when  the  concentration  is  great  enough  that  the  impurity  atoms  inter¬ 
act.  The  onset  of  the  widening  has  been  reported  to  be  as  low  as  at  approxi¬ 
mately  lo'7  donor  impurities  per  cubic  centimeter  in  silicon  [12],  The 


widening  of  the  impurity  level  to  a  band  changes  the  effective  activation 
ene  gy  of  the  impurity  and  affects  the  carrier  density  accordingly.  Further¬ 
more,  the  density  of  states  at  the  band  edges  deviates  from  the  parabolic 
dependence  on  energy  which  we  have  assumed  {k)  and  states  become  available 
in  the  previously  forbidden  energy  gap.  This  deviation,  often  called  "band 
tails",  also  increases  with  increasing  impurity  concentrations  and  further 
affects  the  activation  energy  of  the  impurities.  In  fact,  the  impurity  band 
and  the  allowed  band  merge  at  high  impurity  concentrations  and  the  activation 

energy  becomes  zero.  Bonch-Bruyevich  [1]  reports  that  this  occurs  at  donor 
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impurity  concentration  of  approximately  10  cm  in  silicon  and  approximately 
3  x  lo'^cm  ^  in  germanium.  The  effect  of  impurity  bands  and  bard  tails  on 
the  results  derived  in  this  interim  report  are  being  investigated. 

In  addition  to  the  major  assumption  concerning  the  quantum  density  of 
states,  we  have  tacitly  made  two  other  assumptions.  The  validity  of  these 
two  merits  discussion.  First,  the  existence  of  the  tunnelling  phenomena  is 
well-known  and  its  occurrence  in  semiconductor  device  junctions  is  the  basis 
of  the  Esaki  diode.  The  development  reported  in  this  report  has  proceeded  on 
the  assumption  that  no  tunnelling  occurs  across  the  junction.  The  fact  is 
that  the  importance  of  tunnelling  increases  as  the  impurity  concentration 
increases  into  the  upper  part  of  the  range  with  which  we  purport  to  deal.  The 
effects  of  tunnelling  on  the  present  results  are  not  significant  for  diffused 
silicon  junctions  [15],  but  should  be  considered  for  germanium  junctions. 

Second,  the  assumption  that  the  quasi-Fermi  levels  are  constant  in  the 
junction  region  allowed  us  to  relate  the  splitting  of  the  quasi-Fermi  levels 
to  the  applied  voltage.  Sah's  development  [8]  dealt  with  the  special  cases 
of  a  symmetric  step  junction  and  a  linearly  graded  junction.  The  results  for 
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forward  biased  junctions  of  either  profile  are  equivalent  in  that  the  quasi- 
Fermi  levels  are  predicted  to  be  approximately  constant  in  the  space  charge 
region.  Sah's  development,  however,  depends  on  the  assumption  of  non-degen¬ 
eracy,  although  he  has  suggested  that  the  results  are  valid  in  general  [15]. 

Despite  the  limitations  imposed  by  the  underlying  assumptions  concerning 
the  density  of  states,  the  new  equations  developed  in  the  interim  report  pro¬ 
vide,  in  some  practical  instances  of  pn  junction  design,  information  superior 
to  that  supplied  by  conventional  expressions.  The  usefulness  of  the  new  equa¬ 
tions  can  be  appreciated  through  the  following  reasoning. 

Both  the  conventional  and  the  new  equations  depend  on  these  assumptions: 
negligible  tunneling,  negligible  impurity  concentration  effects  or.  the  densities 
of  states,  and  negligible  variation  of  the  quasi-Fermi  levels  in  junction  regions. 
The  new  equations,  reviewed  below,  do  not  depend  on  the  assumption  of  non-degen¬ 
eracy  whereas  the  convernional  equations  do.  But  the  conventional  equations  are 
often  used  to  make  calculations  in  which  the  non-degeneracy  assumption  is  vio¬ 
lated.  Hence  the  new  equations  are  clearly  more  valid  in  these  situations. 

Equation  (15),  with  (17)  or  (18),  permits  the  calculation  of  the  product  of 
the  number  densities  of  charge  carriers  at  any  point  in  a  semiconductor  in  equili¬ 
brium  or  with  applied  voltage.  The  evaluation  is  particularly  easy  at  the  edges 
of  the  space  charge  region  of  a  pn  junction.  Also,  if  the  depletion  approximation 
is  used  to  find  the  spatial  variation  of  the  edges  of  the  allowed  bands  through¬ 
out  the  space  charge  region,  Equation  (15)  allows  the  calculation  of  the  np  prod¬ 
uct  in  this  region  for  use  in  recombination  calculations.  Figure  1  illustrates 
that  the  np  product  differs  s i gn i f i cant l v  from  the  conventional  prediction  at 
high  impurity  concentrations  for  voltages  between  zero  and  approximately  97% 
of  the  low  injection  voltage  limit. 

The  position  of  the  Fermi-level  relative  to  the  band  edges  may  be  calculated 
from  (27)  or  (30)  in  n-type  material  and  (28)  or  (31)  in  p-type  material.  Equa¬ 
tions  (27)  and  (28)  indicate  the  dependence  of  the  Fermi  levels  on  the  impurity 
concentration,  the  semiconductor  material  properties,  and  the  impurity  activation 
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energy.  Equations  (30)  and  (31)  provide  an  auxiliary  means  of  calculation. 

Figure  2  illustrates  the  predictions  of  (27)  and  (28)  for  boron  or  phosphorous 
doped  silicon.  The  Fermi-level  approaches  the  band  edge  with  increasing  impurity 
concentration  as  expected. 

The  contact  potential  of  a  pn  junction  is  given  by  (^1)  or  (^M .  Equation 
(M  )  shows  that  correction  terms  must  be  added  to  the  conventional  equation  for 
use  at  high  impurity  concentrations.  Equation  (bb)  requires  that  the  same  energy 
differences  be  known  as  are  needed  it,  (M  )  but  (^*0  provides  for  ease  of  calcu¬ 
lation.  Figure  3  shows  the  behavior  predicted  for  junctions  between  boron  and 
phosphorous  doped  silicon.  Figure  illustrates  the  correction  terms  in  (4 1 ) 
and  thereby  the  errors  inherent  in  using  conventional  equations  ,  (A2)  or  (^3) 

Equations  (^9)  and  (91)  express  low  injection  minority  carrier  densities 
at  the  edge  of  the  space  charge  region.  Figure  8  indicates  that  the  predictions 

from  (1*9)  (and  (51)]  differ  from  those  of  the  non-degenerate  equation  (50)  at 

high  impurity  concentrations  and  at  voltages  near  the  low  injection  limit. 

A  quantitative  definition  of  t  he  voltage  I  i  r'<  i  t  for  low  injection  is  given 
by  (55)  for  injection  into  n-type  material  and  by  (56)  for  injection  into  p-type 
material.  Figure  6  indicates  the  predictions  of  (55)  6y  presenting  qV^  versus 
(Ep  "  with  various  values  of  F  as  a  parameter.  F  is  the  parameter  that 

defines  the  limit  of  low  injection.  The  figure  illustrates  the  expected  result 
that  the  more  injection  one  is  willing  to  allow  the  more  voltage  can  be  applied. 

Finally,  the  equations  for  current  density  due  to  electrons  and  holes  are 

presented  in  (68)  and  (691.  These  equations  have  the  same  form  as  the  non-degen¬ 
erate  equations  but  they  have  been  derived  by  accounting  for  the  dependence  of 
the  Einstein  relation  on  impurity  concentration. 

8.  COHMEIIT 

It  should  be  noted  that,  wh i I e  the  discussion  has  been  cast  in  terms  of 
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APPENDIX  A 


Figure  10  provides  a  reference  diagram  of  the  energy  bands  of  a  semi¬ 
conductor  in  the  region  near  a  pn  junction.  The  list  of  symbols  defines  those 
used  in  the  paper  and  the  number  in  parenthesis  indicates  the  equation  in  which 
the  symbol  is  first  found.  Throughout  the  paper  a  superscript  "prime"  on 

energies  indicates  that  the  energies  have  been  normalized  by  kT ,  i , 

..  E 

t  3  .  Subscripts  n  or  p  denote  quantities  dealing  with  electrons  or  holes. 

The  same  subscripts,  n  and  p,  arc  used  to  denote  quantities  to  be  calculated 

on  the  n-  or  p-side  of  the  junction.  For  example,  (E‘  -  E'  )  means  that  the 

F  C  n 

energy  difference  between  the  Fermi-level  and  the  conduction  band  edge  is  to 
be  calculated  on  the  n-type  side  of  the  junction. 


_  1  Q1. 


SYMBOLS 


-1/2 

g(E) 


defined  in  (B-l) 

a  fraction  applied  to  the  low  injection  voltage  limit  (57) 
carrier  diffusivity  (58) 
general  energy  (I) 

the  energy  of  the  acceptor  level  in  the  forbidden  gap  (28) 
the  energy  of  the  conduction  band  edge  (3) 
the  energy  of  the  donor  level  in  the  forbidden  gap  (22) 
the  Fermi-level  (I) 

the  quasi-Fermi  level  for  electrons  (12) 
tht  quasi-Fermi  level  for  holes  (13) 
the  forbidden  gap  width  =  E^-Ey  0*0 
the  energy  of  the  valence  band  edge  (5) 
the  x-component  of  electric  field  (58) 
defined  in  (52) 

the  distribution  function  for  electrons,  the  probability  of  an  electron 
occupying  a  state  at  the  energy  E  (I) 

the  Fermi-Dirac  integral  of  order  1/2  (6) 

the  Fermi-Dirac  integral  of  order  -1/2  (62) 

the  quantum  density  of  states  at  the  energy  E  (M 

Planck's  constant  (3) 

electrical  current  density  (58) 

defined  in  (15) 

Boltzmann's  constant  (I) 

the  dens  i  ty-of-s tates  effective  mass  of  electrons  in  the  conduction  band 

the  dens  i  ty-of-states  effective  mass  of  holes  (5) 

the  number  density  of  electrons  in  the  conduction  band  (12) 

the  excess  number  density  of  electrons  at  the  edge  of  the  space  charge 
region  in  p-type  material  *  np  “  npo  (51) 


nl 


rc 


u 


the  intrinsic  number  density  of  electrons  (15) 

the  number  density  of  acceptor  atoms  (28) 

the  effective  density  of  states  in  t fie  conduction  band  (6) 

the  number  density  of  donor  atoms  (22) 

the  number  density  of  ioni/ed  donor  atoms  (21) 

the  effective  density  of  states  in  the  valence  band  (II) 

the  number  density  of  holes  ( I  '3 1 

the  excess  number  tensity  of  holes  at  the  edge  of  the  space  charge 
region  in  n-type  material  -  p^^  (^9) 

the  magnitude  o<  the  electronic  charge  (14) 

the  temperature  in  °K  (I) 


genera  I  vo I tage  ( I  k) 

the  contact  potential  of  a  pn  junction  (<tl) 

the  low  injection  voltage  li  it  (5M 

def  i  ned  in  (2M 

def  i ned  in  (2M 

def  i  ned  in  (2M 

carrier  drift  nobility  (58) 
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APPENDIX  B 


CONSTANTS  USED  IN  SAMPLE  CALCULATIONS 
(See  Appendix  A  for  Definition  of  Symbols) 

T  -  300°K 

q/kT  *  38.685  ev 
mc  =  I  •  08  ( 6 ) 

mv  =  0-59  m0  (6] 

Nc  =  2.816  x  IO,3cn'3 

Nv  =  1 . 1 37  X  I  0 1 9Cm" 3 

Eg*  A3. 3  kT  units  (  1.12  ev)  (6,9) 

Eq  -  Eq  -  1.7^08  kT  units  p  in  Si  (3) 

E^  -  Ey  =  I.7A08  kT  units  B  in  Si  [3] 


FIGURE  CAPTIONS 


Figure  I 

Figure  2 
Figure  3 
Figure  k 

Figure  5 

Figure  6 

Figure  7 

Figure  8 

Figure  9 


The  behavior  of  the  np  product  at  the  edge  of  the  space  charge 

region  of  a  step  junction  in  silicon  at  300°K  as  a  function  of 

the  relative  Fermi-level  (impurity  concentration)  depicted  as  the 

ratio  of  equation  (15)  to  the  conventional  expression,  np  =  nj  exp(qV'). 

The  development  in  Section  V  shows  that  this  curve  is  applicable  to 

cases  where  the  applied  voltage  is  less  than  approximately  97%  of 

the  voltage  limit  of  low  injection. 

The  behavior  of  the  relative  Fermi-level  as  a  function  of  impurity 
concentration  in  phosphorous  and  boron  doped  silicon  at  300°K. 

.  The  contact  potential  of  a  step  junction  (phosphorous-boron)  in 
s i I i con  at  301°K. 

.  The  magnitude  o'  the  correction  terms  for  the  contact  potential  of 
a  step  junction  cs  a  function  of  phosphorous  and  boron  impurity 
concentration  in  silicon  at  300°K.  The  "100%  ionization"  curve 
gives  the  amount  by  which  the  non-degenerate  equation,  qVp  =  IntNpN^/n?] 
overestimates  the  contact  potential  of  a  symmetric  step  junction.  1 

The  remaining  curves  give  the  amount  by  which  the  non-degenerate 
equation  qV^  =  I n f nnoPpc/n? ]  underestimates  the  contact  potential. 

.  The  carrier  density  versus  impurity  concentration  for  phosphorous 
and  boron  doped  silicon  at  300°K.  The  actual  carrier  densities 
are  contrasted  with  the  prediction  based  on  the  assumption  of 
100%  ionization. 

.  The  voltage  limit  for  low  injection  (holes  into  n-type  silicon  at 
300°K)  across  a  step  junction  versus  the  relative  Fermi-level. 

The  parameter  F  defines  the  low  injection  limit  by  pn  =  Fnno. 

.  The  ratio  of  the  voltage  limit  for  low  injection  (of  holes  into 
n-type  silicon)  to  the  junction  contact  potential  versus  impurity 
concentration  at  300°K. 

.  The  behavior  of  the  minority  carrier  density  (holes)  at  the  edge  of 
the  space  charge  region  of  a  step  junction  in  silicon  at  300°K  as 
a  function  of  the  relative  Fermi-level  (impurity  concentration) 
depicted  as  the  ratio  of  equation  (^9)  to  the  conventional  equation 
P|  =  pnQ  [exp(qV')  -  I].  The  parameter  c  is  defined  by  qV  =  cqV^. 

The  solid  curves  are  for  F  =  0.1  whereas  the  dotted  curve  is  for 
F  =  0.05. 

The  ratio  of  the  np  product  at  certain  applied  voltages  to  the 
equilibrium  np  product  versus  donor  concentration.  These  curves 
show  that  the  np  product  at  the  n-type  edge  of  the  space  charge 
region  in  silicon  does  not  depend  significantly  on  the  applied 
voltages  up  to  approximately  0.97  ■ 
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C.  GENERATION  AND  RECOMBINATION  NOISE  IN  SOLIDS  (J.  Salz) 

1 .  INTRODUCTION 

The  noise  behavior  of  solids  has  been  extensively  studied  theoreti¬ 
cally  as  well  as  experimentally.  An  excellent  treatise  of  this  subject  and 
an  extensive  bibliograph  can  be  found  1\  reference  (1).  As  may  be  expected 
the  theoretical  approach  to  this  problem  must  be  statistical  since  noise  by 
definition  is  a  stochastic  process.  In  this  report  we  treat  generation- 
recombination  noise  from  a  fresh  mathematical  point  of  view.  Here  we  re¬ 
strict  the  analysis  to  homogeneous  materials  and  consequently  our  descrip¬ 
tive  statistical  parameters  will  be  independent  of  spatial  coordinates.  We 
focus  attention  on  the  most  important  parameters  of  the  noise  such  as  the 
mean  value,  mean-square,  covariance  and  spectral  density. 

As  discussed  in  reference  (1)  the  characteristics  of  the  noise  in  homo¬ 
geneous  solids  can  be  completely  determined  by  the  statistical  behavior  of  a 
vector  stochastic  process  V(t)  *  (v^t),  V2(t)  .  .  .  vN(t)  whose  components 
represent  the  number  of  carriers  at  time  t  in  energy  level  Ej,  j  *  1,  2  .  .  .  N 
respectively  subject  to  the  constraint  that  Zv^  is  fixed. 

In  previous  treatments  of  this  problem  it  was  assumed  that  the  vector 
stochastic  process  V(t)  is  Markovian.  This  assumption  opens  up  a  whole  body 
of  mathematical  works  devoted  to  the  study  of  Markov  processes.  In  particular 
the  theory  of  Fokker  and  Planck  (2),  (3)  is  most  applicable  and  the  derivations 
of  the  noise  statistics  under  these  assumptions  reduce  to  the  calculation  of 
various  conditional  moments. 

Our  approach  is  more  direct.  We  attack  the  problem  by  evaluating  par¬ 
tial  derivatives  of  the  characteristic  function  for  the  vector  stochastic 
process  from  which  we  extract  systems  of  differential  equations  in  the  de¬ 
sired  statistical  parameters.  Our  results  then,  are  solutions  to  linear 
state-equations.  It  appears  that  our  approach  has  greater  applicability  and 
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is  not  restricted  to  only  Markov  processes.  It  has  however  been  observed  (4) 
that  even  the  Fokker-Planck-Kolmogorov  equations  appear  to  be  of  wider  appli¬ 
cability  than  just  to  Markov  processes. 

In  Section  I  we  define  the  problem  and  describe  the  mathematical  model. 

In  Section  II  we  derive  the  covariance  and  the  spectral  density  matrices 
while  in  Section  III  we  apply  the  general  results  to  a  specific  situation. 

2.  THE  MATHEMATICAL  MODEL 

Our  basic  objective  in  this  report  is  to  derive  formulas  for  the  perti¬ 
nent  statistics  of  generation-recombination  noise  arising  in  homogeneous 
solids.  Our  aim  then,  is  first  to  formulate  a  phenomenological  model  based 
on  fundamental  physical  concepts  and  then  to  analyze  this  model.  A  most 
useful  statistical  function  that  can  be  empirically  determined  is  the  spectral 
density  of  the  noise.  This  function  describes  the  distribution  of  the  average 
power  as  a  function  of  frequency.  The  area  under  this  curve  gives  the  total 
average  noise  power  —  a  most  fundamental  parameter  of  devices.  Another  im¬ 
portant  statistical  function  associated  with  noise  is  the  correlation  func¬ 
tion  which,  as  will  be  seen,  is  related  to  the  spectral  density  by  Fourier 
transformation.  We  shall  be  principally  interested  in  deriving  formulas  for 
these  functions  in  terms  of  fundamental  properties  of  solids. 

Our  working  physical  model  is  the  following:  A  dc  voltage  is  applied 
to  a  sample  of  homogeneous -so lid.  As  a  result  a  current  IQ(t)  can  be  observed. 
We  presume  this  current  to  be  modeled  as  a  stochastic  process.  The  random 
aspects  arise  from  carrier  fluctuations.  The  current  flowing  in  this  circuit 
is  assumed  to  be  given  by  a  linear  combination  of  more  elementary  counting 
processes  as: 

N 

I0(t>  ■  £  «n*„(t)  (1) 

n  ■  1 
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where  the  a  '  s  are  constants  determined  by  fundamental  parameters  of  the 
n 

material  such  as  mobilities,  carrier  charges  and  applied  potential.  The 

set  of  functions  v  (t),  n  ■  1,  2  .  .  .  N  represent  the  number  of  carriers 
n 

at  time  t  associated  with  corresponding  energy  levels  En>  n  -  1,  2  ...  N. 
These  functions  are  counting  stochastic  processes  whose  differentials  Av^, 
n  ■  1,  2  .  .  .  N  can  take  or  only  ±  1  or  0.  The  statistics  associated  with 
these  Increments  completely  determine  the  behavior  of  I0(t)*  F°r  instance, 
the  average  value  of  IQ(t)  is 

<Io(t)>  -  E  an  <Vn(t)>  (2) 

n 

where  <  .  >  denotes  the  ensemble  average.  The  fluctuating  portion  of  IQ(t) 
is 

I(t)  -  I  (t)  -  <I(t)>  (3) 

o 

and  the  total  average  noise  power  is 

02-  <I2(t)>  (4) 


In  the  next  section  we  shall  evaluate  these  averages  more  explicitly 
The  spectral  density  of  I(t),  or  the  spectrum  is  defined  (5)  as 


FW  -  £  fe  <  |WT(w) 


(5) 


where  W_(“)  * 
T 


I 


-iwt 

I(t)  e  dt.  In  terms  of  the  original  process  (2) 


-T 


<|WT  ( | > 2  -  E  anam  <An(w»t)  A  B(a).t)> 


n,m 
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where  A^w.t)  -  f  Un(t)  e"iu!t  dt  with  Ufl(t)  -  vn(t)  -  <vn(t)> 

-T 

and  *  denotes  the  complex  conjugate.  Thus  we  see  that  the  spectrum  of  I(t) 
is  given  by  a  quadratic  form. 


-  -  -  T 

F(oi)  -  A  G(u>)  A  (6) 

where  the  vector  A  *  (a^  -  -  -  a^)  t  p?  denotes  the  transpose  of  A,  and 

G(<*>)  is  an  NxN  matrix  with  elements  G^fa)  given  by 

p  /  \  m  11®  -i—  <An^a)’'^  Am(w»T)>  (7) 
Gnm(ai)  T-x»  2T  n  111 


Eqn.  (7)  can  also  be  obtained  by  the  Fourier  transform  of  the  corre¬ 
lation  coefficients-*- 


R  (t)  -  <U(t)  U  (t4t)> 
nm  n  m 


(8) 


To  see  how  this  is  done  write  Eqn.  (7)  as 

rT  T 

lim  1 


Gnm(a))  _  T*°°  2T 


-L-  [  f  iu)(t'-t)<jj  (t)  u  (t')>dt  dt» 

2T  J  J  n  ® 


-T  -T 


lira  1 

■p-xo  2T 


T  T 


—  [  [  R 

2T  J  J  nm 


(t-t')  eiu)(t'-t)  dt'  dt 


-T  -T 


lim  1 

p-xao  2T 


p  2T 

w  [  \  (I-*)  R- 


2T 


ia)x 


-J 


(x)  e  dx  -  (x-T)  R  (-x) 

nm  l  11111 


+  We  assume  that  the  stochastic  processes  Un(t) ,  n  -  1  .  .  .  N  are  ergodic 


00 


‘V<-> 


lotx 

e  + 


Rnm(-X> 


iwx 
e  ] 


dx 


(9) 


Using  matrix  notation  we  can  express  (9)  in  a  compact  form. 


G(w 


[R(x)  eiux  +  rT(x)  e  iUX]  dx  > 

t 


(10) 


where  R(x)  is  the  covariance  matrix  with  elements  R^te).  ln  (10)  we  used 
the  fact  that  R(-x)  ■  RT(x).  This  can  be  established  by  noting  that  for 
stationary  processes 

Rnm(x)  -  <Un(o)  Um(x)  > 

-  <U  (-x)  U  (o)> 

n  m 

-  <Um(o)  Un(-x)>  -  Rmn(-x)  (11) 

The  "real"  part  in  (10)  must  be  taken  since  both  G(to)  and  R(x)  are  real 
matrices. 

In  the  next  section  we  shall  derive  in  detail  the  form  of  the  covariance 
and  the  spectral  density  matrices. 
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3.  Tnts  COVARIANCE  AND  THE  SPECTRAL  DENSITY  MATRICES 


The  covariance  matrix  of  a  zero-mean  vector  stochastic  process  is 
defined  as 

R(t1(t2)  =  <U  T(t1)  U(t2)>  <12> 

When  the  process  is  stationary  (12)  becomes 

R(t)  -  <UT(o)U(t)>  (13) 

where  t  *  t-  ~  tj_«  Making  use  of  the  sub-conditioning*  property  of  expec¬ 
tations  we  can  write  (13)  as 

R(t)  -  <  U  T(o)  E  (0(t)  I  U(o)}>  (14) 

where  E  {•(  }  denotes  conditional  expectation. 

We  shall  now  derive  a  differential  equation  for  the  conditional  mean- 
vector  E  { U ( t )  |  U (o) }  through  its  characteristic  function.  In  what 
follows  we  shall  work  with  the  original  vector  process  V(t)  **  (v^(t),  v2(t)  . 
.  vN(t)  defined  in  Eqn.  (1).  Recall  that  U(t)  =  V(t)  -  <V(t)>. 

For  the  V  process  the  conditional  characteristic  function  is  defined  as 

C(o),t)  =  E  (exp  i<i>V(t)  |  Vq  }  (15) 

where  <*>  ■  (w^  .  .  .  <*>  )  and  V  =  V(o).  The  partial  derivative  of  C  with 
respect  to  t  is  interpreted  in  the  following  manner 


+  The  sub-conditioning  property  of  expectations  is  simply  a  statement  of  the 
fact  that  averages  can  be  carried  out  by  first  computing  a  conditional  aver¬ 
age  and  then  averaging  with  respect  to  the  conditioned  variables. 
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3C  .  lim  i  v  {  e1^  V  T(t+6)  -  e*  “  ^(t)  I  V  } 
3t  6~*o  5  ° 


li»  1  E{ei"VT(t)  [ei»4vT(t).i)|?o 


(15) 


where  ^V(t)  ■  V(t+6)  -  V(t)  and,  we  assume  that  the  above  limit  exists. 

By  again  invoking  the  sub-conditioning  property  of  expectations,  we  can 
put  (15)  into  the  form 

X 

1C  _  lim  1  E  {gi  »  V  T(t)  £  {gi  »  AV  (t)_  1  |y  v(t)}  |  v  }  (16) 

3t  6->o  6  o  o 


We  now  focus  attention  on  the  inner-conditional  expectation. 


E  {  e1  u  A"vT(t)  -  1  |  V0,  V(t)  } 


i  u  v  (t) 


-  1]  P(AV  j V  ,  V(t)  )  d  (AV) 


where  P(*{)  is  the  conditional  probability  density  of  AV. 

Since  the  components  of  the  process  V(t)  can  change  only  by  one  unit  at 
any  t,  the  incremental  components  Avj(t),  j  ■  1,  2  .  .  .  N  can  assume  only 
the  three  possible  values  1,  -1,  or  0.  It  is  reasonable  to  assume  that  the 
vector  increment  AV(t)  conditioned  on  V(t)  and  V(o)  does  not  depend  on  both 
the  process  at  time  t  and  an  earlier  time  t  ■  o.  This  essentially  presumes 
that  the  process  V(t)  is  Markovian.  As  a  consequence,  therefore. 


P  (AV(t)  |  V0,  V(t)  )  -  P  (AV(t)  |  V(t)  )  (18) 


Ue  now  define  the  following  probabilities 


(19) 


Pr  [AV(t')«  Ak  I  V(t'),  t ' e[t ,  t+fi]  ]  -  gk(V(t))6 

where  Ak  ■  (akj  •  •  •  akN)  and  ak£  is  either  1,  -1  or  o,  k,  j  ■  1,  . . 
A  typical  Ak  vector  might  look  like  (1,  1,  -  1,  0,  0,  1,  .  .  .1).  There 
can  be  at  most  K  -  3n  such  vectors  but  because  of  physical  constrains  the 
permissible  number  of  vectors  may  be  less.  The  scalar  functions  g^(V)  are 
the  probabilities  per  unit  time  that  the  increment  Av  given  V  takes  on 
the  vector  A^  .  Upon  substituting  (19)  into  (17)  we  obtain 


E  ( 


iu)  AV 

e  -  1 


V(t)J 


gk(V)  [e 


ia)  AkT 


-1  1  « 


(20) 


and  the  basic  equation  (16)  becomes 


3C 

at 


[e 


_  _  T 

iui  Ak 


-  1]  S  {  e 


iw  v  T(t)  gk(V) 


V0  } 


(21) 


We  now  derive  the  conditional  mean  of  V  needed  in  the  evaluation  of  the  co- 
variance  matrix. 

Towards  this  and  let 

3C 

m£(t)  =  E  {v£(t)  |  VD}  -  -i  3^  lu-o  ’  *  "  1,2  *  *  *' 

From  Eqn.  (11)  we  obtain 

(t)  r->  -  I  - 

ST  ■  A  VE  <8k<v)  '  vo  > 

k 

-  E  <£  aM  *k(V)  I  V0  } 

k 

-  E  {  f£(V)  I  Vo  }  ,  l  -  1,  2  .  .  .  N 
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where 


f*<V)  “  L  akjlgk(V) 

k 


In  order  to  make  further  progress  we  approximate  f^(V)  by  the  linear 


form 


where 


f£(V)  ~  f n  (N  )  +  E  9-£i(No).  (v  -  N  ) 

Z  °  n-1  3vm  ®  °m 


no  -  <v(t)>  -  (no1,  no2  .  .  .  noN) 


is  the  steady-state  mean  vector.  The  mean  vector  components  Nq^  can  also 
be  obtained  from 


oil  t-**> 


m£(t) ,  £  *  1,2.  •  .  N. 


When  we  substitute  (24)  into  (23)  we  obtain 


~  "  MNJ  +  E  ^ 


E  (  v_  -  N  Vn  } 


3t  ,  3V  1  m  om  '  o 

m«l  m 


By  taking  the  limit  as  t  00  on  both  sides  of  (26)  we  conclude  that 


f  Jl  (Nq)  "  0,  £  -  1,  2  ...  N 


since 


E  (vm  -  Nom  |  v0  }  -  0 


Because  N  ,  n  ■  1,  2  .  .  .  N  does  not  depend  on  time  we  can  write  (26)  in 
on 

the  form 


IHi  .  E 

at 


b*mUm  ,  a  -  1,  2  ...  N 
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3f  (V) 


V  =  N 


A  more  compact  representation  of  (28)  is  possible  by  making  the  following 
observations.  Define  the  conditional  mean  vector  by 


M(t)  =  E  {  U(t)  |  U(o) } 


whose  elements  are  Un(t)>  n  =  1,  2  .  .  .  N  .  Furthermore  define  the  matrix  F 
whose  elements  are  bt  ,  then  we  can  write  (28)  as 


dM(t) 

dt 


-  M(t)  F 


Post  multiplying  (14)  by  F  we  obtain 


R(C)F  -  <V  T(o)  M(t)  F> 


•<vT<°>  Mil > 

dt 


_  —  <  V  T(o)  M(t)> 
dt 


3y  undoing  the  sub-conditions  in  (31)  we  obtain  the  desired  state-differential 


equation  in  R(t)  as 


dRjt 1 

dt 


-  R(t)  F 


with  the  well  known  solution 


R(t)  «=  R(o)  e' 
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'mmmmmmp&rn*- 


urnmn****  '' 


For  a  complete  specification  of  the  covariance  matrix  we  need  R(o) 
the  steady-state  variance  matrix.  This  matrix  can  be  obtained  by  evaluating 
second  derivatives  of  the  characteristic  function  as  follows.  Observe  that 


3  C—  -  E  {v  (t)  v  (t)  IV  }  -  X  ,  ft,  s  -  1,  2  ...  N  (34) 

dui.b)  “  £  S  l  O  XS 

ft  S  0)  *  O 


By  twice  differentiating  Eqn.  (21)  and  evaluating  the  result  at  oi  =  o, 
we  obtain 


(aks  E  {  Vft8k  (V)  I  V) 


+  aw  E  {  vs  gk(V)  )  VQ  } 


+  akft  aks  E  <«k0r>|vo  } 


E  ( v  fs  (V)  |VQ  }  +  E{vs  ft(V)|  VQ  } 


+  ~  akft  aks  E  {8k(V)  |  Vo  >  ,  ft  ,  s  -  1,  2  .  .  .  N 

k 


Making  use  of  (24)  we  find  that 


E  { v . f  (V)  |V  } 
*  s  o 


■  -E  b..(v»  -  Nom)  I  Vo  ) 

m 


E  (v.  h  <V)  |v„  ) 


-E  {vs  L  b«m  <v„  -  "J  I  Vo  ) 
m 
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Since  we  are  only  interested  in  the  steady-state  situation,  the  left- 
hand  side  of  Eqn.  (35)  can  be  set  to  zero.  (A.  cannot  depend  on  time  because 
of  stationerity . )  Also,  the  conditional  averages  on  the  right  side  become 
ordinary  expectations  which  do  not  depend  on  Vq.  Thus  we  obtain  from  (35) 
with  the  use  of  (36)  and  (37)  the  relationship 


we  obtain  from  (38)  the  matrix  equation 

C  -  F  R  (o)  +  R  (o)  F  T  (40) 

where  C  is  an  NxN  matrix  with  elements 

Co  ■  £  a  a  g  (N  )  (41) 

As  m£  ms  6m  o 

m 

and  the  elements’1"  of  R(o)  are  clearly  <  U  U  >  .  1°  obtain  R(o)  explicitly 

m  x» 

we  must  solve  Eqn.  (40).  As  it  turns  out  the  spectral  density  matrix  involves 
the  matrix  C  only  and  therefore  we  shall  not  attempt  to  solve  (40) . 

The  calculation  of  the  spectral  density  matrix  is  now  straightforward. 

+  The  transpose  of  R(o)  is  equal  to  itself. 
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rom  Eqn.  (10) 


oo 


6(w)  “ 


J  [R(t)  +  R 


R  (t)]  cos  ut  dt 


sing  (33)  In  (42)  we  obtain 


-  ,  „  i  f  ' /  v  -Ft  +  iwt 

G(w)  ■  Re  ij  R(o)  e  dt 

o 


f  e"l’Tt  “  iu)t 

*  n 


R(o)  dt  ) 


■V  -V  -][  —  J.  -  . 

Re  {  R(o)  {F  -  io-'I]  +  [F1  +  iu)I]  R(o)> 


where  [•]  denotes  the  inverse  and  I  denotes  the  identity  matrix.  Since 
G(oi)  is  symmetric  we  can  put  (43)  in  the  simple  form 


G(u>)  -  Re  {[  FT  +  ioJl)"1  l  (FT  +  iu)I)  R(o)  +  R(o)  (F  -  iuI)  ]  {F  -  iwij  ±) 


-  Re  { (FT  +  iwl)"1  [FT  R(o)  +  R(o)  F]  (F  -  iiol)' 


-  T  —1  -  -  —1 

Re  { (F1  +  iiol)  C  (F  -  iuI) 


where  C  -  FT  &(o)  +  R(o)  F  is  a  symmetric  matrix  given  in  Eqn.  (40). 

We  now  wish  to  represent  G(w)  in  normal-modes,  i.e.,  in  terms  of  the 
various  corner-frequencies.  This  can  best  be  accomplished  by  diagonalizing  F. 


Now,  we  suppose  that  F  has  real"*"  eigenvalues  =  1,  2 


N  and  let 


F  -  Q  1  D  Q 


+  As  will  be  apparent  later  these  conditions  must  be  satisfied  in  order  for 
the  mathematics  to  make  physical  sense. 
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where  D  is  a  diagonal  matrix  with  elements  A  ,  £  =  1,  2  ...  N.  Using 
this  representation  for  F  puts  44  in  the  form 

~_1  -1  -  -1 

G(uj)  =  Re  {Q  A  (D  +  iwl)  W  (D  -  iool)  }Q 

(46) 

=  7  Q"1  {L(ai)  w  L(-u)  +  L(-u)  W  L  (oi)  }  Q~T 


where 


'  '  T  -1 

W  =  Q  C  Q  ,  L(co)  =  (D  +  iu>I) 


and  Q  1  is  the  inverse-transpose  matrix.  Since  W  is  a  symmetric  matrix 
we  find  that  the  elements  of  the  matrix 


fl  =  2  L(<^)  W  L. ( — co )  +  2  L ( — co )  W  l (io) 


can  be  represented  by 


Wij(XiXj  +  > 

(A.^2  +  us2)  (Aj2  +  w2) 


Xi  +  Xs  L  *i“  +  w2 


*  ..»j  1 

X  j  2  +  a/2  J 


and  W^'s  are  the  elements  of  W.  Now  if  we  let  the  elements  of  Q-T 
be  cliJT»  the  elements,  a  of  the  matrix  H  Q-T  are 

“U  '  ^  hik  qkJ 


and  therefore  the  elements  of  the  spectral  density  matrix  in  normal  modes 


become 


G  (<*>)  =  Y. 


'‘ns  ^  s” 


^  q; 


s  k 


nk  ^sk 


I  hnlf  r 


nk  km 
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In  the  forthcoming  section  we  shall  apply  these  results  to  a  special 

case. 

A.  NOISE  STATISTICS  IN  A  TWO-LEVEL  SEMICONDUCTOR 
Noise  in  a  semiconductor  may  be  regarded  as  arising  from  carrier  tran¬ 
sitions  that  occur  between  two  energy  levels,  e.g.,  the  conduction  band  and 
a  set  of  impurity  levels.  Because  the  number  of  charge  carriers  is  conserved 
we  need  only  consider  a  single  counting  stochastic  process  v(t). 

Fig.  1  shows  a  schematic  representation  of  the  given  situation. 


*7 


Fig.  1  Two-level  system 
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In  this  application  the  vectors  A^  in  (19)  are  one-dimensional, 
taking  on  the  values  1,  -1  and  0.  Their  respective  probabilities  are 

Pr[Av(t')  -  1  |  v(t');  t'e  [t,  t  +  6]  J  =  r(v)  6 
Pr[Av(t’)  -  1  |v(t’);  t'  [t,  t  +  <5]  ]  =  g (v )  6  (53) 

Pr[Av(t')  =  o  | v(t ' ) ;  t'e  [ t ,  t  +  6]  ]  =  (1  -  r(v)  -  g(v))6 

From  Eqn.  (23)  and  (24)  with  Z  =  1,  N  =  1  we  obtain 


' 

% 


fx(v) 


E  r(v)  -  g(r) 

'  r(No)  -  g(No)  + 


(54) 


dg(v) 

dv 


(v-Nq) 

v=N„ 


where  N  =  <v(t)>  . 

o 

An  important  result  from  (27)  is  that 

r(N0)  »  g(N0)  (55) 


from  which  NQ  may  be  determined. 

Observe  that  in  this  application,  F  has  only  one  element  and  is 

therefore  a  scalar  equal  to 


9fx(v) 


3V 


N 


8'<»o>  -  r’<V 


(56) 


where  the  prime  denotes  differentiation, 
simply 

R(t>  - 


The  covariance  function  now  is 


<v(o)  v(t)> 
R(o)  e"bltl 


\57) 
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Solving  Eqn.  (40)  for  R(o)  now  yields 


r'(N0)  -  g'<N  ) 


The  Fourier  transfer  of  (57),  which  is  the  spectral  density,  is  simply 


G(U>  -  2  8<"o> 


,2  2 

b  -f  oj 
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D.  LOW  FREQUENCY  NOISE  IN  INTEGRATED  CIRCUIT  TRANSISTORS  (E.  R.  Chenette, 

A.  J.  Brodersen,  R.  C.  Jaeger) 

The  low  frequency  noise  performance  of  all  semiconductor  devices  is 
dominated  by  sources  of  excess  noise  with  a  spectral  density  which  varies 
approximately  as  1/f.  Although  much  work  has  been  done  in  the  study  of  this 
1/f  noise  it  remains  an  unsolved  problem^".  It  is  impossible,  at  the  present 
time,  to  predict  the  1/f  noise  performance  of  a  bipolar  transistor  (or  that 
of  any  other  semiconductor  device)  on  the  basis  of  a  knowledge  of  device  and 
circuit  parameters.  At  higher  frequencies  where  the  sources  of  1/f  noise  are 
negligible,  it  is  possible  to  predict  noise  performance  with  precision  for  a 
wide  range  of  operating  conditions. 

Recently  progress  has  been  made  in  relating  the  1/f  noise  observed  in 

o 

MOS  field-effect  transistors  to  the  surface  charge  density  .  Earlier  experi¬ 
mental  studies  have  shown  that  a  major  source  of  1/f  noise  in  otherwise  well- 
behaved  low-noise  bipolar  transistors  is  related  to  fluctuations  in  the  surface 
recombination  velocity  near  the  emitter  base  transition  region^»^»5»6.  It  has 
been  shown  that  controlling  the  surface  potential  near  the  emitter-base  tran¬ 
sition  region  of  a  planar  transistor  can  effect  the  surface  recombination  rate 
and  hence  the  d.c.  current  gain.  It  Is  reasonable  to  expect  that  this  control 
of  the  surface  potential  must  also  effect  the  1/f  noise  since  this  excess 
noise  in  these  otherwise  well-behaved  transistors  Is  proportional  to  the  recom¬ 
bination  current  (base  current). 

The  brief  discussion  which  follows  gives  a  summary  of  some  of  the  prelimi¬ 
nary  work  being  done  to  better  understand  this  problem.  The  results  of  some 
preliminary  measurements  of  noise  in  Integrated  circuit  transistors  will  be 
presented.  Equivalent  circuits  describing  the  noise  performance  of  transistors 
in  the  1/f  noise  region  are  compared. 
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1.  MEASUREMENT  OF  LOW  FREQUENCY  NOISE  IN  INTEGRATED  CIRCUIT 

TRANSISTORS 

Preliminary  measurements  have  been  made  of  the  equivalent  noise  resis¬ 
tance  of  several  types  of  planar  transistors.  The  measurements  are  made  as 
a  function  of  frequency  with  operating  point,  source  resistance  and  environ¬ 
ment  as  parameter.  The  CA  3018  has  been  a  convenient  source  of  typical  inte¬ 
grated  circuit  planar  transistors  for  these  early  measurements.  For  comparison, 
measurements  were  also  made  on  low  noise  2N930  planar  transistors.  Special 
transistors  pairs  with  carefully  matched  geometries  and  including  field  plates 
for  control  of  surface  potential  are  being  fabricated  in  our  laboratory  for 
additional  study. 

Figures  1  and  2  show  the  results  of  some  of  these  measurements.  As  dis¬ 
cussed  below,  these  results  are  in  agreement  with  the  expected  behavior.  At 
this  point  they  serve  mainly  to  build  confidence  in  our  measurement  technique 
and  the  operation  of  the  spectrum  analyzers. 

The  most  significant  feature  noticed  to  this  point  is  considered  to  be 
the  small  spread  in  the  magnitude  of  the  noise  in  this  1/f  noise  region.  This 
is  in  contrast  to  results  reported  earlier  of  large  variations  in  noise  from 
transistor  to  transistor  even  of  the  same  type.  It  is  hoped  that  this  unifor¬ 
mity  is  an  indication  that  a  well-defined  cause  for  the  apparent  modulation  of 
the  surface  recombination  rate  may  soon  be  identified. 

2.  COMPANION  OF  LOW-FREQUENCY  NOISE  EQUIVALENT  CIRCUITS 

Figure  3  shows  three  different  noise  equivalent  circuits  which  have  been 
used  to  describe  the  noise  performance  of  bipolar  transistors  when  flicker 
noise  or  1/f  noise  predominates.  Figures  3A  and  3B  are  basically  common-base 
T-equivalent  circuits.  3A  shows  two  flicker  noise  current  generators;  i^, 
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Noise  Resistance,  Ohms 


FIGURE  1:  Noise  Resistance  of  Typical  Integrated  Circuit 

Bipolar  Transistors 


Frequency,  Hertz 
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Noise  Resistance,  Ohms 


«*3*»s5rp jwpvw  *»»**««*• 


FIGUR.E  2:  Noise  Resistance  of  Low  Noise  2N930  Silicon 

Planar  Transistors 


Frequency,  Hertz 
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in  parallel  with  rg  and  i^  in  parallel  with  rc  and  aIe*  3B  shows  a  flicker 
noise  emf ,  e^  series  with  rg  and  a  current  generator  i^o  in  parallel  with  rc. 
These  two  circuits  are  identical  if: 


Lfo 


=  i 


f  2 


+  ai 


fl 


and 

ef  “  ifl  '  re 
4  6 

Work  by  Chenette  and  Plumb  ’  has  shown  that,  for  a  large  class  of  bi¬ 
polar  transistors  —  with  negligible  flicker  noise  generated  by  collector 
leakage  current  —  the  flicker  generators  of  Figure  3A  or  3B  are  almost  com¬ 
pletely  correlated.  Since  this  is  the  case  the  simplified  representation  of 
Figure  3C  is  often  adequate.  Here  a  single  flicker  noise  generator  is  placed 
in  parallel  with  the  emitter-base  resistance  in  a  hybrid-pi  equivalent  circuit. 
This  circuit  can  be  used  to  predict  an  equivalent  noise  resistance  or  an  ex¬ 
pression  for  noise  figure.  The  results  are: 

2 

R  “  R  4-  g,(Rg  +  r  ) 
n  s  °f  s  x 

or 

F  =  1  +  —  (R„  +  r  )2 

R  s  X 


Here  g^  is  the  equivalent  noise  conductance  of  i^,  i-e* 

if  -  4  kTgf  df 


Empirically  we  know  that 


«f 


Klgr 

fa 


where  the  exponents  a  and  r  are  often  near  unity  —  yielding  the  common  1/f 


dependence  and  a  direct  dependence  on  Ig.  The  results  of  the  preliminary 


measurements  shown  on  Figure  1  and  2  are  in  essential  agreement  with  this 
model.  Some  evidence  of  an  anomaly  at  high  resistances  has  been  observed 
but  more  careful  investigation  is  required. 

Work  now  in  progress  includes  a  more  detailed  study  of  the  effect  of 
source  resistance  and  an  investigation  of  the  possible  effect  of  thermal 
environment  on  the  1/f  noise  source  as  well  as  the  study  of  the  matched  pairs 
with  field  plates  described  earlier.  The  results  of  that  work  will  be  pre¬ 
sented  in  a  later  report. 
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E.  ADDITION  TO  LABORATORY  FACILITIES:  MEASUREMENT  OF  ELECTRONIC  AND 
QUANTUM  PARAMETERS  (S.  S.  Li) 

This  section  outlines  the  progress  made  toward  the  measurement  of  the 
electronic  and  quantum  parameters  of  materials. 

The  development  and  set-up  of  facilities  for  electrical  and  quantum  parameters 
study: 

During  this  period,  a  system  for  measuring  of  Hall-effect,  conductivity, 
magneto-resistance,  thermoelectric,  thermomagnetic  and  PEM  effect  at  tempera¬ 
tures  between  4.2°K  and  room  temperature  has  been  designed  and  built.  The 
basic  equipment  in  this  set-up  consists  of  a  6"  Varian  electromagnet.  Air 
Products  Helium  Cryotip  refrigeration  system,  high  sensitivity  nano-volt 
potentiometer,  sampling  oscilloscope,  nanosecond  light  source  and  chart  re¬ 
corder.  This  set-up  has  been  successfully  tested  using  gold  diffused  silicon 
specimens.  Several  measurements  on  silicon  specimens  are  being  undertaken 
using  these  facilities.  A  liquid  nitrogen  cryostat  has  also  been  built  for 
minority  carrier  lifetime  experiments  using  the  photoconductivity  decay  method. 
Preparation  of  Samples: 

Another  important  aspect  of  progress  is  the  preparation  of  samples.  A 
diffusion  furnace  has  been  adapted  for  the  diffusion  of  gold  in  a  neutral  or 
reducing  atmosphere;  this  equipment  has  been  successfully  tested  by  diffusing 
gold  at  1200°C  into  phosphorous-doped  silicon.  Personnel)  have  been  trained  in 
the  preparation  of  the  surfaces  of  the  samples;  this  preparation  includes 
mechanical  polishing,  lapping  and  chemical  etching  of  the  specimens.  Experi¬ 
mentation  has  yielded  a  reliable  technique  for  making  ohmic  contacts;  the 
technique  makes  use  of  electroless  nickel  plating  and  soldering. 

All  these  preliminary  tasks  have  been  completed  and  samples  are  now 
ready  for  measurements. 
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V.  Discussion 


The  research  of  this  second  semi-annual  period  has  enlarged  our  compe¬ 
tence  in  solid-state  materials  and  devices.  It  has  yielded  results  of 
technical  significance,  and  it  has  laid  the  foundations  for  future  study 
and  growth. 

From  the  research  concerning  glasses  has  come  a  fundamental  understanding 
of  the  role  that  structural  disorder  plays  in  determining  electronic  proper¬ 
ties.  This  knowledge  will  undergird  much  of  the  research  to  come  —  certainly 
in  glasses,  perhaps  in  more  ordered  materials  as  well. 

From  the  research  concerning  measurement  techniques  has  emerged  capa¬ 
bilities  that  will  prove  invaluable  in  many  of  our  envisaged  materi ils 
studies.  The  x-ray  research  has  yielded  a  method  for  determining  particle 
size  distributions;  one  future  effort,  for  example,  will  be  the  adaptation 
of  this  method  to  semiconductors,  enabling  the  measurement  of  impurity-atom 
clustering  in  degenerate  silicon  and  germanium.  Analysis  concerning  the 
electron  microprobe  has  opened  the  way  to  the  employment  of  this  tool  in 
determining  the  structure  of  thin  films;  this  work  has  begun,  utilizing  samples 
of  permalloy  films  made  during  this  period  in  our  magnetic-film  laboratory. 
Finally,  the  analysis  of  thickness  error  in  quantitative  transmission  micro¬ 
scopy  has  provided  another  technique  for  the  examination  of  the  structure  of 
ceramics . 

From  the  research  concerning  semiconductors  and  semiconductor  devices  has 
come  progress  collateral  to  that  reported  with  respect  to  structural  analysis. 
This  period  has  seen  the  acquisition  and  installation  of  our  facilities  for 
Hall-effect,  lifetime,  and  conductivity  measurements  at  cryogenic  and  ele¬ 
vated  temperatures.  Samples  have  been  prepared  and  the  first  measurements 
have  been  made.  The  data  forthcoming  from  these  facilities  will  be  correlated 
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with  that  yielded  from  structural  analysis  and  from  measurements  of  noise 
spectra.  Device-oriented  research  will  continue,  building  upon  the  pre¬ 
liminary  knowledge  reported  here  dealing  with  noise  performance  and  with  the 
influence  of  degenerate  concentrations  of  impurity  atoms. 

Radiation  studies  have  continued  and  have  yielded  knowledge  that  will 
be  described  in  the  report  for  the  next  semi-annual  period. 

The  plan  of  research  will  remain  the  same.  In  the  materials  research, 
correlation  will  be  sought  between  electronic  and  structural  parameters, 
which  will  be  further  correlated  with  radiation  effects.  The  materials 
research  will  both  suggest  collateral  device  studies  and  will  be  guided  by 
them. 

The  progress  described  in  this  report  is  in  accord  with  this  plan. 
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